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A Mixing Theory tor the Interaction 
Between Dissipative Flows and 
Nearly Isentropic Streams: 


LUIGI CROCCO? anp LESTER LEES: 


Princeton Unwersity 


SUMMARY 


By means of a simplified theoretical *‘model,’”’ the present paper 
treats the general class of flow problems characterized by the 
interaction between a viscous or dissipative flow near the surface 
of a solid body, or in its wake, and an “‘outer’’ nearly isentropic 
stream. For the present, the external flow is taken to be a plane, 
steady, supersonic flow, which makes a small angle with a plane 
surface or plane of symmetry, although the methods used can be 
extended to curved surfaces, to axially symmetric supersonic 
flows, and also to subsonic flows. The internal dissipative flow 
is regarded as quasi-one-dimensional and parallel to the surface 
on the average, with a properly defined mean velocity and mean 
temperature. The nonuniformity of the actual velocity distri 
bution is taken into account only approximately by means of a 
Mix 


ing, or the transport of momentum from outer stream to dissipa- 


relation between mean temperature and mean velocity 


tive flow, is considered to be the fundamental physical process 
determining the pressure rise that can be supported by the flow. 
With the aid of this concept, a large number of flow problems is 
shown to be basically similar, such as boundary-layer-shock- 
wave interaction, wake flow behind blunt-based bodies (base 
pressure problem), flow separation in overexpanded supersonic 


nozzles, separation on wings and bodies, et« 
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When the mixing rate is taken to be proportional to the mass 
flux density of the isentropic stream, the equations of motion are 
reduced to a single nonlinear ordinary differential equation that 
can be integrated numerically. An important property of this 
equation is the existence of a “‘critical point”’ for supersonic wake 
flows and, also, under certain conditions, for supersonic flows 
directed toward a solid surface. This critical point acts much 
like the ‘‘throat”’ of a nozzle in determining the base pressure, for 
example, or in some cases the surface pressure distribution in a 
boundary-layer-shock-wave interaction. One important reason 
for the marked difference between laminar and turbulent flows 
is the fact that the turbulent mixing rates are from five to ten 
times larger than the laminar mixing rates. 

By introducing several reasonable physical assumptions, a 
‘simplified’ form of the mixing theory is developed particularly 
for separated and reattaching flows and wake flows. Separating 
flows, as well as reattaching flows, are found to be capable of sup- 
The 
maximum compression, as measured by the isentropic flow de- 
flection, is roughly proportional to the square root of the mixing 
coefficient and also to (.\/? 1)‘, where / is some average 
Mach Number. Thus, the pressure ratio increases rapidly with 
Mach Number, or the base pressure ratio, for example, de- 
creases rapidly with increasing Mach Number. In separated 


porting considerable pressure increases at high velocities. 


flows, the pressure gradient along the surface is a maximum at 
separation and drops off steeply downstream; in reattaching 
flows, or wake flows, the pressure gradient is negligible some dis 
tance upstream of the “reattachment point”’ and increases rapidly 
The inflected 
surface pressure distribution observed in laminar boundary-layer 
shock-wave interactions with flow separation is now understand- 
able, and the distribution itself can be calculated approximately 

When the present mixing theory is applied to the problem of 
determining the base pressure for a supersonic airfoil with a blunt 


toward a maximum as this point is approached 


trailing edge, it gives the correct fluid-mechanical explanation of 
the observed phenomena. Qualitative agreement is found be- 
tween the theoretical calculations of the curve of base pressure 
versus Reynolds Number and the data of Chapman and Bog- 
donoff on bodies of revolution and Chapman’s data on blunt 
trailing-edge airfoils. The theory is now being extended to axially 
symmetric supersonic flows 
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The results obtained in the base pressure problem for a super 
sonic airfoil with a blunt trailing edge open the way for applica 
tion of the mixing theory to boundary-layer-shock-wave interac 
tions, boundary-layer separation, and many other phenomena 
However, the dependence of the mixing rate and the mean 
velocity-mean temperature relation for the dissipative flow region 
on the flow parameters must be carefully investigated experi 
mentally and theoretically and the results incorporated into the 


analysis 


List OF SYMBOLS 


The subscript e denotes quantities in the nearly isentropic flow 
the subscript 1 denotes mean quantities for the dissipative flow; 
the subscript cr signifies values taken at the critical point of the 
basic differential equation. Except where otherwise noted, the 
subscript s denotes stagnation conditions. In addition, the sub 
script 7 applies to quantities of the equivalent incompressiblc 
flow in the Stewartson transformation; the subscript >) denotes 
quantities in the airfoil boundary layer just upstream of the 
trailing edge: the subscript >’ denotes quantities evaluated in 
the mixing zone just aft of the trailing edge, while the subscript 


) applies to the mixing zone in general 


x,y = coordinates parallel and normal to surface, o1 
plane of symmetry 

6 = thickness of dissipative flow regior 

Pp = pressure 

p = density 

r = temperature 

h = enthalpy 

a = speed of sound 

Mu = ordinary coefficient of viscosity 

y = ratio of specific heats, ¢,/« 

R = gas constant 

ue = magnitude of local velocity of nearly isentropic 
stream at y = 6 

VW = Mach Number, u,/a- 

“ = local angle between external stream and x-axis 
atv =6 

v = characteristic angle variable of plane, isentropic 


supersonic flow 


‘“free-stream’’ pressure, Mach Number, re 


vy, etc = 
duced velocity, etc., far from body 
u = component of velocity parallel to surface in in 
ternal flow 
*5 
mt = mass flux, | pu dy 
J0 . 
*5 
I = momentum flux, { pu> dy 
J0 . 
= ‘ r pu 
5* = displacement thickness, 1 dy 
0 7 
Pe lle 
a F Pop uf u ; 
5* = momentum thickness, l dyt 
J9 peu, tu, 
u = mean velocity of internal flow, // 
T; = mean temperature of internal flow, defined by 
m = (pm6/RT;) 
é = t/a 
We = U,/Ag 
m = ma 
Teo 4 | Y 1) /2 ree 
Ler = eed = 
Ts Ye YW, 
¢1 = (71/T,)(1/yvu 
kK = U,/uU, = W/W, 


t The symbol 6** is used in order to avoid confusion with the 


local angle, @ 
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= parameter appearing m mean teniperatu 


mean velocity relation 7)/7 = 


(4 1) /2 jury? 
F = (f/x?) — 1 
k = mixing coefficient, defined by relation dm /da 
Rp, uy 

TH = viscous shear stress at solid surface 

( = local skin-friction coefficient, ry /(1/2)p-4u-* 
if 

Cy = total skin-friction coefficient, ~ cr aX 

o = (f 2k (1 A 

a6. xX.. 

br,, Ox,, Xr = thickness of dissipative tlow, local angle of 


inclination of external stream, and coordi 
nate of separation point (s) or reattachment 


point (R 


\ = X/dp or x/6 
6 Kj(6 6 ri) Aj(6 6 
0 = = or 
6; 6 
6r Or 6 = | K;{1 (o'/0 
Ay, = flow deflection across incident oblique shock 
= airfoil chord 
h = airfoil thickness at trailing edge 
Re, = Reynolds Number, p,i1,c/i, 
Re; = transition Reynolds Number, p.4.x7;/ ue, wher 


‘yr is defined in the text where needed 
Xz = PelteXx Me, Reynolds Number based on location 


of separation point 


L) INTRODUCTION 


O* OF THE LEAST UNDERSTOOD PHENOMENA itl 
fluid mechanics is the interaction between a vis 
cous or dissipative flow near the surface of a solid body, 
or in its wake, and an “‘outer”’ nearly isentropic stream 
with varying static pressure. This phenomenon is ob 
served in all types of boundary-layer-shock-wave inter 
actions, in the wake behind blunt-based bodies and air 
foils with blunt trailing edges at supersonic speeds 
(base pressure problem), in flow separation associated 
with recompression in an overexpanded supersonic 
nozzle, in flow separation from the surface of airfoils 
and bodies at high angles of incidence, ete. Apart 
from its theoretical interest, the problem is of consider 
able technical importance because of the widespread 
occurrence of such interactions in compressors, in 
rocket nozzles, in supersonic and subsonic diffusers, 
and on supersonic wings and bodies. 

In all these interactions, the ‘“‘external flow” cannot 
be regarded as a known datum for the calculation of 
the “internal” dissipative flow. In contrast to the 
usual Prandtl boundary-layer theory and its extension 
to wakes and jets, the development of the dissipative 
flow itself helps to determine the external flow, and it is 
now generally recognized that this interplay must be 
made the basis of any valid theoretical treatment. It 
is the aim of this paper to bring out the importance of 
the transport of momentum from outer stream to dis- 
sipative flow (turbulent “mixing” or laminar diffusion 
in determining the flow pattern and pressure rise in these 
interactions and to formulate this concept in quantita 


tive terms. By means of this concept a large number 
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of seemingly diverse flow problems are shown to be 
basically similar and therefore capable of being treated 
by one general theory. 

Previous theoretical treatments of the interaction 
between an internal dissipative flow and an outer 
nearly isentropic stream do not seem to have recog- 
nized the importance of mixing. For example, in the 
case of the boundary-layer-shock-wave interaction 
problem, Tsien and Finston,' Lighthill,? and others re- 
gard the interaction between the subsonic flow field 
adjacent to the surface and the outer supersonic stream 
as paramount and treat the problem by means of a 
linearized nonviscous flow theory. Their calculations 
show only a negligible upstream spreading of the inter 
action near the surface, and this result corresponds to 
the facts only for a turbulent boundary layer and a 
weak shock. Even in this case the mechanism is not 
correct, because it does not include the mixing between 
external and internal flows and does not take into ac 
count the nonlinear character of the external flow and 
of the interaction. The importance of mixing in de 
termining the pressure rise that can be supported by 
the turbulent boundary layer in an adverse pressure 
gradient has already been recognized by Schubauer 
and Klebanoff* for low-speed flows. 

Treatments of the boundary-layer-shock-wave inter 
action problem which regard viscous effects as impor 
tant, such as those of the junior author,’ utilize the or 
dinary Prandtl boundary-layer theory to determine the 
rate of growth of the displacement thickness, d6* dx, 
in terms of the (unknown) pressure variation in the ex- 
(Here, mixing is taken into account 
approximately.) |The wall is then replaced by a stream 
line of slope dé*/dx, and the pressure variation in the 
external stream is obtained in terms of the curvature 
d*6* ‘dx* by means of the Prandtl-Meyer relation. The 
solution obtained for the laminar boundary layer shows 
an extensive region of upstream influence and an ex- 
ponential decay of the upstream pressure rise along the 
surface away from the shock, and predicts separation of 
the boundary layer ahead of the shock. While these 
predictions now seem to be qualitatively correct, the 
theory is unable to deal with the flow downstream of 
separation and tells us nothing about the subsequent 
reattachment of the flow downstream of the reflected 
expansion fan, where most of the pressure rise along the 
surface apparently occurs. Ritter’ recently attempted 
to improve upon this treatment, but the basic difficulty 
has not been overcome because of the inability of the 
usual boundary-layer theory to deal with this type of 
Since the characteristics of the com- 


ternal stream. 


dissipative flow. 
pressible turbulent boundary layer in an adverse pres 
sure gradient are largely unknown, this technique has 
not been applied to the turbulent layer at all. 

In the case of the base pressure problem, it has long 
been recognized that mixing in the wake just aft of the 
body must determine the amount of recompression 
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that can be supported by the flow and, hence, the base 
pressure, but this concept has never been precisely for 
mulated. Chapman® recently reviewed the entire 
problem and found that existing data on base pressure 
could be correlated in terms of the ratio of the boundary 
layer thickness at the base to the base diameter. No 
mechanism is advanced to explain the correlation. The 
ordinary theory of wakes and jets, which is usually re 
stricted to uniform pressure, is of little assistance in 


this problem. 


Another important phenomenon involving the inter- 
action between a dissipative flow near a surface and an 
“external” flow with increasing pressure in the flow 
direction is observed in an overexpanded supersonic 
nozzle. When the ratio of chamber pressure to re 
ceiver pressure is reduced below a certain critical value, 
it is observed experimentally” * that the oblique shocks 
at the nozzle exit move upstream into the nozzle. The 
stream lines near the nozzle wall are deflected through 
these shocks, and the flow separates from the wall. 
For the conditions investigated, the flow deflection 
angle is remarkably constant (=~ 14°), and the pressure 
ratio across the shock is about 3.0, virtually independ 
ent of nozzle divergence half-angle, in the range from 
5° to 30 The flow deflection and pressure ratio lie 
well below the possible maxima corresponding to the 
local Mach Number ahead of the shock. Nonviscous 
flow theory establishes certain limits for the flow de 
flection angle and shock pressure ratio, but otherwise 
the problem is indeterminate without the introduction 
of viscosity. No theoretical explanation of this phe 
nomenon has as yet been advanced. 


It seems hopeless to attempt to solve the complete 
Navier-Stokes equations in detail, or these equations 
including the Reynolds stresses, for these interactions 
between an internal dissipative flow and an external 
nearly isentropic stream. What is required is a general- 
ization of the concept of the von Karman momentum 
integral for the dissipative flow region, where, however, 
this internal flow is treated as quasi-one-dimensional 
with properly defined mean velocity and mean tempera- 
ture. It is essential to retain the nonlinear character 
of the external flow and of the interaction itself. The 
nonuniformity of the actual velocity distribution across 
the internal flow is taken into account only approxi- 
mately by means of a relation between mean tempera- 
ture and mean velocity. Mixing between the external 
and internal streams is admitted as a fundamental 
physical process, which will furnish the mechanism for 
the pressure rise. While some of the details of the flow 
are undoubtedly lost, it is believed that this simplified 
theoretical ‘‘model’’ preserves the main features of the 
interaction with the external stream. The loss of de- 
tail is outweighed for the present by the generality of 
the results and by the fact that flows are treated which 
cannot be analyzed at all by the classical boundary- 


layer theory 
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(2) GENERAL THEORY 


(2.1) Basic Equations of Motion 


Suppose that the flow along a plane surface or a 
plane of symmetry (x-axis) is characterized by two more 
or less distinct regions (Fig. 1):T Region (1), an inter 
nal dissipative region of local thickness 6, in which the 
nonuniform velocity u is essentially parallel to the x 
axis, on the average; and Region (2), an external isen- 
tropic flow, making the local angle @ with the x-axis at 
For the moment, @ is assumed to be sufficiently 


y = 6. 
Because of the 


small so that tan 6~@ and cos @= 
mass transfer, or ‘“‘mixing,’’ between the external and 
internal streams, dé/dx > 0. 

In order to bring out the main aspects of the theory 
without becoming too involved in mathematical dif 
ficulties, the treatment is restricted for the present to 
plane, steady, isentropic, supersonic external flows, 
for which the Prandtl-Mever relation 


dp/p = (yM2/V M? — 1) de 


holds. Actually, the analysis can be extended to 
axially symmetric supersonic flows by the method of 
characteristics, and shock-wave losses in the external 
flow can also be taken into account. In principle, 
subsonic external flows can also be treated by conformal 
mapping techniques. 

The following assumptions are introduced: 

(i) The static pressure p is constant across the inter 
nal flow and equal to the local pressure /, in the external 
stream at y = 6. 

(11) Heat transfer between the fluid and suriace (if a 
solid surface is present) is zero, and the stagnation en 
thalpy A, is constant across the flow and equal to /i,, at 
y = 6. The applicability of these assumptions will be 
discussed later in brief (see Section 6). 

For the internal stream, the mass flux is given bv 


m= / pu dv 
a] 
the momentum flux is 


7 The treatment can be extended to the more general case of 


a curved surface 


and the flux of enthalpy is 


H = | h,pu dy = mh,g, 


By referring to Fig. 1, one sees that the angle between 
the stream line and the outer boundary of the dissipa- 
tive flow at vy = 6is (dé/dx) — 6; therefore, the rate at 


which mass is transported to the internal flow is 
dm/dx = p,u,|(di/dx) — 6 2.1) 
A momentum balance vields 


dl dm _dp 
— Hu — - 


2) b5:— 6 
dx dx dx 


= Te Ze 


where ty is the frictional stress at the surface (if a sur 
face 1s present Since /i,, is constant, the energy equa 
tion is automatically satisfied under the assumption (ii 
above. 

The momentum equation |Eq. (2.2)| represents the 
integral of the exact equations of motion across the flow, 
except that the gradients of the viscous or Reynolds 
stresses in the flow direction are considered to be negli- 
gible in comparison with the pressure gradient. Up to 
this point the approach does not differ essentially from 
the von Karman momentum integral method of ordi 
nary boundary-layer theory. 

Let an average velocity of the internal flow be defined 


by the relation 


1.e., the ratio of the momentum flux to the mass flux. 
Let us also define, without attributing any particular 
physical meaning to the definition, a mean density py, 
or a mean temperature 7), through the mass flow rela 


tion 
mM = pid = (puj6)/RT7) 2.4 


Introducing the displacement and momentum thick 


nesses of ordinary boundary-laver theory + 
i pu ; 
6* ( ] ) dv = ¢0 
J peu 
° pl 4 ™ m 7 
dv = 6 2.9a 
Jo peu p,u 
i "6 pu u f : 
oO” / ( | dy = 6 6* 
J0 Pelt U, 


‘ pu . i t I rs 
/ : dy = 0 é* 2 5b 
J Pette’ pP.Ue 


From these equations, # and / can be derived and re- 
placed in Eqs. (2.3) and (2.4) to obtain the following 


t 6** has been used here to indicate the momentum thickness 


in order to avoid confusion with the angle of flow, @ 
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expressions for the mean quantities in terms of the 


thicknesses: 


- 6**) 


U/ue=« = (6 —6 (6 — 6*) (2.6a) 


pi/pe = T/T, = (6 — 6*)?/6(6 — 6* — 6**) (2.6b 


where x = //mu, is an important parameter always less 
than 1. The definition of the thickness 6 will be dis- 
cussed in Section (2.2). 

Nondimensional reduced velocities w = u/a, are now 
dividing through by the stagnation 


The auxiliary quantities 


introduced by 
speed of sound ay. 


qi = T; ¥ if l yw) | 
rT, | Ll — [(y — 1)/2)}e,? (2.6¢) 
ee ae = 
‘ Ty, yw, 


are defined, in terms of which the Bernoulli equation 


can be written as 


dw . 
dp/p = - F (2.6d) 
¢ 
Useful relations are, also, 
> : pw 
puAa, = : pu. = (2.6e) 
Ce Pe 
The reduced mass flow is then given by 
m = ma, = pi gq = 1/w (2.7) 
and Eqs. (2.1) and (2.2) become 
dm dé p {dé 
= p,U,as —O#)}= — 6 (2.8) 
dx dy ¢g, \dx 
d dm _dp 
(Mw) = We —_*_ — Ty = 
dx dx dx 
dm _dp pw, C 
W, — ig — (2.9 
dx dx ¢g 2 


where the wall friction ry = (1 2)p,u,*c, has been ex- 
pressed through the friction coefficient c,. 
now that c, can be related directly to the unknowns 
Suppose also that 7) 7’, can be ex- 


Suppose 


m, 6, @,, and @). 
pressed explicitly in terms of the unknowns through a 
more general formula than the ordinary one-dimensional 


energy relation 7,/7, = 1 — [(y — 1)/2]w,? (correct 
lor a uniform velocity distribution on/y). Then, since 
b= p(w,) and 6 = @(w,) are known by the Prandtl- 


Meyer relation, Eqs. (2.7), (2.8), and (2.9) furnish 
three relations between the four variables m, 6, w,, and 
w. An additional relation between these variables is 
introduced by means of an assumption regarding the 
mixing rate dm dx, which will be discussed in Section 
2.3). It is at this point, with the introduction of an 
approximate expression for the mean temperature and 
lor the mixing rate, that the present approach departs 
irom the usual boundary-layer theory. 
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2.2) Expression for the Mean Temperature: The F-k 

Curve 

In any given flow the mean temperature can be com 
puted from the defining relations, Eqs. (2.6) and (2.7). 
However, since the solution is not known, the exact 
evaluation cannot be obtained a priori. Therefore, 
if any progress is to be made toward a solution of the 
equations, an approximate relationship must be intro 
duced, which depends on some average properties of 
the flow and involves the minimum possible number of 
parameters. 

Laminar Flow.—In the laminar case, such a relation 
ship can be established from the following considera- 
tions: Stewartson® has shown that, for Prandtl num- 
ber unity, 4 ~ 7, and an insulated wall, any compres- 
sible boundary-layer flow with prescribed variation of 
external velocity can be reduced to an equivalent in 
compressible boundary-layer flow with a transformed 
external velocity distribution. This means that, from 
a given incompressible flow, a family of compressible 
flows can be generated whose solutions are related in a 
simple way to the incompressible solution. The Stew 
artson transformation is represented completely by the 


relations 
pity dn = pa, dy; pyay dt = pa, dx | 
(2.10) 
p p 
w= P= ves m= vy t= —Vel 
po Po 


where the subscript 0 denotes some fixed reference 
condition and the subscript 7 denotes the quantities of 
the equivalent incompressible flow with the constant 
density po, with coordinates £, 7, and with the same 
stream function y as in the compressible case. From 
the relations (2.10), the two following additional rela- 
tions can be written: 

Mi; u 


= -—; pot; dn = pu dy 
ay a, 


(2.10a) 


It follows that at corresponding stations the relative 
velocities in the compressible and incompressible flows 


coincide. 


U;/Uie = U/U, (2.11) 


By utilizing the first two relations in Eq. (2.10), the 


expressions corresponding to Eqs. (2.5a) and (2.5b) 
can be written as follows: 


6; 0 
4 r u; pA, ® pu 
6, — 6,* = dyn = [ dy = 
Jo uy polo Jo pelt, 


Pelle .. " 
(6 — 6*) (2.12 
poly 
ied pel ba pu- 
: 5 5? .. ere 
6; — 6,* — 6,** = J — dy = / -dy = 
0 Uj, polynJ 0 pu,” 
Palle _. . . . 
(6 — 6* — 6**) (2.13) 


Polo 
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py ul 


’ é : 0 nto a) 
0 = / ay = / dn =" 
J 0 FO pu pa, 0 é 
oo T, "6; T 
ee — dy (2.14 
pa, is 0 1 


However, from the constaney of the stagnation en 


dn 


mj) 


thalpy, 
r/T, = 1 — [y — (1/2) ew’ 


Therefore, utilizing Eqs. (2.11) and (2.13), one finds 


fi 
oe F, 


°6 ml _ | ? 
dn = | (1 w.2— ) dy = 
J0 2 i,” 


(6 6* 0 (2.19 


Finally, by employing Eqs. (2.12) and (2.13), Eq. 
(2.6a) can be written as follows: 
ity 6; — 6 © ae : 
kK = = : oe 2.16 
i, 0 <a. 


By means of Eqs. (2.14) and (2.15), Eq. (2.6b) is trans 


formed into 


T, 6 = 6 * 6 ¥* ¥y a l 
S S S S *K* 
_ = We." \0; —0O eae | 
i ie (6; — 6,*)? 2 
or, utilizing Eq. (2.16), into 
T y— 1 y-— 1 
» T " > / ) 
T = J = ; KW," = I ; Ww" iZ.13 
where 
f = (6; — 6,;* — 6,**)6,/(6; — 6,*)? (2.18 


This mean temperature-mean velocity relationship 
differs from that for a one-dimensional isoenergetic flow 
only in that the quantity f appears rather than unity. 
It is to be observed that, for every incompressible 
boundary-layer flow, f [Eq. (2.18)] and « [Eq. (2.16) ] 
can be related to each other, so that the given incom- 
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pressible flow and the corresponding complete family of 
compressible boundary layer flows obtained through th 
Stewartson transformation are characterized by a cer 
tain f(x) relation. All the effects of Mach Number for 
the compressible family are contained in the second 
term of the right-hand side of Eq. (2.17). Of course, 
in general, this f(x) relation is different for every flow 
However, at this point it should be observed that cer 
tain classes of low-speed laminar boundary-layer flows, 
such as those investigated by Falkner and Skan 
ui, = Cé), are described completely in terms of 
single parameter. Each flow corresponds to a_ point 
in the f-« plane, and the whole class of flows is repre 
sented by a single f-x curve. By means of the Stewart 
son transformation this one-parameter description is 
extended to the corresponding compressible flows as 
well. 


The f(x) relation can be evaluated for known incom 
pressible solutions once the thickness 6; of the dissipa 
tive region is suitably defined. Now the thickness of a 
laminar boundary layer is not a sharply defined quan 
tity, because of the asymptotic transition from dissipa 
tive to nondissipative conditions. Too large a value oj 
6, would result in the inclusion of too much of the nearly 
isentropic flow in the dissipative region, while for too 
small 6; appreciable viscous effects could be lost 
Clearly, the choice of 6; is influenced by the relative im 
portance of the transversal pressure variations in the 
flow {see Section (6)]. For comparison two different 
values of 6; are employed in what follows, correspond 
0.99 and 0.95. 


Once 6; 1s defined, « 


(2.16) and (2.18) by 


ing to u(6)/u, = 
and f can be obtained from Eqs. 
** computed 


utilizing the ordinary values of 6,* and 6 


for y—> which corresponds to concentrating all of 
the actual mass and momentum defect of the boundary 
layer in the conventionally defined dissipative region 
Observe that the definition of 6 follows immediately 


from that of 6; because of Eq. (2.11) 


In Fig. 2a the f-x curves are shown for two laminar 
incompressible boundary-layer solutions: the How- 
arth'* solution for a linearly decelerated velocity distri 
bution, and the Falkner-Skan'" |! family (a;, = C&) ol 
solutions. The first curve corresponds to the develop- 
ment of a specific boundary-laver flow, with « going 
final 


value to a minimum 


this 


from an initial maximum 


curve 


value where separation occurs; passes 
through the f-« point corresponding to the 
solution on a flat plate, which is characterized by con 


On the other hand, each pomt 


Blasius 


stant values of f and x. 
of the second curve corresponds to a different boundary 
layer flow, evolving in such a way that f and « remail 
constant and the “‘position”’ of the flow with respect to 
Blasius 


separation remains fixed. Here, again, the 


solution is included; all the points on the right of the 


Blasius point correspond to acceleration; those on the 
left, to deceleration. The minimum « corresponds to 4 


flow that is continuously on the verge of separation. 
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DISSIPATIVE FLOWS 


\ll of the curves show a common qualitative behav- 
ior: going toward separation « decreases and f increases 
steadily, the slope becoming steeper and steeper until, 
at separation, the slope is practically vertical—.e., x 
varies but little as f increases. This fact has a clear 
physical interpretation. The increase in thickness near 
separation is caused mainly by an increase in the extent 
of the low-velocity portion of the boundary layer, the 
high-velocity region being pushed away from the wall, 
while the velocity distribution in the high-velocity re- 
gion remains approximately similar to itself. In the 
extreme case, when the velocities near the wall are very 
low, the effect of this region on J and m is negligible, 
and x is nearly constant at the value computed from the 
high-velocity region alone. On the other hand, f in- 
creases steadily because of the increase in total thickness 
of the dissipative flow. This fact can also be expressed 
in the following way: With approaching separation, 
the effect of the wall on the high-velocity portion be- 
comes less and less marked, and the conditions in this 
portion, which determines the value of x, become more 
and more similar to those of a “‘free’’ half-jet with a 
constant « value. This physical picture is essential, 
since it gives a clue to the qualitative prediction of what 
happens after separation, where no theoretical solutions 
are known. It is to be expected that in the separated 
region the jetlike behavior is still more pronounced, 
and therefore x undergoes only small variations, while 
The 


f-« curve will have therefore a nearly vertical exten 


all the variations in the flow will be reflected in f. 


sion beyond the separation point. 

An interesting case is one of a flow that, after having 
undergone separation in the way just described, is 
brought to reattach itself to the wall and to build up a 
regular boundary layer again. A compressible flow of 
this type is observed, for example, in the laminar 
boundary-layer-shock-wave interaction. Again, de- 
spite the absolute lack of theoretical information on 
flows of this kind, the qualitative behavior of a ‘‘re- 
attaching’ flow in the f-« plane can be predicted to 
follow a path analogous to the one observed for sepa- 
rating and separated flows, but in the opposite direction. 
The f-« curve, after having described some loop with 
small x variation and large f-variation, reaches a re- 
attachment point, generally distinct from the sepa 
If, for 
instance, the conditions are uniform after reattachment, 


ration point, after which « starts increasing again. 


a Blasius flow will finally be reached, and the f-« curve 
igain passes through the Blasius point. 

Consider finally the case of a wake formed behind a 
body, where the velocity distribution across the dissi 
pative flow region passes from a typical boundary-layer 
profile just upstream of the base of the body to a neces- 
sarily unjform distribution in the wake far behind the 
body (e.g., Fig. 7). Far beyond the body, 6;* = 6,** = 
YVand« = f = 1. In the f-« plane the path of the flow 
will be represented by a smooth line connecting the 
initial f-« point (say the Blasius point) with the point 
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x = f= 1. Fora blunt-based body, as in the compres- 
sible base pressure problem, conditions in the wake just 
aft of the body are similar to those of a separated flow, 
the low-velocity region being in this case represented 
by the ‘‘dead-water”’ region. The path in the f-« plane 
will be developed at first as in the case of reattachment 
(where, however, the wall influence is negligible 
through a line of nearly constant x. 


increasing, until far beyond the body x and f—~1. The 


Then « will start 


f-« curve for this flow does not necessarily pass through 
the Blasius point, but it is not likely to pass too far from 
it. 

From this discussion it follows that every incompres 
sible dissipative flow involving separating, separated, 
reattached regions and wakes can be described as a par 
ticular curve in the f-« plane having qualitatively the 
shape just discussed. Supposing that the Stewartson 
transformation still holds in the separated regions, the 
family of corresponding compressible dissipative flows 
is described by the same curve. The exact determina 
tion of this curve requires the complete solution of the 
flow, which can be obtained only in the case of unsepa- 
rated flows, since no known method is available for the 
other cases. The next step is therefore to replace the 
exact curve with a suitable approximation. 

Based on the preceding qualitative considerations, it 
appears that the loop characterizing the separation 
reattachment case can be replaced with a single repre- 
sentative mean curve, extended vertically (x = con- 
stant) beyond the separation point (jetlike behavior in 
separated region). For instance, the single mean 
curve can be the theoretical f-« curve for the Falkner- 
Skan solutions. Moreover, this curve can be extended 
up to f = « = 1, and in this way it includes also the 
wakelike flows. However, the value of « at separation 
given by the Falkner-Skan solutions does not coincide 
with the x; value for the free jet given by Chapman's 
solution,!® as shown in Fig. 2a. One could neglect this 
small difference in x and extend the Falkner-Skan curve 
vertically or, as we have done in Fig. 2a for u(6)/u, = 
0.95, utilize an extrapolation between the Falkner- 
Skan f(x) curve and the line « = «; given by Chapman’s 
solution. [In order not to complicate the figure, we 
have not included the curve for u(6) u, = 0.99. ] 

With the assumption of this general relation between 
f and x, the character of the flow, even if slightly dis- 
torted, is certainly preserved. Thus, f(«) in Eq. (2.17) 
becomes a general function for all types of incompres- 
sible flows and therefore for all types of compressible 
flows covering the cases of separation, reattachment, 
and wake flow. In other words, the assumption is made 
that a single parameter, « (or, better, f), is sufficient 
to characterize the state of the flow through all its pos- 
sible development. Of course, all the assumptions and 


+ This procedure may recall to a certain extent the method of 
Thwaites,'* who used a single mean function derived by com- 
parison between known solutions in his successful treatment of 
incompressible laminar boundary layers with a pressure gradient 
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quantitative considerations regarding the f-« curve 
should be thought of as preliminary and may be sub- 
ject to some revisions after analysis of existing (and 
new) experimental data. 

It proves to be more convenient later to work with 
the quantity F, which is related to f by 


+ 6 
; (2.19) 
In Fig. 2b the corresponding F-« curves for laminar 
flows are shown. (Note that, in wakelike flows, when 
fandx«—~>1, F—~0. 

Turbulent Flows.- While the laminar case gives the 
possibility of determining the important relation (2.17 
and the approximate behavior of the quantities involved 
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from theoretical considerations, nothing of this kind 
can be done in the turbulent case because of the almost 
complete lack of theoretical knowledge on the effect 
of compressibility and because even the empirical in- 
formation in the incompressible case is meager. ‘The 
possibility of considering the velocity profiles for an 
incompressible boundary layer as a one-parameter 
family has been discussed many times, but it is not vet 
clear whether this assumption is even approximately 
correct in general. If the velocity profile is not chang 
ing too fast, then this assumption may be correct. In 
that case the quantity «x defined by Eq. (2.16) can again 
be chosen as the parameter representing the state of 
the flow, and the function f is computed from Eq. 


* 


(2.19) by using known incompressible values of 6,;* and 


6 baht 


the low-speed data of Schubauer and Klebanoff. It 


The curve of Fig. 3a has been computed from 


should be observed that in the turbulent case it is pos 
sible to take 6 such that u(6),u, = 1.0, because the 
thickness is a sharply defined quantity. | 
culated from the von Doenhoff-Tetervin semiempirical 
The behavior of this f-s 


Points cal 


relation’? are also shown. 
curve, which corresponds to boundary layers going to 
ward separation, is similar to that of the laminar case, 
and the same qualitative discussion can be applied 
The f-« path corresponding to any given low-speed cas¢ 
involving separation, reattachment, or wakes can _ be 
represented approximately by a single curve, consisting 
of the curve just described extended vertically beyond 
extrapolated toward f = 


separation (x = «,),! and 


x = | to cover the case of wakes. The curve passes 
through the point corresponding to zero pressure gradi 
ent (flat plate). 

Finally, the additional bold assumption is made that 
the effect of compressibility can be taken into account 
exactly as in the laminar case by utilizing Eq. (2.17 
As a preliminary test of this assumption, good qualita 
tive correlation is observed between the Schubauer low 
speed f-« curve and the points calculated from turbulent 
boundary-layer profiles measured through a shock-wave 
interaction region at ./ = 3.0 by Bogdonoff and Sol 
arski'® in this laboratory (Fig. 3a). Some of these 
points correspond to the region preceding the incident 
shock wave, where «x is decreasing, and some to the 
region downstream of the incident shock, where the 
velocity profile is gradually returning to its undisturbed 


form. 


2.3) Mixing Rate and Friction Coefficient 

Again, the laminar case can be treated by eliminating 
compressibility effects with the aid of the Stewartson 
transformation and examining known incompressible 


solutions. By integrating the second relation of Eq 


t For a semi-infinite ‘‘half-jet’’ at constant pressut 
0.720 by direct calculation from the theoretical (Tollmien 


So little is known about turbulent flows after 


veloc: 
ity distribution 
separation that the vertical extension of the f-« curve at this value 


of x should be regarded only as a first approximation 
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DISSIPATIVE FLOWS AN 


(2.10a) across the boundary layer, it appears that 
m;, = m. Therefore, using the second relation of Eq. 
(2.10), the compressible and the incompressible mixing 


rates are connected by the following equation : 


dm dx = (a,.p, dopo) (dm,/dé) (2.20 
Now, for an incompressible flow [Eq. (2.5a) ], 
mM; = pott;,(6; — 6,* (2.21) 


and on a flat plate 


IC, l /o 

: P 2 Co 

6, — 6," = ( ) (2.2la 
polt ;, 


where Cis a constant whose value depends on the defini 


Str 


tion of 6; The expression for dm, d& can be found 
from Eqs. (2.21) and (2.2la) and substituted into Eq 
(2.20) to obtain the value of dm dx. If a mixing coef 


ficient is introduced such that 
dm dx = kp,.u, (2.22) 


The 


analogous expression for the incompressible case is 


the result for the flat plate is k C(u, Wt). 
k C(uo m;), and the nondimensional expression 
m, 4 here replaces the square root of the Reynolds 
Number {see Eqs. (2.21) and (2.2la)]. The quantity 
k, can also be evaluated for various other solutions of 
the laminar boundary layer; the form of the relation 
giving k; remains the same, except that C is now a func- 
tion of the parameter determining the conditions of the 
boundary layer. 
C= 
for the corresponding compressible family, and there- 


If the parameter « is chosen, then 
C(x) is the same for the incompressible case and 


fore the mixing rate is given by Eq. (2.22) with 


m)C(K) (2.23) 


e = (ML, 


The values of C(x) depend, of course, on the defini 
tion of 6,;; for the theoretical solutions already dis 
cussed, C(x) is given in Fig. 4 for the two definitions 
already used in the f-« computations. 

A similar treatment can be given for the friction coef 
ficient. From Eq. (2.10a), the following relation is 
derived between the viscous stresses: 
pe 1 pb. T, 

r. poli, = bo To T 


T = wu, = 
where the assumption «» ~ 7 has been used, On the 
other hand, 


Pelt,” = (Pe/ Po) poli? 
therefore, 


; 2r ; le 2 Me 
peu,” To potie? Lo 


Now, for the flat plate, remembering the relations 
between #7; and Reynolds Number, one finds that 
(i; = D(uo/m;,), with constant D. For the other in- 


compressible solutions, the same relation can be writ- 
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ten, with D = D(x), the function being generally differ 


ent for different solutions. Finally, the relation 


cr = (u,/m)D(x) (2.24) 
gives the friction coefficient for the corresponding com- 
pressible families. D(x) is plotted in Fig. 4 for the 
Howarth and Falkner-Skan solutions. 

The next step now is to replace the particular func- 
tions C(x), D(x) representing particular solutions with 
two general functions, which have to be extended be- 
yond separation and into the wake region and must also 
Provisionally, C(x) and 
Falkner-Skan 


cover the reattachment case. 


D(x) obtained from the solutions are 
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chosen as “universal” functions. However, further 
study is required to determine the mixing coefficient 
for the jetlike flow beyond separation and prior to re- 
attachment and also for wake flows with a pressure 
gradient. Of course, D(x) = 0 beyond separation, up- 
stream of reattachment, and in wake flows. 

Again the turbulent case can be treated only by 
analogy. For the flat plate in the incompressible case, 
the nondimensional expression 777;/uo is again a function 
only of the Reynolds Number. If it is assumed that 
the '/;th power law holds for the velocity distribution, 


then two relations can be written: 


Mo ; 1 Mo , t 
k; = c( ) : Ch; = D( 2 ) 
My m 


with constant values for C and D. Equations analo 
gous to Eqs. (2.23) and (2.24) could be written, where 
C and D could be considered as functions of x and yu, m 
would appear to the ' yth power. However, it is 
doubtful whether these relations could give a correct 
answer, including the effects of approaching separation 
and of Mach Number. For the present, considering 
that the mass flow variation seems to have but little 
effect on the turbulent values and that insufficient ex- 
perimental data are available on the spreading of tur 
bulence with varying pressure, it seems justified to in 
troduce the approximation that & is constant through 
out the turbulent region. The skin-friction coefficient 
c, is neglected in separated flows or is taken as a simple 
function of « in separating or reattached flows. Of 
course, the same rough assumption could be made in 
the laminar case, especially when only narrow regions 
of flow are considered. In physical terms the approxi 
mation that k is constant amounts to the assumption 
that the rate of mass transport from isentropic to dissi 
pative flows is proportional to the mass flux density of 
the isentropic stream [see Eq. (2.22) ]. 

So far as the order of magnitude of k is concerned, for 
the boundary layer on a flat plate with zero pressure 
gradient, Eq. (2.9) yields the relation 


Ro Ce 2(1 Ko 


where the subscript zero here denotes conditions at 
zero pressure gradient. At a local Reynolds Number 
of 10° and a Mach Number of 2.0, one finds ky = 1.6 X 
10~* for the laminar boundary layer with u(6)/u, = 
0.99, while k) = 1.5 K 107° forthe turbulent layer. In 
turbulent separated flows and jets somewhat larger 
values of k are expected. From the Schubauer and 
Klebanoff data* one finds that k — 0.03 near turbulent 


First, with the aid of Eqs. (2.6c), (2.17), and (2.19), the mass flow equation [Eq. (2.7) 


mk) F+ 1 — 


By differentiating the value of 6 obtained from this equation with respect to x, dividing 
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boundary-layer separation, and this value is utilized 
for turbulent flows in the present paper as a first ap- 
proximation. The fact that the turbulent mixing coef 
ficient k is of the order of ten times the laminar value is 
largely responsible for the marked difference between 
turbulent and laminar interaction phenomena. 

The constant mixing coefficient introduced in this 
way should be regarded only as an average value, the 
introduction of which makes it possible to decide 
whether mixing furnishes the main mechanism for the 
pressure rise in interactions between dissipative and 
isentropic gas streams. Refinements can be introduced 
later when this mechanism is firmly established and 
when more experimental data are obtained on the mix 


ing process with varying pressure. 


2.4) Reduction of Equations of Motion to a Single Equa- 
tion for Constant Mixing Coefficient 


By introducing the relation (2.17) between mean 
temperature and mean velocity and the additional Eq 
(2.22) for the mixing rate, which can be also written as 


(2.2: 


= k(p/¢,) ye 5" 


dm/dx 


the equations describing the flow are completely defined, 
provided the coefficients k and c,; are known. In what 
follows, the discussion is restricted to the simplest case 
of constant k mentioned in the preceding section, which 
permits a great simplification in the discussion of the 
important properties of the equations and, moreover, 
at the present time seems to be the most reasonable as 
sumption in the case of turbulent flow. On the other 


hand, c, can be replaced by the quantity 


¢ = ¢; PRI kK 

with the advantage that, even when & and C, are as 
sumed to depend on m (as, for instance, in the laminar 
case [Eqs. (2.23) and (2.24)]} the effect of m is elimi- 
nated in o, which is therefore a function of « alone. For 
separated flows or wake flows, c, is negligible or zero 
For the flat plate, o = 1. 

is reduced to the 


and therefore o = 0. 
By utilizing Eq. (2.25), Eq. (2.1 


suggestive form 


d6 dx =O0+k 2.26 
The distance x along the surface is now eliminated from 
Eqs. (2.7), (2.8), and (2.9) with the help of Eq. (2.25), 


and the momentum and continuity relations are re- 
duced to two simultaneous equations for dx ‘d In m and 


dw,/d\1nm. These equations furnish a single equation 


for dx dw, or for dF dw,. 


is rewritten in the form 


we | = ypwd 


] } 
through by dm/dx, and 
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using Eq. (2.26), the following equation is obtained : 


i 4 : y-1 ,N - y— d (? ) 
<x] F+ 1 WW,” — F+i1- W," (In pw.) = +1} 
1 In 44 2 i} «| 2? dinm me k 
¥~1 — 
(1 - —— w,?) lr bi~ = we | (2.28) 
& 


: _ ieoall 
x) F+1—- Ww, 
d 2 din w, | (2.29) 


(KW, } — ae 


dinm Y¢- din m 


- 


to 


where the Bernoulli equation [Eq. (2.6d)] has been employed to eliminate the pressure. 
So long as the determinant of the simultaneous equations {Eqs. (2.28) and (2.29)] for dx/din mand dw,/d In m 
does not vanish, these equations yield a single nonlinear first-order equation for 


(dx dw.) = (dk din m) (dw, dinm) 
or for dF dw, if the relation 
dk /dw, dx dF)(dF dw, 


is employed, 


| im (2.350) 
u 7 : 
} 1 w, ; 
+ 2) a) I dk 
W, l wt’, } — 
2 «(il —« k K( I h dk 
em g 
i 
F+1 i? 
2 dk 
l+oa 7 l 
iN dk 
Since k, o(x), and F = F (x) are presumed known and @(w,) is given by the Prandtl-Meyer relation (or by any 
other procedure for nonisentropic flows), this equation is integrable by numerical methods. 
When « = constant and o = 0, as in separated flows or wakes, Eq. (2.30) reduces to 
r y¥—! 
I Ww, 
di y+ 1 I k 2 6 2.31 
w = -" Ww, — |] os yw, I T 
dw, 2 1 | l k k k 
l Ww," 
» 
which is a Riccati equation. A ‘simplified’ theory has been developed in Section (4) for this important special 
case. 
When the determinant of Eqs. (2.28) and (2.29) vanishes, the numerator and denominator of the right-hand side 


ff Eq. (2.30) also vanish, and the equation for dF dw, hasa “‘critical point’’ that turns out to be significant for the 


determination of the proper solution. 


25) Critical Point of the Equation for dF dw, and Its jnator of Eq. (2.30) equal to zero and solving the re- 
Significance . ' ‘ ‘ : 
sultant simultaneous equations for F,, and w,,, with 
Suppose, for simplicity, that only wake flows or sepa- P . , ; . 
tapi nes pl . : : the aid of the Prandtl-Meyer relation, the location of 
tated flows for which o = O are considered here. a ae: z ; ‘ a = 
. eee ok the critical point is determined as a function of w, or .V/, 
Flows in which skin friction is not negligible can be in- oe ; ; 
duded in the analysis and will form the subject of a where .\/ is the Mach Number of the free stream far 
later paper.) Suppose also that the external flow is from the body where 6 = 0. When the denominator 


isentropic. By setting both the numerator and denom- vanishes, 
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eg Oe meen ae Ke(l — rod(<) + F., = (1 oe w2) 
om ak J e 2 k 
9 25) 
95, 24 (38 
2.0 provided dk/dF # 0. If dk dF = O or «x = constant, 
9-5 f the denominator never vanishes and no critical point 
‘ - exists. When the numerator also vanishes, ] 
ie! ae dF 
8 14 ' 
+ —_ C = in? — i)(. ee twee?) (2.33) ; 
~ (eI : 2 2 
-| P 
10 where Eq. (2.52) has already been used to simplify , 
7LO. (te wee 3 i a the result. Eqs. (2.32) and (2.33) define 6,,/2 and w,.,, | 
lo ia 14] 16 18 [20 22 24 W : ; - t 
10 15 20 30 40 ° M as functions of x,,.. Since @,,,k = O(1), when R is very 
small, as in the laminar case, 6., ~ O and w,, ~ W,. In 
FIGURE 5. VALUES OF 2€, 6, We AT CRITICAL the turbulent case, however, k is not small. For ex- it 
POINT FOR TURBULENT FLOW. k=0.03 ample, choosing k = 0.03 on the basis of the Schu- b 
bauer-Klebanoff data near turbulent separation’ and ‘i 
utilizing the F-« curve of Fig. 3a and the dx/dF curve P! 
3.0 of Fig. 3b, the values of 6,.,, x,-, and w, shown in Fig. 5 v 
are obtained as functions of @,. It appears that 6,, is tl 
“a always negative. th 
Now, what is the significance of the critical point? of 
20} Since 6,, < 0, it is necessary to consider only ‘“‘reattach- la 
| @.< OE ing flows,’ such as the flow aft of a supersonic airfoil 
‘ al » r , aiid with a blunt trailing edge or the reattachment of the 2 
Gi. = O separated flow downstream of the reflected expansion 
‘ol i. ane “fan” in the laminar boundary-layer-shock-wave inter- 
action. 
| / - In the first case, for example, x must increase, or F gel 
5] must decrease and w, must decrease steadily during the na 
1 a] — recompression that occurs in the wake after the expan- ble 
cr, me 190 os sion around the sharp ‘corner’ at the trailing edge Co 
We In other words, dF dw, is always positive (see Fig. 6). pal 
FIGURE 6. SIGNIFICANCE OF CRITICAL POINT, Now, in the separated or ‘‘free-jet’’ region just aft of the 
SHOWING INTEGRAL CURVE FOR We50 the base, x is constant and dx/dF = 0, so that the “et 
a ee Tee denominator in the expression for dF dw, in Eq. (2.30) Ma 
is initially negative. Therefore, the numerator must reg 
DETAIL: FLOW also be initially negative. If 6;, or the flow inclination stre 
AROUND "CORNER" angle after the expansion around the corner, is /arger K al 
. - than the proper value @,, then the denominator will mir 
ly es / vanish before the numerator during the recompression Int 
eH & tf / (since dk/dF < 0), and the slope dF dw, will become por 
mR | infinite at some point and negative thereafter. In that ups 
/ / case, the F(w,) integral curve can never represent the nun 
J actual flow (Fig. 6). On the other hand, if 6; is smaller criti 
_ BOUNDARY LAYER ~ f Sen than the proper value, the numerator vanishes before 
kK 6 yi —<- = — the denominator, and dF dw, passes through zero and Fr 
KL hed ee ree becomes negative thereafter. Again, the F(w,) curve 
. vine cannot be the correct representation of the flow. Only 
\ POINT one value of 6;, say ,-, permits the integral curve to pass 
\ \ through the critical point where numerator and de- +] 
\ nominator vanish simultaneously. Thus, the critical catio’ 
\ point acts somewhat like the throat of a nozzle in de- revol 
Fic. 7. Schematic representation of flow over supersonic airfoil! termining 9, and the base pressure for the given bound- sd 
with blunt trailing-edge. arv conditions at the base [Section (3) ]. a 
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It is an interesting fact that the location of the critical 
point does not, in general, coincide with the sinallest 
cross section of the mixing region—that is, with the 
point at which di dx = 6+ k =0Oor@k=— 1. By 
referring to Fig. 5, one sees that in the turbulent case 
§.. k > — | for 1 <@®, < 1.61, so that the critical point 
lies downstream of the narrowest cross section, while 
§.. Rk < — 1 for ®, > 1.61, so that the critical point in 
this range lies upstream of the narrowest cross section. 
In any case, most of the recompression occurs before 
the critical point is reached. 

Evidently, the mixing plays an essential role here be- 
cause, fork = 0, w,(dF/dw,) is negative and no critical 
point exists. An analysis of compression flows of this 
type which ignores the mixing process cannot possibly 
represent the actual conditions. 

The discussion of the role of the critical point (when 
it exists) applies also to the reattaching flow in the 
boundary-layer-shock-wave interaction phenomenon, 
even when o ¥ O, although the discussion is more com- 
plicated and will be given in a later paper. In the case 
of a turbulent boundary-layer-shock-wave interaction, 
the flow problem is undoubtedly modified not only by 
the action of skin friction but also by the penetration 
of the incident and reflected shock waves through the 


layer [see Section (6) ]. 


(3) APPLICATION TO THE BASE PRESSURE PROBLEM IN 
PLANE SUPERSONIC FLow 


In order to illustrate the concepts and methods, the 
general theory is applied to the problem of the determi- 
nation of the base pressure for a supersonic airfoil with a 
blunt trailing edge and a turbulent wake (Fig. 7).T 
Consider the simplest case in which the airfoil surface is 
parallel to the flight direction at the trailing edge, and 
the Mach Number just before the expansion around the 
“corner” is approximately equal to the free-stream 
Mach Number. 
regarded as isentropic, the Mach Number far down- 
The values of 


If the recompression in the wake is 


stream is again the free-stream value. 
«and F at the critical point are then completely deter- 
mined once the nondimensional mixing rate, k, is known. 
Integration of Eq. (2.30) may be started at the critical 
point and proceeds both downstream in the wake and 
upstream toward the airfoil trailing edge. Since both 
numerator and denominator in Eq. (2.30) vanish at the 
critical point, a parabola, defined by the equation 


dF 1 ({d°F 
F- F = ( ) wv - Ww.) + ( ) 
dw, , Bad 2! \dw,?/ ., x 


cr 


t Extension of the theory to axially symmetric flows and appli 
cation to the base pressure problem for blunt-based bodies of 
revolution is now in progress at this Laboratory 

t If the oblique shocks in the wake are sufficiently strong, 
shock-wave losses can also be taken into account in a more exact 


calculation. 
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is utilized in the immediate vicinity of this point. [The 


values of (dF dw,),, and (d*F dw,”),, are calculated by 
substituting the parabolic representation of F(w,) into 
Eq. (2.30).] 


curve of the form 


An analytical representation of the F - x 


kK = 0.92954 — 0.103S2F — 0.36842F? + 0.44575F* — 
0.09582 F* 0.09948 F° + 0.04505 F* 


is employed near the critical point, and appropriate 
polynomial approximations involving fewer terms are 
utilized away from the critical point, where the integral 
curve is less sensitive to the function x = «(/). Numer 
ical integration is begun as soon as the slope dF/dw, 
can be calculated numerically with sufficient accuracy 
from Eq. (2.30).fT 

Far downstream of the airfoil, as already discussed, 
« — | and F + 0. The proper formulation of the 
boundary condition at the trailing edge, on the other 
hand, requires a knowledge of the process by which the 
boundary-layer flow at the body surface just upstream 
of the trailing edge is transformed to the wake flow 
downstream of the trailingedge. No similarity param 
eter exists for this viscous flow, and the problem is 
further complicated by the interaction between the vis- 
cous flow and the “external” isentropic supersonic ex- 
pansion around the “corner.” In the F - « diagram, 
the actual transition 1-2-3, indicated by a dotted 
curve in Sketch A, consists first of the expansion 1-2, in 
which « increases because of the uniformizing effect of 
the expansion and F increases because the dead-water 
region is now included in the dissipative flow. A sub 
sequent smooth adjustment 2-3 occurs to a condi- 
tion close to that of a free jet. To be consistent with 
previous approximations [Section (2.2)], it is assumed 
that point 3 lies on the vertical extension of the mean 


F-« curve (Fig. 3a). 


tt It would also be possible to begin the integration just after 
the expansion at the trailing edge, with various trial values of 6, 
(or we;), and thus determine the value of 6; for which the integral 
curve passes through the critical point (see Fig. 6 In fact, this 
method is the only feasible one for the general case when k is a 
function of m and x, or for axially symmetric supersonic flows 
where the relation between w, and @ is not known a priori. 
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A rough calculation indicates that the ‘‘free-jet’’ flow 
is established in a distance less than one boundary 
layer thickness, 6,. Provided that the ratio of bound 
ary-layer thickness to airfoil thickness at the trailing 
edge is small, it seems justified to treat the problem ap 
proximately by assuming that the transition to the 
“free-jet’’ flow occurs over a negligible portion of the 
wake just behind the body and that negligible mixing 
occurs between stations 2 and 3—1.€., My) ~ my, The 
actual transition 1—2—3 is replaced by a “‘jump” 1-3 to 
the free-jet condition. In integrating Eq. (2.30) from 
the critical point toward the airfoil, the F 
followed until « equals «;, or the free-jet value, and the 
value of « is taken to be constant thereafter and equal 


to x; Of course, Eq. (2.30) is replaced by Eq. (2.51 


kK curve 1S 


in this region. 
To this approximation the boundary condition at the 


trailing edge is 6 = (h 2) + d,or 
6/6, = (h/26,) + (d/6p) (3.1 


where / is the airfoil thickness at the trailing edge and 
the value of d lies somewhere between 6,, the boundary 
layer thickness at the trailing edge, and 6,,, the thick 
ness of the ‘‘free-jet’’ mixing region (point 3) estab 
lished just aft of the isentropic expansion around the 
trailing edge (Fig. 7). Since d/6, is generally small 
compared with /1/26,, it does not need to be calculated 
with the same degree of accuracy as 6 6,. Therefore, 
the approximation d = 6,-, which seemed to be a reason- 
able one, was actually employed in the present calcula 
is obtained from the mass flow 


tions. The ratio 6, 6,- 


relation as follows: 


Myr = (Pelle)n Oy [1 — (6*/6) o> = My = 
(pelte)nOn{1 — (6*/6)], (3.2a 
or 
& [1 — (6*/8)- GM, 
= : ; (3.2b 
5» [1 (6* 6)], GCM, 
where 
1/2){ 
. Y = 5 ™ 
GCM) = (i + uv") 
and Jf, M, the Mach Number at large distances 


from the body. For a turbulent boundary layer with 


a '/;th power velocity distribution,* zero heat transfer 
at the surface, and constant stagnation enthalpy, the 


values of [1 — (6* 6)], are tabulated below 
V/ l (6*/6 
1.5 0.806 
2.0 0.768 
3.0 0.669 
@ 0 583 
5.0 0 506 
The quantity [1 (6* 6)|, is assumed to have the 


. 
value corresponding to a turbulent boundary laver 


¢ The meager data available indicate that the '/;th power veloc 


itv distribution is not too far from the truth 
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right up to the point where transition reaches the 


trailing edge. 

The value of [1 — (6* 6) |,’, calculated for a low-speed 
isothermal constant-pressure turbulent half-jet accord 
ing to Tollmien’s theory, is 0.45, and this value is 
utilized in Eq. (3.2b). 

With the aid of the assumption that m, = my, the 
ratio 6/6, appearing in the boundary condition [Eq. 


(3.1) | is obtained as a function of w.,, from Eq. (2.27 


as follows: 


Here, x = 0.720. To be consistent with the previous 
assumption regarding the effect of compressibility in 
turbulent flow [Eq. (2.17)], the values of x and F, 
at the trailing edge must be taken as identical to the 
corresponding values for zero pressure gradient at low 
speeds. Actually, it is believed that the variation of the 
zero pressure gradient values of F and x with Mach 
Number is quite small. 

By utilizing the results of the numerical integration 
of Eq. (2.30) for F(w,), together with Eqs. (3.2b) and 
(3.3), both 6/6, and 6,- 6, are determined as functions 
of Wey. Therefore, to each value of we, or 6,/, there 
corresponds a value of 6, 4 calculated from Eq. (3.1 


so that the base pressure ratio 


Po p f(0,+, We 

is determined as a function of 6,// and the mixing coef- 
ficient, R. Chapman’s correlation” of base pressure 
data in terms of 6,/h for an airfoil with a blunt trailing 
edge and a turbulent wake is now completely under- 
standable. From the nature of the integral curve 
F(w,), the relation between 6/6, and F [Eq. (3.3) ], and 
the boundary condition at the airfoil trailing edge |Eq. 
3.1)], one sees also that, for a given mixing rate in the 
wake, w+, is larger and the base pressure ratio p,’ p 1s 
lower the smaller the value of 6, /. 


For any given body shape 


6 6», ¢ ( - Ty 
- = =e" ¢( Re, Mt, — 
h Cn H r, 


where the functional relationship for 6, ¢ depends on 
the Reynolds Number range for laminar-turbulent 
transition on the body. ‘The base pressure is therefore 
determined as a function of Reynolds Number, Mach 
Number, Prandtl Number, and thermal conditions at 
the airfoil surface once the Reynolds Number range for 
boundary-layer transition on the body is known and 
once k in the turbulent wake is known. A detailed dis 
cussion of the various flow régimes in the base pressure 
problem will be given in connection with the applica 
tion of the ‘“‘simplified’’ mixing theory in Section (5). 
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lhe dependence of boundary-layer thickness 6, ¢ on 
Reynolds Number and Mach Number can be pre 
scribed as follows: 

\t sufficiently high Reynolds Numbers, boundary- 
laver transition occurs well forward of the trailing edge, 
and the turbulent boundary-layer thickness is approxi- 
mated by the following formula for zero heat transfer 


to the airfoil surface: 
0.037 Chy l 
( (Re. 


The ratio 6**/6 is evaluated by assuming a '/;th power 
velocity distribution and isoenergetic flow,+ and the 
“compressibility correction” (Cy,/Cy, is taken from 
(The factor 0.037 
Re is, Of course, the low-speed value of 6**/c. 
Values of Cyy/Cfy<, and 6” 

In this so-called “fully turbulent”’ régime it is clear 


Fig. 19 of Van Priest’s paper.'S 


* 6 are shown in Table | 
that the base pressure ratio depends only upon the 
parameter Re,’*(h/c) at a given Mach Number, as 
found by Chapman.* 7 

At sufficiently low Reynolds Numbers, on the other 
hand, the airfoil boundary layer is completely laminar 
and, for zero heat transfer at the airfoil surface, Pr = | 
andu~ 7,{ 


6» . 


= [5.60 + 1.20(y — 1) 7]/V Re, (+ 


( 


Two different types of transition between laminar and 
turbulent flow in the airfoil boundary layer are as 
sumed at intermediate Reynolds Numbers (Fig. 8). 

(1) Transition for a smooth body, with the transition 
point located at the trailing edge when Re, = 4 X 10°, 
and a gradual approach to fully turbulent flow, which 
is assuined to be attained at the trailing edge when 
Re. = 12 X 10°. Such gradual transitions have been 
observed experimentally by Wilson,'’ for example. 

(Il) Transition for a body with some sinall rough- 
ness or external disturbances, with the transition point 
located at the airfoil trailing edge at Re. = 3 X 105, 
and fully turbulent conditions approached much more 
rapidly than in (I), at Re, = 6 X 10°. 

The transition Reynolds Numbers are assumed to be 
independent of Mach Number. Values of F, and x, 
for the boundary layer at the airfoil trailing edge are 
taken to be those corresponding to a fully turbulent 

t Of course, it is possible to compute 6**/é also from Van 
Driest’s formulas,'* but, as no evidence exists that permits us to 
decide between the Van Driest velocity distribution and the 
‘5th power distribution, we have used the latter until reliable 


data are obtained 

t The value 5.60 appearing in Eq. (3.5) corresponds to u(6)/u, 
= 0.9975. To be consistent this value would have to be ad- 
justed to conform to the value of u(6)/u, chosen in calculating the 
lantinar f-x curve. For example, for u(6)/u, = 0.99, the constant 
is5.0. Calculations based on 5.60 had already been made before 
this point was noticed; however, the resulting differences are 
small, particularly in view of the other approximations introduced 


in the computations 
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TABLE 1 
WV 6** 6 Cru / Chi 
1.5 0.0842 0.870 
2 0 0.0773 0.810 
3.0 0.0635 0. 680 
$0 0.0495 0.600 
5.0 0.0895 0.510 


boundary layer until the transition point moves aft to 
the trailing edge. 

Since the wake is turbulent, it is justified to regard k 
as constant in first approximation [Section (2.3) ]. 
From the Schubauer-Klebanoff data,* k — 0.03 near 
turbulent separation, and this value is utilized in the 
present computations. With this approximation, the 
integral curve F(w,) obtained by numerical integration 
of Eq. (2.30) is independent of Reynolds Number at a 
given Mach Number, and this curve is utilized for the 
calculation of base pressure ratio in the entire Reynolds 
Nuinber range Re, 2 4+ X 10° [curve (I)], or Re. = 
3 X 10° [curve (II) ]. 
curve F(w,) for 17 = 3.0 is given in Fig. 6, and the vari- 


An example of such an integral 


ation of base pressure ratio with Reynolds Number cal- 
culated by the methods described above is plotted in 
Fig. 9 for transition curve (I) and c/h = 3.0. {The 
pvs. Re, obtained by the*simplified”’ mixing 
Section (4) ] is also shown. } 


curve of p,, 





theory to be described later 

A comparison between the theoretically predicted 
variation of base pressure with Reynolds Number in 
Fig. 9 and the recent experimental data of Bogdonoff”’ 
on a body of revolution shows excellent qualitative 
agreement, and it is now clear that the assumption that 
mixing is the dominant factor is correct. Evidently, 
the actual values of base pressure depend critically on 
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the mixing rate, and experimental determinations of the 
mixing coefficient are now urgently required. 

The general theory and methods outlined so far are 
applicable to any flow régime occurring in the base pres 
sure problem. However, when transition occurs in 
the wake [Re, < 4 X 10° for curve (I), Fig. S] a numeri 
cal integration of Eq. (2.30) is required for each value 
of the Reynolds Number, because the laminar mixing 
rate and location of transition in the wake both depend 
on the airfoil Reynolds Number. Calculations for this 
flow régime, and for all flow régimes for Mach Numbers 
other than 3.0, were therefore carried out by means of 
a ‘simplified’ form of the mixing theory developed for 


separated and reattaching flows in general. 


(4) “SIMPLIFIED MIxING THEORY FOR SEPARATED AND 
REATTACHING FLOWS AND WAKES 


4.1) Basic Assumptions and Equations 


Certain outstanding characteristics of separated and 
reattaching flows and wakes behind blunt-based bodies 
suggest that the treatment of these flows can be some 
what simplified with the aid of a few physical approxi 
mations. The aim of sucha “simplified”? mixing theorv 
is to furnish an approximate description of the main 
characteristics of separated and reattaching flows, in 
such a form that some general conclusions can be drawn 
and a clearer insight gained into the role of the mixing 
process in these flows. Of course, the results of the 
siinplified approach must be compared with those ob- 
tained by the more exact analysis in a few critical cases 
to test the validity of the approximations employed. 

Consider first the flow in the wake behind a super 
sonic airfoil with a blunt trailing edge (Fig. 7). Two 
flow regions may be distinguished: (1) a region in 
which the ‘‘external’’ gas streams from the upper and 
lower airfoil surfaces are approaching one another, and 
the ‘“‘dead-water”’ area just aft of the body is being ob- 
literated by mixing (most of the flow deflection and 
recompression occur in this region), (2) a region down- 
stream. of (1) in which the ‘“‘dead-water’’ region has 
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vanished and the external velocity u, varies only 
slightly. Because of the ‘‘mixing,’’ the velocity profile 
across the internal dissipative flow approaches a uni 
form distribution far behind the body -i.e., u; > u,, or 
«— 1, while u, > a, and 6— 0. 

In Region (1), « is nearly constant while f is varying; 
in Region (2), both f and « vary in general. Suppose 
now that the actual f-« curve in Fig. 2a or 3a is replaced 
by a broken curve, consisting of the extension of the 
vertical line x = constant down to f = | and the hori 
zontal segment f = | up tox = 1. The corresponding 
curve in the F-« plane is composed of a vertical line 
and an are of the curve F = (1x?) — las shown, for ex- 
ample, in Fig. 3a. Especially in the turbulent case, 
this approximation is not too far from the actual F-x 
curve. 

The two flow regions will now be distinguished more 
sharply by the conditions « = constant = x, in Region 
(1), and f = 1,«; S « S 1 in Region (2). The differ 
ential equation corresponding to the first region is the 
Riccati equation [Eq. (2.31)], the integration of which 
gives F(w,) in Region (J). By integrating Eq. (2.29 
with «x = «x, and o = O, the mass flux m is obtained asa 
function of w,. Finally, Eq. (2.27) gives 6 as a fune 
tion of w,, and by utilizing Eq. (2.26) both 6 and w, are 
obtained as functions of x in Region (1) 

In Region (2), the governing equation (to be given 
later) is that corresponding to Eq. (2.30) with f = 1.0 
and ¢ = 0. By properly connecting the solutions for 
the two regions at the junction point f = 1, x = «,, the 
complete (approximate) solution is obtained. 

Actually, a different and somewhat more direct pro 
cedure was employed, as follows: 

Region (1). 
proximately by 


In this region, the mass flux is given ap 


m = pi ¢ = (p ¢,)(6; — 6;* 1 
where 6; is the thickness of the “‘jet’’ or mixing region 
and 6,* is the displacement thickness, corresponding to 
the velocity and temperature distribution across the 
jet. In this region «x = constant = «;, and by utilizing 
Eq. (4.1) the momentum equation [Eq. (2.9) ] becomes 


1 dw, Ww, . 


gy, dm p 


where the Bernoulli equation [Eq. (2.6d)] is also em- 
ployed. Since di/dx = 06+ kanddm dx = k(p/¢,), tt 


follows that 
di/dm = |(¢/k) + 1]/(p/¢) $5 
When Eq. (4.3) is divided by Eq. (4.2) and the Prandtl- 
Meyer relation dw,,w, = —d@ V \J* — 1 for the ex- 
ternal flow is introduced, one obtains the relation 
dé l (0/k) +1 


= 5 — (5, — §.*)] $4) 
dé i—« warns” ia ~s 
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TABLE 2 
lhe Function /(.\/) in Radian Measure 
VW J(M V J(M 
1.0 0 2 70 0. 16989 
1. 40 0. O1L587 2.80 0.18118 
1.50 0.02461 2 90 0. 19209 
1 60 0. 03486 3.00 0.20262 
170 0.04604 3.10 0.21277 
1.80 0.05795 3.20 0. 22253 
1.90 0.070385 3.00 0.23191 
2 OO 0.083805 3.40 0.24091 
2 10 0.09588 3.30 0.24955 
2 2% 0. 10872 3.60 0. 25783 
2 30 0.12144 3.70 0.26577 
2? 40) 0.13397 3.80 0.27338 
2.0 0. 14626 3.90 0. 28067 
» 60 0. 15824 $ O00 0. 28766 


As a first approximation, it seems justified to regard 
the variation in the quantity «,(6; — 6,*) appearing in 
Eq. (4.4) as of secondary importance compared to the 
variation in 6. Actually, x; ~ 0.720 and [1 — (6* 6)]| 

0.45 for a low-speed turbulent jet, so that the quan 
— 6;*) is never greater than about (1 3)6,, 
* 


tity «,(6 
while 6 2 6). 
a constant, which corresponds to taking a constant 


Suppose that «;(6; — 6,;*) is regarded as 


mean value of 6;—for instance, the value 6g at the junc 
tion point with Region (2), where the effect of variations 


in 6, is most important. Then we may define 


> = (6 — «(d; — 6;*)] br 


So 


Eq. (4.4) is rewritten in the form 


1 dé (@k) + 1 | 


= (4.5 
6 dé V ji? -—11 K 
The solution of Eq. (4.5) is simply 
b/ér = @ MI — «dN, Ox (4.6a 
where 
” (@k) + ] 
1(0, Or) = dé (4.6b 
Jor VY Mr -— 1 
and 
be = 1 — «;{l — (6* 56)] 
By utilizing the characteristic relation 6 = vp — », 
where 


| 7 
\ MP? - | —_ 


tan-' yY V?- 1 


v(M) = \. tan ik 
q T 


and > = v(17/), the function /(@, 0g) is evaluated as 

follows: ° 

‘ l V W, = 

I(0, Oz) = [J(.M) — J(Mr)] - + 1) 1n (4.7) 
k k We p 

where 


- vdM 
J(iM) = 
Ji Mii +d(y -— 1)/2)47} 
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The function /(.\/) is tabulated in Table 2 and plotted 
against Mach Number in Fig. 10. 

By utilizing Eqs. (4.5), (4.6a), (4.6b), and (2.24), 
the variation of 6, or p, p with axial distance x, is read- 


ily calculated. One finds that 


d6 dt = (1 — « IRN Mr - 1] 6 (4.8 
or 
#— 2p | my * dé 
= = - (4.9 
OR 1 — x,Jor k \Or V Me 1 
where ¢ = x/dp. 


When transition occurs in the wake, the mixing rate 
parameter k is approximated by a step function, as 
shown in Sketch B (6,, denotes the value of 6 at which 
abrupt laminar-turbulent transition is assumed to oc- 


cur). For the range |@r;' S 6 SS |6,,|,the function 
1(0, Oe) in Eqs. (4.6b) and (4.7) is evaluated exactly as be- 
fore, with k = ki. In the laminar region, |6,,, S 


@\< 6,-\, the integral in Eq. (4.6b) is split into two 


parts, and 
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1(0, Or) = [(J(M) — J(M,,)] + 
‘lam 
l = l We), 
[J( Mi) — J(Mr)] — yp In + 
-turb Reurt Wer 
l W, W, 
- In - — In— (4.10 
Ria m Wetr Wep 


The integral appearing in the relation between @ and x 
Details of 


[Eq. (4.9)] is also evaluated in two parts. 
the calculation are discussed later in connection with 


The critical point of Eq. (4.11) occurs when w, 


V 2/(y + 11), or AL, 


1.0, and, by analogy with the 
exact solution [Section (3)], the integral curve in Re 
gion 2 is made to pass through this point, although it 
lies outside the region of real flows. The two regions 
are then joined by requiring that «x, equal «x; at the 
junction point. Also the dead-water region is assumed 
to vanish at this point, and therefore 6 coincides with 
5). 


Eq. (4.11) is integrated by two different methods: 


(1) When M < 


small that the equation can be integrated by an itera 


2.0 (approximately), dw,/dw, is so 


tion scheme in which dw,/dw, = 0 in first approxima 
tion. By solving the equation obtained by setting the 
numerator of the term on the right-hand side of Eq. 
(4.11) equal to zero, one finds w, and therefore dw, dw, 
as functions of w,;. The value of dw,/dw, so derived is 
then introduced on the left-hand side of Eq. (4.11), and 
the resulting equation is again solved for w, = f(w,), 
etc. 

(2) When ‘7 > 2.0 (approximately), the first scheme 
did not converge sufficiently fast. The integral curve 
is then obtained by passing a parabola with the correct 
slope and curvature through the critical point at w, = 
V 2/(y + 1) and utilizing numerical integration there 
after. This method is similar to the scheme utilized 
for the ‘‘exact’’ equation [Eq. (2.30) ]. 

The result of a typical integration for the flow in 
Regions (1) and (2) is shown in the w,-w, plane in 
Fig. 11 for Wf = 
The dashed curve represents the more exact solution 


3.0 and a turbulent wake (k = 0.03). 


obtained by the methods described in Section (3). 
(Clearly, for the curve passing through the point w, = 


22 CAL 
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the application of the ‘‘simplified’’ theory to the base 


Section (5) ]. 





pressure problem 
Region (2).—In Section (2.5) it was pointed out that a 
critical point of Eq. (2.30) does not exist when « 
constant, as in Region (1 The condition determin; 
the integral curve must therefore be associated with 


Region (2). In this region, x; S «x S 1 and f = 1, or 


» 9) | 
=r we | | 
y+ | 7) 
») w | 
WeV ¥Y | 
WV ¥ J 
V2, (y + 1), only the segment on whichx; S «x S< 1.0 


has any reality. 
Values obtained for Op at ko kK 0.720 are tabu 


lated as follows 


Mf On (Deg 
1.5 0.38 
2.0 0.54 
3.0 —1.24 
1 0 3.38 


The value of w, at the junction point does not differ 
appreciably from @,, so that the character of the solu 
tion in Region (1) is not very sensitive to small errors 
in the value of w, as determined by this approximate 


method. 


4.2) Characteristics of Separated and Reattaching Flows 
and Recompression in Wakes (‘‘Simplified’’ Theory 


Reattaching Flows.—Once it is known that @z and the 
recompression in Region (2) are small, the main charac- 
teristics of the recompression in wakes, or of reattaching 
flows in general, can be predicted from the form of the 
momentum equation and the flow solution in Region 
and Eqs. (4.6) and (4.8)]. By 


referring to the momentum equation [Eq. (2.9)] in a 


1) [Eqs. (2.9) or (4.2), 
jetlike flow (7, = 0, x = «;), one sees that the value of 
the pressure rise Ap times a mean 6 is equal to the 
momentum transported to the dissipative flow, fu, dm, 
which is of the order of magnitude of #, Am, minus the 
net change in momentum of the dissipative flow, which 
is of the order of magnitude xi, Am + mx,Au,. Now 
Au, ~ Ap/p,ti,, and Am is not much larger in order of 
magnitude than mm itself, or 


Am = p,u,6;[1 — (6*/6)], 
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By making use of these approximate relations one finds 
that (roughly) a eye — y | 
Ap /p.u,.* =~ 6; 6 
When 6 > 6;, as in the more familiar case of a “half } 
jet,’ then the pressure is practically constant at the #=0.720—- 
ambient value. But when two gas streams converge 19, ig | 
or when a gas ‘stream is directed toward a surface, Ss : 
then 6 — 6;, and large increases in static pressure are | & =—€=1.0 
made possible by the transport of mass (and momen- We K * 
tum) to the internal dissipative flow. This mechanism >" | 
is not confined to supersonic flows; it is undoubtedly > ~f] 
the main factor in the rapid pressure rise observed on an 18H \ wt a W > 
airfoil in low-speed flow when the separated laminar | REGION 2 we 
flow reattaches itself as a turbulent region downstream 
of the “‘separation bubble”’ (see, for example, reference | W,= We 
21 Incidentally, for a given mixing rate one would 
expect a much larger increase in pressure ratio at high 
velocities than at low velocities, because a 12 14 ie 18 20 
Ap Ap \ put? P en 
i" (- ©) * M,’ FIGURE 11. RECOMPRESSION IN TURBULENT 
WAKE. M=3.0, k=.03 
This conjecture is borne out by the numerical calcula 
tions to be discussed later and by experimental evidence. a 
(ualitatively the fori of the pressure distribution in 
reattaching flows or wake flows can be predicted from 7 
the relations (4.6)—(4.8). It is clear that dé/dx, and ™ 8 ” 
therefore (1/p)(dp/dx), is small when 6 > 6; and 
reaches its maximum value rapidly when 6 — dz be- 30. M=30 6. 
cause of the exponential character of the 6-6 relation. : } ~ Mei5 
The variation of 6/ér with 6 for Mlach Numbers of 1.5, % 20 _ 4 
2.0, and 3.0 is illustrated for a turbulent wake (k = tie i ~_M=2.0 
0.03) in Fig. 12a, while Fig. 12b shows the corresponding lO ° 
axial static pressure distribution for these cases. A esti 2 M=3.0 
typical case with transition in the wake 1s shown in 7 = =20 -30 ¥ “20-40 -60 -80 
Fig. 13a, and the corresponding pressure distribution 1s 6 X-Xa 
plotted in Fig. 13b. Here W = 3.0, Riury, = 0.03, and FIGURE 12. VARIATION OF THICKNESS B. DISSIPATIVE 
Re. = 2 X 10% For simplicity kjgm. has been taken FLOW REGION AND THE STATIC PRESSURE 
constant [see discussion in Section (2.3)] and equal to ALONG THE AXIS IN TURBULENT WAKE 
the flat plate value 0.00124 at Re. = 2 K 108. (These k =.03 SIMPLIFIED THEORY 
numerical examples are taken from the base pressure 
calculations to be discussed later. ~ a 
From Eqs. (4.6a) and (4.6b) and Figs. 12 and 13, it is 10, 10 
clear that 6 de increases rapidly with |@|. As pre (a) (b) 
dicted from the qualitative discussion of reattaching 8|. 8|) 
flows, most of the compression occurs near reattach 
ment, in the range 1 < 6 dg < 5, say. Im fact, by 6} LAMINAR -+ * 
Eqs. (4.6a) and (4.6b), e . TURBULENT 
5/5p =e (4.6c) . 
roughly, where 2 TURBULENT; . oA 
Or 
h=2VIP-—1(1 —« * 5 ee ee eee ee ee 
eS X-Xr 
and \/ is some average Mach Number. In the base FIGURE 13. VARIATION OF THICKNESS OF DISSIPATIVE FLOW 


REGION AND STATIC PRESSURE ALONG AXIS. 


pressure problem, for example, when & is constant, the 
TRANSITION IN WAKE 


base pressure ratio drops as h ‘6, or 6 dz increases but at 


a decreasing rate (high Reynolds Nuinber régime). M=3.0, kyype*.03, Re,=2 X 10%, kyay?.00124, 
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FIGURE 14. VARIATION OF THICKNESS OF 
SEPARATED LAMINAR FLOW REGION 


AND STATIC PRESSURE ALONG SURFACE. 
M=2.0, Rey =110,000, ky 4y=.0055, 
@,=1.23° 


Several important conclusions can be drawn from 
the approximate relation above [Eq. (4.6c) |. 

(1) For a specified value of 6 6g one finds that @ ~ 
Vk, or the compression, as measured by the isentropic 
flow deflection, increases with increase in mixing rate. 
Since the pressure ratio is a rather sensitive function of 
9, the marked difference between laminar and turbulent 
reattaching flows is now placed on a quantitative basis. 
Figs. 13a and 13b bear out the conclusion that the lam- 
inar portion of a flow with transition generally con 


tributes little to the recompression. 


(2) Since @ ~ (.\/? — 1)‘, the pressure ratio p p 
in a reattaching flow increases rapidly with Mach 
Number [Fig. 12b; see also Section (5.2) |. 


(3) @ is sensitive to the value of «x; assumed in the 
‘‘jet."" The factor (1 — x;) is a measure of the rate at 
which the internal dissipative flow ‘“‘destroys’’ momen- 
tum received from the external flow. Experimental 
data are required to determine this parameter, which in 
this paper is taken to be 0.720—1.e., the value for a 
‘free’ jet at constant pressure. 

Separated Flows.--The approximations introduced 
in Region (1) of the wake behind a supersonic airfoil 
with a blunt trailing edge [Section (4.1) ] are also appli- 
cable to separated flows with a supersonic external 
stream, where the influence of the viscous stress at the 
surface is negligibly small. An important example of 
such a flow occurs in the laminar boundary-layer- 
shock-wave interaction, where the laminar layer sepa- 
rates upstream of the incident shock. 

In the separated flow region, 


where 
(4.12) 
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and where 6, is the thickness of the dissipative flow and 
§, is the flow inclination ‘‘at separation.”’ [Of course, 
the value of @, cannot be predicted by the simplified 
mixing theory; it must be obtained by an extension of 
the usual boundary-layer theory’ or by integration of 
Eqs. (2.28) and (2.29) with o ¥ 0, & variable, and suit 
able boundary conditions. This problem will be dis 
cussed in a later paper.| By referring to the expression 
for d@ dx |Eq. (4.8)], one sees that (1, p)(dp/dx) has 
its maximum value in this case at the separation point 
and drops off rapidly downstream because of the ex 
ponential character of the 6-4 relation [Equation (4.12 
An example of such a laminar separated flow is illus 
trated in Figs. 14a and |4b. 

Here, Re,, = 1.10 X 10°, 6, = 1.23°,T IT = 2.0, 
Riam. = 0.0055.T A pressure distribution of this form 
is typical of laminar separated or ‘‘detached”’ flows 
and also of turbulent separated flows, although, in the 
latter case, because of the much larger mixing rate, the 
pressure increase is correspondingly larger. A pressure 
distribution similar to that shown in Fig. 14b 1s found 
downstream of flow detachment in an overexpanded 
supersonic rocket nozzle.’ 

In the case of the laminar boundary-layer-shock 
wave interaction, the region upstream of the incident 
shock is a separated flow, while the flow downstream 
of the reflected expansion fan is a reattaching flow 
Thus the S-shaped pressure distribution on the surface 
observed experimentally is understandable. Provided 
that friction is neglected, which is strictly correct only 
if the plate is cut off at the beginning of Region 2, the 
reattaching flow can be treated exactly as in the base 
pressure problem. Therefore, the laminar boundary 
layer-shock-wave interaction problem can be solved, at 
least approximately, by means of the simplified mixing 
theory, provided a suitable reflection condition for the 
incident shock wave is introduced. If the incident 
shock wave is reflected from the ‘‘jet’’ 
boundary, then the pressure drop across the reflected 


as from a free 
expansion fan exactly cancels the pressure rise across the 
incident wave. In that case 6. = 6; — 26,,, where 4; 1s 
the flow inclination angle in the separated region just 
upstream of the incident shock, #2 is the flow inclination 
angle in the reattaching flow just downstream of the 
reflected expansion fan, and 6, is the flow deflection 
across the incident shock. The other boundary condi 
tion at the junction of the separated and reattaching 
flow regions 1s, of course, 6; = 6. [6g for the reattach 
ing flow is determined from .\/,.., see Section (4.1).| 
In Fig. 15 the calculated static pressure distribution on 
the wall is shown for JJ = M = 2.0, 6, = 4 
Re., = 2.5 X 10°, ky, = 0:0037. 
pressure increase on both sides of the shock is about 


equal, and the flow pattern is nearly symmetrical. The 


For this case the 


at separation were computed by the 


approximate method of reference 4. The value of Rjgm. is taken 
at the flat-plate value for a Reynolds Number 


+ The flow conditions 


to be constant 


equal to Rex, 
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DISSIPATIVE FLOWS AN 
asyinmetric flow pattern and pressure distribution 
found at higher Reynolds Numbers, or for stronger in- 
cident shocks, appears to be associated with transition 
in the reattaching flow, which greatly increases the 
This flow 


régime can be treated by methods similar to those em 


mixing rate and pressure rise in this region. 


ployed in the base pressure problem for the case when 
transition occurs in the wake [Section (5)]. A detailed 
treatinent is reserved for a later paper. 

For the turbulent boundary layer the situation is 
more complicated. In this case, the deep penetration 
basic assumptions of the present paper at least locally, 
particularly the assumption of constant pressure across 
the internal dissipative flow [see Section (6) ]. 

5) APPLICATION OF 
THE BASE PRESSURE 


SIMPLIFIED MIXING THEORY TO 
PROBLEM IN PLANE SUPERSONIC 
FLow 


5.1) Method and Assumptions 


This problem has already been formulated in Section 
When 
the wake is fully turbulent, the procedure for the simpli- 


3) as an application of the more exact theory. 


fied theory is similar to that previously described in 
Section (3), but when transition occurs in the wake, the 
scheme is somewhat modified. 

(a) Turbulent Wake—Transition in Airfoil Boundary 
Layer.—-When the wake is fully turbulent, it is again 
assumned that k = constant, and 6/ 6x is then completely 
determined as a function of 6 for every Mach Number 
1 by Eqs. (4.6) and (4.7). 


the boundary condition 6 = (h 


At the airfoil trailing edge, 
2) + 6, must be satis 


fied, or 


(5.1 ) 
Approximately, x; = 0.720 and {1 — (6* 6)], = 0.45 
for a turbulent low-speed isothermal jet. To every 


value of 


[(h/2 


) + by, On: h 6, 
: ae eo 
Or Orn 26, by: 


there corresponds a value of 6 dr and therefore a value 
f 6, by Eq. (5.1) and, finally, a value of base pressure 
ratio py / pb = f(, M1). In order to relate this value 
of the base pressure ratio to 6, h, it is necessary to cal 
culate 6,-/dp and 6,-/6,. The ratio é6p/6,- is obtained as 
follows: 
By definition, 
Mr = (pr Ger) or] — (§ 
Mityy = (Por Gey)dy[1 — (6*/6) ], 


so that 
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FIGURE 15. STATIC PRESSURE DISTRIBUTION ALONG 


SURFACE FOR LAMINAR BOUNDARY 
LAYER —SHOCK WAVE _ INTERACTION. 
M.o*M=2.0, @e=4, Re,=250,000 
k, a= .0037 
SIMPLIFIED THEORY 


Mr (Pr Cer) Or 
My: (Por Gen) Or 


where it 1s assumed, as in Section (3), that 


6* 6* 6* 
; = ( -) = ( : - ) (1 _ ~ | = const. 
0 ne 0 R 0 P 


By integrating Eq. (2.25) for the mixing rate and divid- 


ing through by m,,, one obtains 





os 
Mr k en 
—=1+— — dx 
1» My t.F 7b" g, 
k "” p (db dé 
l+— — ——_—— (0 
Myr J’ og, (db/dx) 
or 
Mr Pr/ Ser Or) | 
— ] - is - cx som jaar a { D 
mM, Por Sqn’ byt [1 — (8* 6)] l—k« 
(5.3 
where 
“on / ¥\ 
0 (Pp ¢&, dg 
D = | (=) —————- — 
Be Or] (Pr/ Ger) VM? —-1 
Finally, by Eqs. (5.2) and (5.3), 
(Po’ Pep’) 


#(1 Ss | — D ) 
ss Ul — (6*/8)], Si — KJ 


For the completely turbulent wake, the integral D is a 
unique function of 6,, for every Mach Number /, ac 


( Pr ¥ Rr) 
(0. 


i 


cording to Eqs. (4.6) and (4.7). 

As explained in Section (3), the ratio 6,/6,, as given 
by Eg. (3.2b) is obtained by applying the condition of 
conservation of mass across the expansion around the 
‘“‘corner’’ at the airfoil trailing edge. 

By utilizing Eqs. (5.1), (5.4), and (3.2b), the base 


pressure ratio is related to 6,/h = (6,/c)(c/h). The 
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variation of base pressure with Reynolds Number at a 
given Mach Number is obtained by introducing the 
same postulated variation of 6,/c with Reynolds Num 
ber and Mach Number which is described in Section 
(3) and illustrated in Fig. 8 for WV = 38.0. 

(b) Transition in Wake—-Airfoil Boundary Layer 
Laminar.—-When transition, in the airfoil boundary 
layer is described in terms of curve (1) of Fig. 8, for 
example, then for Re, | X 10° the boundary layer is 
completely laminar, and the transition ‘point’? moves 
off the body and into the wake. Since the mixing coef 
ficient k in the laminar portion of the wake is a fune 
tion of Reynolds Number, the relation between 6/5, 
and @ is no longer independent of Reynolds Number, 
and the procedure outlined in (2a) above must be modi 
fied. 

In the first place, a statement must be made regard 
ing the transition Reynolds Number in the wake, about 
which little is known at present. Transition in the 
wake must depend to a large extent on the previous his 
tory of the laminar oscillations in the airfoil boundary 
layer, as well as on the growth of disturbances along the 
laminar portion of the wake itself. Based on the 
meager data available, it is assumed that a ‘‘free’’ high- 
velocity laminar half-jet, undisturbed by oscillations 
fed into the jet from an upstream boundary layer, 
would have a transition Reynolds Number of 4 X 10°. 
Therefore, when Re, < 4 X 10°, the transition Reyn 
olds Number in the wake must approach a value of 
4 X 10°, while Re,,,,. = 0 when Re. = 4 X 10° and 
transition just reaches the airfoil trailing edge. Taking 
a linear variation of Re;,,,,,. with Re, between these two 


conditions, roughly. 


Rew ware 1 X 10° — Re,/10 (5.5) 


[For transition curve (II), Reér,,,. 
= 3 X 10° (Re./10).]| A smooth movement of 
transition from airfoil to wake is presumed. All that 
can be claimed for this approximation is that it gives 
the correct form of the dependence of transition on 
Re, and provides for an accelerated forward movement 
of transition in the wake with increasing Reynolds 
Number. 

On the basis of this assumption regarding transition 
in the wake, the variation of base pressure with Reyn- 
olds Number at a given Mach Number in this flow 
régime is calculated by an iteration method. Suppose 
that c/h is given. For some value of Re, < 4 X 10° 
a series of trial values of ,, or the flow deflection angle 
at the transition point, is selected. The corresponding 
locations of the transition ‘‘point’’ and values of 
Rew... are calculated by an iteration scheme, in which 
the value of @,, and the base pressure are also calcu- 
lated. Now the correct value of Rey... at the pre- 
scribed value of Re, is given by Eq. (5.5), so the proper 
value of 6, can be determined from a plot of Re... VS. 
6,, and the corresponding base pressure can be calcu- 
lated by interpolation. 


at high velocities. 
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The details of the method are as follows: Once y 
trial value of @, is selected, the relation between by 
and @ required in the statement of the boundary coy 
dition [Eq. (5.1) ] is obtained from Eqs. (4.6) and ( 1.10), 
However, the value of the mixing coefficient in the 
laminar portion of the wake must be known or chosen 
For our present purposes, it is sufficient to take &,,,, 
equal to the mixing coefficient in the laminar boundary 
layer just at the airfoil trailing edge —i.e., 7 


P c;,/2 g(M) 0.332 
fins l1—« WRe1 — x 


Che Mach Number factor g(.\/) is obtained from Fig 
of Van Driest’s paper’ and is tabulated below 


VW g( MJ) 
1.55 0.98 
20 0.95 
3.0 0.9] 
$0 0. &7 
5.0 0.8] 


The value of xo utilized in calculating k),,. 1s 0.827, cor 
responding to u(6)/u, 0.9975, and the corresponding 
value of x; employed is 0.720.5 <A. typical plot of 
6/ be vs. 6 with transition in the wake is illustrated in 
Fig. Ida. 

The right-hand side of Eq. (5.1) is also determined 
as a function of @ by utilizing Eqs. (5.4) and (3.2b), 
where the value cf [1 (6*/6) |, for the laminar bound 
ary layer with zero heat transfer at the airfoil surface 
is given byt 


3.87/[5.60 + 1.20(y7 1)M,?] 


The value of [1 (6*/5) ],’ is taken to be 0.45, exactly 
as in the case of the fully turbulent jet. 
5y/h = (6)/c)(c/h) is already known once c/h and Re, 
The cor 


The value of 


are prescribed, since 6,/c is given by Eq. (3.5). 
rect value of 6, for the particular trial value of 4, is 
then found from the intersection of the curves of 6/6 
and the right-hand side of Eq. (5.1) plotted as functions 
of 6. The location of the transition point in the wake 
aft of the trailing edge is then calculated from the rela 


_— oe by/bp ee 
~ = = = d (0.0 
OR OR  (5/5R)tr o + Riam OR 


and the corresponding value of Rej;,,,, is given by the 


tion 


l Pu *s M, -- 
Reso... = - Re, F 7 (5.7 


relation 


where 


if u~ 7”, and 


tIfu~ 7, g(M) = 1.0 


t See footnotes on page 6638 
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DISSIPATIVE FLOWS AND 


/ l 6* Or Bn: Oy 
—s | k,l | ae 
( OR 0 j On On 
0.64. 
from Eq. (5.7) 


For air, ¥ 1.40 and m 0.76, so n 

Since the value of Rey,,,, obtained 
does not, in general, coincide with the assumed value 
given by Eq. (5.5), the whole procedure must be re 
peated for three or four trial values of @,, in order to 
determine the correct Re,,,,, and base pressure by 
interpolation with sufficient accuracy. The scheme is 
repeated for several representative values of Reynolds 
Number, Re,. Fortunately, after some experience only 
a few trials are necessary. 

Evidently, for sufficiently low Reynolds Numbers, 
§,, > 92 and the flow in Region (1) becomes completely 
laminar. For Reynolds Numbers below the value at 
which 4, is calculated to be approximately equal to 4p, 
the conditions are determined by analogy with the more 
exact theory. In this theory, conditions upstream of 
the critical point are unaffected by conditions down 
stream of this point and are therefore independent of 
the movement of transition once transition passes the 
critical point. By analogy, it is assumed now that, 
once transition passes into Region (2), conditions in 
Region (1) remain independent of Region (2). The 
value of 6, at a given Mach Number J/ is regarded as 
a constant and the character of the flow solution in 
Region (1) is determined completely by the variation 
of laminar mixing rate with Reynolds Number {see 
Section (5.2)]. At extremely low Reynolds Numbers, 
the present mixing theory must be modified because 
4,/h is no longer small. 


5.2) Effect of Reynolds Number on Base Pressure: The 


Various Flow Regimes 


A sufficient number of nuinerical examples have been 
calculated for all flow régimes to indicate the important 
trends in the variation of base pressure with Reynolds 
number, Mach Number, c/h, and type of transition in 
the airfoil boundary layer. Fig. 16 shows the variation 
at Wf = 3.00f py: /p with Reynolds Number for c/h = 
10.0, with the two types of assumed boundary-layer 
[In all the other figures only the calculated 
In Fig. 9 the 


transition. 
base pressure for curve (I) is shown. | 
curve computed from the simplified theory for 17 = 
3.0 and c/h - 
obtained by the 


3.0 is shown in comparison with the curve 
“exact’’ theory [Section (3)]. These 
curves show a similar qualitative behavior. In view of 
the assumptions made in the simplified theory and the 
probable sensitivity of the results of the ‘‘exact’’ theory 
to the assumed F-« curve [Sections (2) and (3)], the dif 
ferences in the numerical values of base pressure should 
not be taken too seriously at present. The numerical 
results are properly regarded as preliminary. 

The effect of Reynolds Number and c/h on base 
pressure at 17 = 2.0 is illustrated in Fig. 17 for values 
of c/h of 3.0 and 10.0. The low Reynolds Number 
region at J7 = 


2.0 is emphasized in Fig. 18 for c/h = 
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9 Sections (3) and (4.2)]. These qualitative considera old 
__.CHAPMAN tions are borne out by the results of the computations che 
8 arr THEORY In order of increasing Reynolds Number, the fou; ath 
: flow régimes are as follows: cat 
7| ' (1) At sufficiently low Reynolds Numbers, not onh ig 
M=15 the body boundary layer is laminar but also that region Ice 

6 4 of the wake in which the major portion of the recom 
a . 3 pression occurs [Region (1)].. Within this régime, as tio 
Re, increases, the laminar mixing rate decreases and tra 
4 % a M-2.0 the base pressure increases slowly with Reynolds Num stil 
\ ~ ; - ber (Fig. 18). The decrease in Rjam with increasing litt 
3 a —s Re, is apparently more than enough to offset the ‘oil 

- crease in thickness of the boundary layer at the airfoj] 
2 . " trailing edge or the decrease in thickness of the mixing ni 
M=3 | region in the wake. Because of the low laminar mixing rég 
a 2 4 6 86 10 1t2 14 16 18 20 Fate, the base pressure ratio is relatively high str 
Re. x % (2) As the Reynolds Number is increased, eventually lec 
FIGURE 19. BASE PRESSURE IN FULLY TURBULENT transition begins to move upstream in the wake, and wil 
REGION the accompanying large increase in local mixing rat thi 
(Rearp, ~ D-10 Riam ) counteracts the effects discussed in Th 
10.0, and in Fig. 20 the base pressure with airfoil bound- (1) above. After reaching a local maximum (Fig. 18), the 
ary layer laminar to the trailing edge is plotted against © the base pressure drops precipitously with increase in Fi 

the parameter Re, *h/c employed by Chapman.” Re, (Figs. 8, 16, 17, and 20). The large increase in 
Data for two of the airfoils with blunt trailing edges local mixing rate is a much more important factor than rea 
tested by Chapman in this flow régime are also shown — the accompanying increase in the thickness of the mix is ¢ 
for comparison. ing region, which, by itself, would have just the opposite slo 
In the ‘fully turbulent’’ régimes at high Reynolds _ effect on base pressure. ber 
Numbers, the base pressure ratio is plotted against After transition in the wake has moved close to the lay 
the parameter Re, ‘h/c for M = 1.5, 2.0, and 3.0 in airfoil trailing edge, the base pressure ratio continues ba: 
Fig. 19. Again Chapman's data in this flow régime to drop with increasing Reynolds Number, so long as Re 
are shown for comparison; the curves are reproduced _ the airfoil boundary layer remains laminar, because of ing 
from Fig. 28 of reference 17, in which the ‘“‘mean’”’ the decrease in boundary-layer thickness at the trailing Re 
curves for the entire “thickness group” of airfoils are edge with increasing Reynolds Number. This de 
presented. The thickness ratios for the group are 0.05, crease is expected to be much less rapid than the first tra 
0.075, and 0.10, and the ratio of airfoil trailing-edge steep drop that occurs when transition moves forward are 

thickness to maximum thickness has the discrete values jin the wake. If the boundary-layer transition Reyn 
0.25, 0.50, 0.75, and 1.0, so that c/h ranges from 10.0 tra 
to 80.0 in these tests. cola: wus tauie the 
It is remarkable that, in spite of the approximate SIMPLIFIED THEORY cul 
character of the present mixing theory, it appears to Bi at 
give the correct fluid-mechanical explanation of the aff 
observed phenomena and flow régimes in the base pres- 7 sur 
sure problem. Unfortunately, Chapman's data for \ irre 
the blunt trailing edge airfoils do not include the im- P. 6 4, = 3.0 len 
portant intermediate régime, in which transition moves r is € 
upstream in the airfoil boundary layer, but base pres- 5} " 710 old 
sure data on bodies of revolution,® *° which should be 426.7 pre 
at least qualitatively similar to the two-dimensional 4 fee 
case, illustrate perfectly the régimes determined by our vee pre 
calculations. Four such flow régimes can be distin 3} thi 
guished, in which the contrasting effects of mixing rate air 
variations and variations of 6, determine the quite dif 50 100 50200 250 300 ( 
ferent dependence of the base pressure on Reynolds Re % x hy eX] 
Number in the different flow régimes. In preceding rat 

: . : : . : FIGURE 20 TRANSITION IN WAKE, 

oe it has been shown that an increase in k ed AIRFOIL BOUNDARY LAYER Mi 
creases the pressure ratio across a compression an | HOW! NCE ray 
therefore decreases the base pressure [Section (4.2) ], aa Posse a Th 
while an increase in 6, increases the base pressure M=2.0 ph 
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DISSIPATIVE FLOWS 


olds Number is sufficiently high and if the ratio of 
chord to airfoil trailing edge thickness is low, but nol 
atherwise, the base pressure ratio depends only on the 
ratio 6, or on the parameter Re,*h/c over an ap 
preciable Reynolds Number range in this flow régime 
see later discussion, Section (5.2) (72) |. 

3) At a sufficiently high Reynolds Number, transi 
tion will occur in the airfoil boundary layer just at the 
trailing edge. When the Reynolds Number is increased 
still further, presumably the mixing rate in the wake is 
little affected, but transition moves forward on the air 
foil, and the boundary layer at the trailing edge thick 
ens accordingly. Thus, the base pressure ratio at first 
increases With increasing Reynolds Number in this flow 
régime. However, when transition has moved up 
stream sufficiently far from the trailing edge, the normal 
lecrease in local turbulent boundary-layer thickness 
with increasing Reynolds Number begins to offset the 
thickening effect of the forward movement of transition. 
The base pressure curve reaches a local maximuin and 
then decreases again with increasing Reynolds Number 
Figs. 8, 16, and 17). 

1) Finally, a high 
reached in which transition in the airfoil boundary layer 
is essentially “‘fixed’’ and the base pressure ratio drops 


Reynolds Number régime is 


slowly but noticeably, with increasing Reynolds Num- 
ber, because of the decrease in turbulent boundary- 
layer thickness at the trailing edge. In this régime, the 
base pressure ratio is a function solely of the parameter 
Re, *h/c (Fig. 19) or of some similar parameter involv- 
ing the product of h/c and a logarithmic function of 
Reynolds Number. 

Briefly, the effects of Mach Number, c/h-ratio and 
transition in the airfoil boundary layer on base pressure 
are as follows: 

t) Transition and Mach Number.—-Whatever affects 
transition in the airfoil boundary layer also determines 
the location, shape, and depth of the ‘“‘valley”’ in the 
curve of base pressure ratio versus Reynolds Number 
at a given Mach Number. Thus, the base pressure is 
affected by the heat-transfer rate at the airfoil surface, 
surface pressure gradient, surface roughness, surface 
irregularities such as sharp corners, free-stream turbu- 
lence,etc. Since transition in the airfoil boundary layer 
isexpected also to depend on Mach Number, the Reyn- 
“‘valley’’ in the base 
The ef- 


fect of each of these factors on base pressure can now be 


olds Number range defining the 
pressure curve is a function of Mach Number. 


predicted, at least qualitatively, when its effect on the 
thickness and character of the boundary layer at the 
airfoil trailing edge is known. 

On the basis of the discussion in Section (4.2), it was 
expected that, for a given mixing rate, the pressure 
ratio supported by the wake flow increases rapidly with 
Mach Number 
rapidly with increasing Mach Number (see Fig. 19). 


i.e., the base pressure ratio decreases 


This result is probably also significant for other flow 
phenomena in which mixing is important, such as 


AND 
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boundary-layer separation. Of course, the actual 
numerical values of base pressure depend critically 
upon the mixing coefficient k. By selecting & as a suit 
able function of Mach Number, coincidence between 
the theoretical computations and experimental data 
could be obtained. By referring to Fig. 19 one may 
conclude tentatively that the mixing rate at low super 
sonic Mach Numbers is probably of the order of 0.05 
(somewhat larger than the constant value of 0.03 as 
sumed in this paper) and that k apparently decreases 
with increasing Mach Number. Some support for this 
tentative conclusion comes from recent experimental 
work on the rate of spreading of a supersonic turbulent 
stream into quiescent air 

(1t) Chord to Trailing-Edge Thickness Ratio (Trans 
ition in Wake).—-In general, the higher the ratio of 
chord to trailing-edge thickness, the larger is 6,/4 and 
the higher the base pressure ratio at a given Reynolds 
Number and a fixed Mach Number.¢ The effect of 
c/h ratio is particularly interesting in the intermediate 
wake transition region, because the rate of upstream 
movement of transition in the wake with increasing 
Reynolds Number is governed largely by the parameter 
c/h. The transition Reynolds Number Re,,,,,, is de 
termined by Re, [Eq. (5.5)], and therefore so is the 
ratio //c = Re, ware’ Re,, where lis the distance from air 
But the 
important parameter governing the character of the 
When the transition “‘point’ 


foil trailing edge to wake transition ‘“‘point.”’ 


wake flow is the ratio //h. 
just enters the significant recompression region, 

= a DS 
roughly. Two important conclusions can now be 
drawn: 

(a) At a given Reynolds Number Re,, the laminar 
flow region occupies a larger portion of the wake the 
larger the c/h ratio. In other words, the forward 
movement of transition into the major recompression 
region and the accompanying drop in base pressure 
begins at higher Reynolds Numbers the higher the 
c/hratio. For very large c/h the base pressure ratio is 
nearly unity, and 6, is small and varies nearly linearly 
with 4 c. Therefore, by Eq. (5.8), //c approaches a 
limit, and the Reynolds Number marking the transition 
between the two flow régimes also approaches an upper 
limit as ch increases. These conclusions seem to be 
borne out by the data presented in Figs. 14 and 15 of 
reference 17. 

(b) For large c/h, a significant variation in /// occurs 
with only a small variation in //¢ and, therefore, in Re 
Once the transition point does begin to move forward 
into the major recompression region, the base pressure 
drops steeply in a narrow Reynolds Number range. On 


t Thus the effect of increasing body length to obtain higher 
Reynolds Numbers in a series of experiments is certainly much 
different from the effect obtained when c/h is fixed and ambient 


density is increased rhis fact was pointed out by Bogdonoff.** 
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the other hand, for small c/h, the upstream movement 
of transition in the wake extends over a wider Reynolds 
Number range, and the curve of base pressure vs. Re, 
is much ‘‘flatter’’ in this flow régime (Figs. 17 and 20). 
Since this effect of c/h ratio is greatly accentuated at 
extremely large values of c/h, it may be responsible for 
the almost discontinuous change in base pressure ob 
served by Chapman" under certain conditions for air- 
foils with extremely thin trailing edges (c/h = 80). 

In the flow régime characterized by the forward 
movement of transition in the wake with increasing 
Reynolds Number, the base pressure evidently depends 
not only upon 6,/ but also upon the location of transi- 
tion. In géneral, therefore, one would not expect the 
base pressure to be correlated solely against the param- 
eter 6,/h, or Re,’h/c in this régime, but to depend 
also onc/h. Sucha correlation exists only when transi- 
tion in the wake is located close to the airfoil trailing 
edge over a significant range of Reynolds Numbers be- 
low the airfoil boundary-layer transition Reynolds 
Number. By referring to Eq. (5.8), one sees that this 
condition is fulfilled only when: (1) the boundary- 
layer transition Reynolds Number is high and the value 
of Re,,.... Varies only s'owly with Re, (//c varies slowly 
with Re,); In Fig. 20, neither the 
theoretically calculated base pressure nor Chapinan’s 
data are correlated by the parameter Re,’*h/c, largely 
because Re,, in the airfoil boundary layer is certainly 
not greater than 4 X 10° (probably Re,, ~ 2 * 10° in 
these experiments) and the values of ch are already 
At low Reynolds Numbers, the curves of 


(2) c/h is small. 


too large. 
base pressure ratio vs. Re,*h/c for various ¢/h values 
tend to converge and cross as shown in Fig. 20 and in 
Figs. 14 and 15 of reference 17, because the base pres- 
sure is approaching a local maximum of the type illus- 
However, at high Reynolds Num- 
For example, 


trated in Fig. 18. 
bers these curves diverge considerably. 
any attempt to construct a ‘“‘universal’ curve of base 
pressure ratio vs. Re. *h/c by utilizing data obtained 
at one Reynolds Number over a range of c/h values 
could lead to appreciable errors. 

In the case of blunt-based bodies of revolution,® 
transition in the body boundary layer apparently oc- 
curs at high Reynolds Numbers (Re,, > 6 X 10®), and 
the base pressure ratio may be correlated against the 
parameter Re,*D/L over a definite range of L/D 
values (L is body length and D is base diameter). 

The ability of the present theory to predict even such 
details of the base pressure phenomena is another con- 
firmation of the correctness of the basic mechanism as- 


sumed. 


(6) FUTURE INVESTIGATIONS 


The success of the mixing theory in furnishing a quali- 
tative fluid-mechanical explanation of the phenomena 
occurring in the base pressure problem encourages us to 
attempt to apply this theory to other problems involv 
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ing the interaction between a dissipative flow region 
and an adjacent nearly isentropic stream. However. 
in certain cases, either some of the basic assumptions 
of the present paper must be modified or the outer 
boundary of the dissipative flow region must be care 
fully redefined. In some cases, perhaps the most 
serious assumption is that the static pressure is con 
stant across the dissipative flow and equal to the loca] 
value p, in the external stream at y = 6. Ina turbulent 
boundary-layer-shock-wave interaction, for example, 
the incident oblique shock is reflected and refracted as 
it enters the layer, but apparently penetrates to the 
sonic velocity line close to the surface. The main re 
flected shock wave appears to emerge from this same 
region.'® 7? In the important interaction region, the 
static pressure is certainly not uniform across a plane 
transverse to the flow in the supersonic portion of the 
boundary layer, but the pressure is probably nearly con 
stant across the subsonic and transonic zones close to 
the wall. 
most strongly affected by an increasing static pressure 


Presumably, it is also this region that is 
along the surface. It might be correct physically to 
define this region as the new dissipative flow and to 
treat the supersonic portion of the boundary layer as an 
“external” inviscid rotational flow. Experimental and 
theoretical studies are required to settle this point. 
The mixing theory including friction |Eq. (2.30)] is 
expected to furnish information about the turbulent 
boundary layer in an adverse pressure gradient away 
from shock waves. It 
whether the maximum pressure increase that can be 


is important to determine 


supported by a turbulent layer without separation in 
creases rapidly with Mach Number, as we have found 
for the recompression in the wake flow. 


In many respects, laminar boundary-layer-shock 
wave interaction appears to be simpler than the turbu 
lent case, because upstream separation occurs even for 
weak incident shocks and the incident shock wave is re 
flected from the separated flow very nearly as if this 
flow were a free boundary. The static pressure is prob 
ably nearly constant across the low-velocity region be 
tween the jetlike outer boundary and the surface. The 
approximate solution indicated in the present paper 
[Section (4.2)] must be regarded as a first step; in a 
later paper, the important influence of surface friction 
on the reattaching flow and the surface pressure dis 
tribution downstream of the reflected expansion “‘fan 
will be discussed. 

Another important basic assumption of the present 
mixing theory is that the stagnation enthalpy is con 
stant across the flow. For separated flows and wakes, 
there is some question whether the mixing process ma) 
not create regions of low stagnation enthalpy near a 
surface or just aft of a body. This question is best 
settled by experiment. For flows attached to a surface, 
the present theory can probably be extended to cases 
where the surface heat transfer rate is not zero 
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DISSIPATIVE FLOWS 

The present theory is capable not only of giving a 
qualitative explanation of all these interaction phe 
nomena but should furnish realistic quantitative predic 
tions once the F-« curve is known and the dependence of 
mixing coefficient k on the parameters m, w,, and x is 
This information is required particularly 
By specifying the quantities 


determined. 
in the turbulent case. 
that must be known, the theory also serves as a useful 


guide to experimental investigations. 


(7) CONCLUSIONS 


|) An important property of supersonic wake flows 
or of reat‘taching supersonic flows directed toward a 
solid surface is the existence of a singularity in the 
basic differential equation, or ‘‘critical point”’ in the flow 
which acts much like the throat of a nozzle in deter 
mining the base pressure, for example. One important 
reason for the marked difference between laminar and 
turbulent interactions is the fact that the turbulent 
mixing rate is five to ten times larger than the laminar 
mixing rate. 

2) According to the present mixing theory, sepa 
rated flows, as well as reattaching flows, are capable of 
supporting considerable pressure increases at high 
velocities. The maximum compression, as measured 
by the isentropic flow deflection, is roughly propor 
tional to the square root of the mixing rate parameter 
and also to (M? — 1)‘, where MV is some average 
Mach Number. Thus, the pressure ratio increases 
rapidly with Mach Number, or the base pressure ratio, 
for example, decreases rapidly with Mach Number. In 
separated flows, the pressure gradient along the surface 
isa Maximum at separation and drops off steeply down- 
stream; in reattaching flows or wake flows, the pressure 
gradient is negligible some distance upstream of the 
“reattachment point’’ and increases rapidly toward a 
The inflected 
laminar 


maximum as this point is approached. 
surface pressure distribution observed in 
boundary-layer-shock-wave interactions is now under 
standable, and the distribution itself can be calculated 
approximately. 

(3) When the present mixing theory is applied to the 
problem of determining the base pressure for a super- 
sonic airfoil with a blunt trailing edge, it gives the cor- 
rect fluid-mechanical explanation of the observed 

Four different flow régimes can be dis- 
In order of increasing Reynolds Number 


phenomena. 
tinguished. 
they are as follows: 

i) A low Reynolds Number region in which not only 
the airfoil boundary layer is laminar but also the region 
of the wake in which the major portion of the recom- 
pression occurs. In this régime the base pressure ratio 
py P increases slowly with Reynolds Number Re, be 
cause of the decreasing laminar mixing rate. 

#) A régime in which transition to turbulent flow 
moves upstream in the wake and the base pressure curve 
reaches a local maximum and then drops precipitously 
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with increasing Reynolds Number because of the large 
increase in mixing rate. If the ratio of chord to airfoil 
trailing-edge thickness (c//) is low or if the boundary 
layer transition Reynolds Number is high (but not 
otherwise), an appreciable range of Reynolds Numbers 
exists for which transition in the wake is located quite 
near the base, but the airfoil boundary layer remains 
laminar. In this range, p,-/p continues to decrease 
with increasing Reynolds Number because of the de 
creasing boundary-layer thickness at the airfoil trailing 
edge, and the governing parameter is Re, h/c. 

(vit) A régime in which transition, having ‘‘jumped”’ 
from the wake to the airfoil boundary layer, moves for 
ward on the airfoil with increasing Reynolds Number. 
In this régime, p,-/p at first increases with increasing 
Reynolds Number because of the increase in boundary 
layer thickness at the airfoil trailing edge as the transi- 
tion point moves forward. Eventually, the normal de- 
crease in local turbulent boundary-layer thickness with 
increasing Reynolds Number begins to offset the thick 
ening effect of the forward movement of transition, and 
the base pressure curve reaches a local maximum and 
then begins to drop with increasing Reynolds Number. 

(iv) A high Reynolds Number régime in which transi- 
tion on the airfoil is practically “fixed” far from the 
trailing edge and the base pressure drops slowly with 
increasing Reynolds Number because of the decreasing 
turbulent boundary-layer thickness at the airfoil trail- 
ing edge. The governing parameter here is Re, ‘‘h/c. 

Clearly, anything that affects transition in the airfoil 
boundary layer or in the wake will affect the ‘‘valley”’ 
in the curve of base pressure versus Reynolds Number. 
The effect of a given parameter can be predicted, at 
least qualitatively, once its influence on the character 
and thickness of the airfoil boundary layer at the trail 
ing edge is known. 

Qualitative agreement is found between these theo- 
retical predictions and calculations and the base pres- 
sure data of Bogdonoff and Chapman on blunt-based 
bodies of revolution and Chapman’s data on supersonic 
airfoils with blunt trailing edges. 

(4) Considerable theoretical and experimental work 
remains to be done in extending this mixing theory to 
axially symmetric flows and to boundary-layer-shock- 
wave interactions and other flows, and in carefully de 
termining the dependence of the mixing rate and the 
mean temperature-mean velocity relation on the flow 


parameters. 
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A Method for Reducing the Analysis of 
Complex Redundant Structures to a 
Routine Procedure 


L. B. WEHLE, Jr.,* ann WARNER LANSING* 


Grumman clircraft Engineering Corporation 


(1) SUMMARY 


A general method based on minimum strain energy is given for 
the analysis of statically indeterminate structures. Matrix 
methods are used throughout. 

The procedure divides the work into three phases: (1) the writ 
ing of the equilibrium equations that determine the load distribu 
tion in the structure after it has been made statically determinate; 
2) the tabulation of the flexibilities of the individual members 
comprising the structure; and (3) the self-checking matrix opera- 
tions, which can be performed in a routine manner. In connec 
tion with phase (2), formulas and curves are provided for some of 
the elements used most frequently in idealizing aircraft struc 
tures. 

The method has been applied at the design stage to various 
components of several modern aircraft. For purposes of demon 
stration, a simplified example is worked out in this paper. The 
structure representing the root portion of a sweptback wing, hav 
ing taper in plan form and thickness and nonparallel tapered bulk 
heads, is idealized, and its load distribution and deformation 


characteristics are obtained 


(Il) INTRODUCTION 


eg MISSIONS OF high- and low-speed aircraft 
are constantly giving rise to new structural ar- 
rangements that must be designed to minimum weight. 
Some of these structures are highly indeterminate and 
of such an unfamiliar form that neither the commonly 
used short-cut methods nor previous test experience 
provides sufficiently accurate information for structural 
design. The method presented here was developed for 
this type of structure and reduces the analysis to a rou 
tine procedure. 

The basis for the procedure is the commonly used 
method of least work, which Levy! and Lang and Bis 
plinghoff? have recently demonstrated in connection 
with certain swept and unswept wings and which other 
authors* > have applied to similar problems. As _ usu- 
ally employed, this method is cumbersome and some 
what disorganized; however, if airplanes are to be 
completed on schedule, a procedure that is routine and 
yet can be applied to any irregular structure should be 
available to designers. The procedure given here is an 
attempt to fill that need, and it has been successfully 
applied to the design of the unique structures found in 


Presented at the Structures Session, Twentieth Annual Meet- 
ing, I.A.S., New York, January 28—February 1, 1952. 
* Structural Methods Engineer 


several sweptback wings, as well as to other prob- 
lems such as those involving rigid frames and thick 
plates. 

Simplicity is obtained in this method by separating 
the various operations required by any indeterminate 
analysis. This permits the intermediate results to be 
visualized while the calculations are in process and al 
lows most of the operations to be distributed to com 
puters for maximum efficiency. 

In this paper the theoretical development used to 
justify the recommended procedure is first given. 
There follows a summary of those formulas that are 
actually required in the computation, together with a 
discussion of the procedure to be followed. Thereafter, 
a collection is included of the elastic properties of a 
number of the structural elements which may be en- 
countered in an aircraft problem. In order to demon- 
strate the application to an actual structure, an ex 
ample is given of the analysis of a simplified swept two- 


cell box. 
(111) SymBots 
Pet 
Py = applied loads 

Jir J = member loads, any load acting directly upon an 
individual member; this includes any applied 
load or reaction acting directly upon a member 
and all redundant loads 

dre Ise Gt == redundant loads 

Yim» = unit load distribution due to applied loads when 
all redundant loads are set equal to zero; ith 
(jth) member load due to a unit mth (nth) ap 
plied load 

Vir, 7 = unit load distribution due to redundant loads 
when all applied loads are set equal to zero; 
ith (jth) member load due to a unit rth (sth) 
redundant load 

pee = rth (sth) redundant load due to a unit mth (nth) 
applied load 

Dim = unit load distribution in redundant structure; ith 
member load due to a unit wth applied load 

da = strain energy of ath structural member 

l = strain energy of entire structure 

ai flexibility influence coefficient for the ath mem 
ber; displacement at load g; due to a unit 
member load gj, where g; and g; include all 
loads acting on the member except those at as- 
sumed support points 

OG; = sum of ai; terms for all members 
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@mny Ams, 
Grn, Ary = flexibility influence coefficients for the statically 
determinate structure; 

Pm or gr due to a unit load P,, or q 


displacement at load 


) = flexibility influence coefficient for the redundant 
structure; displacement at load P,, due to a 


unit load P,, 


bai = displacement of ath member at load g; with re 
spect to assumed support points 

5m = displacement of redundant structure at load P, 

[ ] = a matrix 

{ } = a column matrix 

PF = the tranpose of a matrix 

{ ] . = the inverse of a matrix 


(IV) THEORY 


The purpose of this section is to justify the procedure, 
discussed in Section (V) for solving structural problems. 


(1) Statically Determinate Load Distribution 


The structure to be analyzed is idealized as consisting 
of a number of structural elements such as axially 
loaded bars, shear panels, beams, or any other conven- 
ient elements whose deformation characteristics are 
either known or can be easily obtained. The structure 
as a whole is acted upon by certain externally applied 
loads that are designated by P’s, each P having assigned 
to it a separate numerical subscript. For the purpose 
of the derivation, these subscripts are described only by 
the letters m,n, oru. The loads acting on the members 
are designated by q’s, each q having a different numeri 
cal subscript which is symbolized by zor7. The g’s may 
represent forces, shear flows, moments, etc., and as such 
may be considered as generalized forces. The extent 
to which the internal loads are broken down into sepa 
rate g’s will depend upon both equilibrium and energy 
considerations. This point is discussed more com 
pletely in the next section. 

With certain of the member loads considered as re 
dundant and therefore as having arbitrary magnitudes 
along with the applied loads, a balanced load distribu 
tion may be obtained by first writing the proper equa 
tions of equilibrium. These equations will be of the 


form 


cugi + T CuQi + Cumin + = (Q (1) 


Designating the numerical subscripts of the redundant 
member loads by r,.s, or ¢, these equations may be 
solved to yield expressions for the member loads in 
terms of the applied loads and redundants 


y= + Vim n + +- Yirdr + } 
(2a) 
= + VYimE'm + r YirGr + 
in which a trivial relation g; = g, or gi = P» is intro- 


duced wherever a member load is also a redundant 
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load or an applied load—that is, where the subscripts 
tandr or 7 and m are identical. This result may be 


written in matrix form as 
f = > 1f 
a = [Yim] {Pm T LYirJU9rs (2b 


In this equation the y’s may be regarded as load distri 
butions in the statically determinate structure due to 
unit applied loads or redundants. The numerical] cal 
culations will be more compact if the terms on the right: 
hand side of Eq. (2b) are combined. It may be r 


written as 
f ] 
Ge! = [VimiVir] | Pies (2c 


in which a dashed line indicates partitioning of a ma 


trix. 


(2) Strain Energy 


The strain energy to be minimized is collected from 
the individual members. Each member deforms under 
the loads g applied to it by adjoining members. For 
any member a, certain of these loads may be regarded 
as reactions; the rest, as causing deflection of the 
structure relative to these reactions. The relative de 
flection 6,; at the load g; depends on the loads applied 


at that point and at other points. 


5a = AqilQi 5 5 Qa i242 - = Qaijfj T (oa 


In matrix notation the set of equations is given by 


( 


f 9 
Sai} = | Qaij Vi I (3b 


It is seen that the a,’s describe the flexibility of member 
a. Thus, ag;;is the member displacement at q; relative 
to the chosen reactions and brought about by applying 
a unit value of the jth load. The energy in each mem 
ber is 


U, = (1/2)(g16a1 + G2da2 + + gibai + ) (da 


or 
U, = (1/2) [gil } Sai} (4b 


Substituting Eq. (8b) into Eq. (4b), the energy in any 
one member is 


U, = (1/2) [ail laa at (4c 
Since few of the members will be acted upon by more 
than one or two loads that do not have to be regarded 
as reactions, each a,;; matrix is made up largely of zero 
elements. The strain energy for the entire structure, 
in terms of the member loads, is the sum of expressions 


of the form of Eq. (4c) for all members. 


U= U+ U2+ + U, + (5a 
= (1/2)[giJlau]} is (5b 

where 
Lai] = lary] + Lavi] + + [aay] + (5c 
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The next requirement is to determine the strain energy for the structure as a function of the applied loads P,, 


and the redundants q,. 





This can be obtained by substituting Eq. (2c) in Eq. (5b). 


U = 1/2[Pmigr][vimiv ir] [orig] [Vin Vie] { Pride} (6a) 
where the primed symbol | |’ is used to denote a matrix in which the rows and columns are interchanged. 
Letting 
, Amni&ms 
Y ' 7 a Yi Vis = (6b) 
Gira | tire 
the energy becomes 
y l Amn! Ams { P ) 
= | Pmidr] Weph-eg Pris (6c) 
2 Cis Hite 
(3) Solution of Member Loads 
If Eq. (6c) is differentiated® with respect to one of the redundants q,, the result is 
ou l 1 . 
= [Pir] {amiars} + - [ary ais] } Pp qs} (7a) 
Od: 2 2 


Since r and s are interchangeable, m and n are interchangeable, and the square array of a’s in Eq. (6c) is symmetri 


cal about its diagonal, the two terms of Eq. (7a) are equivalent and the equation may be rewritten 


Ol’/0qg, = [amiars]{P,' qs} 


Applying the principle of minimum strain energy, 
Eq. (7b) must be equated to zero, and as many equa- 
tions of this form as there are redundants may be writ- 
ten by letting ¢ assume, successively, all the values of r. 


The resulting set of equations, 
(7c) 


may be solved for the redundant member loads to give 





fast — Qrs | orn {Pr} (Sa ] 
By letting 
Fa = |T,,) = — Ors 1] Orn | (Sb) 
this same result may be rewritten 
far} = [Drm] {Pm} (Sc) 


Eq. (Sc) can now be substituted back into Eq. (2b) 
to obtain the values of all the member loads in terms of 
the applied loads only. 


5 ( 1) ( Wr f ( . 
af = [yim] i Pm + [virl(T +m] } Pm} (9a) 


The same result may be written more compactly as 


717 
Las} om [Yim Yir| || {Pm} (9b) 
in which [J] = unit diagonal matrix. 
if 
[T'eeet = [Yemt¥ 1| =| (9c) 
L 1m i im ir. ‘ 
D ym 


then 


(7b) 


fait = [Tim]} Pm 


I';»| is the unit load distribution matrix for the 


(9d) 


and 
indeterminate structure. Each term I,,, is the value 
of the 7th member load when the structure is subjected 


to a unit value of the mth applied load. 


(4) Deflections of Redundant Structure 


If in Eq. (7b) the differentiation is performed with 
respect to an applied load P, instead of a redundant 
load g;, Eq. (7b) becomes 


OU/OP,, = [auniaus]{ Prigs} (10a) 


By Castigliano’s theorem, this is equal to the deflection 
6, at the wth applied load. Letting u assume, succes- 
sively, all the values for m, the deflections 6,, corre- 


sponding to the applied loads are 





{5m} - C= ams] Poids} (10b) 
Substituting Eq. (Sc) in Eq. (10b), 
I 
{5m} _ Cees | {P,.} (11a) 
Letting 
l 
[A -_~ - [amn Ams | r.. (11b) 
{om} = [A witli (11c) 


An element of the A,,,, matrix represents the displace- 
ment at the mth applied load due to a unit mth applied 
load. 
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(V) PROCEDURE 


Since all of the steps in Section (IV) are not required 
in solving a problem a summary is given of those items 
which are needed in obtaining the unit load distribu 
tions and deflections. In addition certain necessary 


details are discussed. 


(1) Summary of Formulas 


Eq. (2c), [Y¥im'vir| (load distribution from the equa 
tions of statics). 

Eq. (4c), [a@aij| (flexibilities of individual members 
to be obtained from formulas and curves). 

Eq. (5c), 

Eq. (6b), 


Qmn'Ams oe 
aah Se = LVemi Vir} Leg | LV int Vs 


Arn | Ars 


a;;| (sums of member flexibilities). 





(flexibility of determinate structure). 

Eq. (8b), [Pm ][T,,] = —[ars|~'La-,] (solution for re 
dundant loads in terms of unit applied loads). 

Eq. (9c), 


I 
[Pim] = [Yimivir] | = 
Evin 


(unit load distributions for indeterminate structure). 


Eq. (11b), 
ted = Fania | 
Amn| = |AmniAms r 


(flexibility influence coefficients of indeterminate struc 





ture). 


(2) Setting Up the Problem 


For the structure as a whole, certain of the external 
loads aresreactions; the rest are applied loads and are 
designated by P,,. If a reaction is also a redundant, it 
is designated by a g, symbol and treated accordingly. 
Any load that acts directly on a member may be desig 
nated by g, and may be included in the analysis; this 
includes applied loads, external reactions, and redun 
dant loads. However, it is of advantage to include only 
(1) those loads that will be of importance in interpreting 
the analysis, (2) those loads that are used to express the 
individual member flexibilities (7.e., at least one load per 
member), and (3) all redundants. 

In connection with the foregoing item (2), a distin 
guishing feature of this analysis procedure is the break 
down of the load distribution into individual member 
loads for each element. Although this may entail the 
inclusion of a large number of loads in the calculations, 
it greatly simplifies the determination of member in 
fluence coefficients, which may then be copied directly 
from the standard texts or from the collection of ele 
mentary cases in Section (VI). 

When the equilibrium equations are set up, they must 
be independent, and the total number of equations 
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must be equal to the number of member loads gi con 
sidered, less the number of redundants g,. Engineering 
judgment is still necessary in selecting the redundant 
member loads, although certain rules are available for 
determining both the degree of redundancy and the par 


ticular redundant loads. 


(3) Determining Member Influence Coefficients 

A collection of influence coefficients for structural 
Some are 
Section (VJ in 
cludes most of the structural elements that occur in in 


elements has been included in Section (VI). 
also available in standard references. 


determinate problems, but it is by no means complete 
and the flexibilities of other elements may be desired 
in some problems. Although the flexibilities of special 
elements may be calculated, it is well to remember that 
it is this aspect of a structure which is most unpredict 
able because of thickness tolerances, connection de 
flections, buckling, shear lag, etc.; consequently, the 
engineering approximation involved in using flexibili 
ties for standard members and loadings may often be 
justified. 

Since one or more of the loads on any member may be 
considered as reactions and since the choice of reactions 
is arbitrary, the energy of a member may usually be ex 
pressed in terms of one of a number of sets of influence 
coefficients. Care in selecting these reactions may allow 
the omission of certain of the qg;’s from the problem 
For this reason, certain alternative sets of influence coef 
ficients are given in Section (VI 

In some cases it may be an advantage to solve the 
problem by parts, obtaining the flexibilities of parts of 
the structure by the foregoing method and then com- 
bining the parts in a simpler final operation. 


Q «cross sectional area of bar 

Le length of bor 

E= Young's modulus 

q = end force or uniform shear flow 
Q, /Q, = ratio of end areas 
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Fic. 1. Flexibilities of bars. 
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ANALYSIS OF COMPLEX 


(4) Matrix Operations 
One of the large labor items is the solving of the equa 
tions of equilibrium |Eq. (1 statically 


determinate load distribution | Eq. (2c) 
have usually used the Gauss elimination method for 


to obtain the 
The authors 


this step, although other good methods are available. 

It is well to emphasize that all matrix operations 
should carry check rows or check columns in order that 
mechanical errors may be detected as soon after they 
are made as possible. This check should be made even 
if the calculations are performed by means of IBM 
equipment. 

Eq. (Sb) of the 
and its post multiplication by a, 


summary of formulas calls for the 


inversion of |a,, 


\lthough 
authors have usually found it more convenient to solve 


this indicates the end result desired, the 


Eq. (7e for [I,,,| directly by the Gauss elimination 
method. 

In a design problem, the applied loads P, 
ina number of different combinations corresponding to 
In this case the P’s may 


may occur 
different analysis conditions. 
better be arranged in a rectangular matrix, where the 
typical element P,,,, is the value of the mth applied load 


in the cth condition. Eq. (9d) is then writien 


| Jie Pen tit (9e 


FOR FLEXIBILITIES OI 
ELEMENTS 


VI) FoRMULAS AND CURVES 
CERTAIN STRUCTURAL 


The separation of a complex structure into its com 
ponents results in a number of elements whose flexi 
bilities may be obtained from simple formulas. This 
section presents the formulas for the simplest commonly 
used elements such as rectangular shear panels and 
constant section bars and beams. However, since 
these idealizations are not usually met with in practice, 
the more generalized cases such as bars and beams of 
varying section and shear panels of parallelogram, 
also included 


trapezoidal, and quadrilateral shape are 


(see Figs. 1, 2, 3, 4, and 5) 


(VII) AN EXAMPLE 

An example that demonstrates the method has the 
idealized structure shown in Figs. 6 and 7. As can be 
seen, it is a two-cell cantilever sweptback box, which is 
symmetrical about its chord plane, is tapered in three 
directions, and has a single elastic tip bulkhead. The 
idealization concentrates the direct load-carrying ma 
terial into bars at the front, main, and rear beams, plus 
a capstrip for the bulkhead, leaving the covers and webs 
loaded by edge shears only. This particular structure is 


selected only for convenience in demonstrating the 


method; in an actual wing the idealization would have 
to be more detailed. 

The loads have been broken down into four applied 
member loads, of which three are re 


Pris. . 21 


loads and 17 


dundant. The loads consist of a 


applied 
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= flexural rigidity of beam 


BEAMS 


E 
| g = length of beam 
q = end moment or end shear force 
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Fic. 2. Flexibilities of beams 
A=surface area of panel 
SHEAR t=thickness of panel 
PAN FLS G= shear modulus 
q= shear flowin Ibs per inch 
Case. se : 3 s A 
3.1 7 q, | T Gt 
___ | 
Case | 9 ove) —_, - | - : A 
3.2 nat | he | a: Gt or ore * oy 
33 ‘ ; /, A | al,,,.8 Gt (due to shear stress only) 


|Shear acting on an oblique edge also 
M4 / 
So, | “2qjl conf gives rise to an axial force which may 
jbe split between the adjacent bars and 
Section showing actual interno! | tweeted as an additional end load for 
leach bar. See Case 12 for bor coeffi- 








loads 
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aie 3 cients. This procedure allows for 
4 q, | the interaction between bending stress- 
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Applied loads and structural element dimensions for 
example structure. 
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ANALYSIS OF COMPLEX 


bending moment and shears at the front, main, and 
rear beams, respectively, all acting at the outboard 
end. The bending moment Pi; is assumed to enter the 
structure distributed according to engineering beam 
theory and is in units of 10 in.lbs. The shear forces 
Py... are in units of | lb. 

The breakdown into structural elements, together 
with member load designations q). . .17 is illustrated in 
Because of symmetry, only the member loads 


In the case 


Fig. 8. 
above the chord plare need be considered. 
of each shear panel, the average shear flow on only one 
of the edges is represented by a qg because the average 
shear flow on each of the other three edges may be 
given by this g times the proper numerical constant, as 
determined by statics. This permits the omission of 
three nonessential internal loads for each panel. 

The member loads 415, gis, and giz shown at the in 
board ends of the three bars represent the axial stress 
components of the shear in the swept covers. Since 
the idealization of the structure places all axially effec 
tive material in the bars, this component must be in 
cluded with the bar loads in the energy calculations 
even though it is not physically applied to the bar ends. 
Hence, the loads are shown dashed. Values may be ob 
tained for these loads by use of shear panel Case 3.4. 

The three member loads selected as redundant in this 
inalysis are the shear flow in the front beam web g; and 
the forces at the inboard ends of the front and rear 
bars gu and qi3. The 14 equilibrium equations, used to 
find the determinate load distribution, may be described 
as follows: 

Eqs. (1), (2), and (3), loads at outboard ends of bars 
due to tip bending moment P;. 

Eqs. (4), (5), and (6), spanwise equilibrium of bars. 

Eq. (7), bending equilibrium at root section. 

Eq. (8), moment equilibrium for tip bulkhead in its 
own plane. 

Eq. (9), vertical equilibrium at attachment of rear 
beam to tip bulkhead. 

Eq. (10), chordwise equilibrium of forward cap of tip 
bulkhead. 

Eq. (11), chordwise equilibrium of rear cap otf tip 
bulkhead. 

Eqs. (12), (13), and (14), determination of axial com- 
ponent of skin shear load. 

The dimensions entering these equations are scaled 
from a '/s-size plan view and therefore represent 
lengths projected onto the chordal plane. The fourth 
equation is typical; it places the front beam capstrip 
in equilibrium. 


qu + 102.3(0.7719¢: — 0.661gs) — gs = 0 


Here, and also at the stage where the member influence 
coefficients are written down, mere slide rule accuracy 
The machine calculations representing 
the subsequent steps require greater precision in order 
that the final results may be as accurate as the initial 
data. 


REDUNDANT 


STRUCTURES 


683 









































) xara 4 ? 2 [ala {alot 7) 
Fic. 11 aij}, sum of member flexibilities 
due to unit load Pm 
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Fic. 12 'I'ym], load distribution in redundant structure. 


The next step is to obtain the flexibility influence 
coeflicients for the individual members of the structure. 
These are written using the formulas and curves of 


Section (VI). 


Case 1.2 is used for the bars even though 


the edge shears in the tapered panels are actually vari 


able rather than constant. 


The difference should be 


negligible here, although in extreme cases it would be 
necessary to take this factor into account 
A typical flexibility calculation is that for the front 


beam capstrip. 


Applying Case 1.2 of Section (VI), 


102.3 X& 0.804 


This bar is shown in Figs. 8 and 9. 


Aa = Au, = A515 = 18.30 
; ; 3x15 
102.3 X 0.860 A 
= ; & 
Mall, Qalb5,8 5 = ‘ 
— S xX i. 
Qqs,s = (102.3 X& 0.928)/(38 XK 1.5) = 21.1 


For convenience, the modulus of elasticity is taken as 


EK = 


1 and Poisson's ratio is 0.3. 


The influence coef- 


ficients for the other members are calculated in a like 
manner, and all are assembled in the a,, matrix, Fig. 11. 
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due to unit load Pp 
































© NUE 4Q_| 20 | 2/ 

© E/18B| 38O!}W/7 | V8.6 125.1 
5 %1/9| 7.6 | 9/0/ | 6488 | 7997 
a &\20|78.6 | 6488 | 698.0 | 676.7 
= 2/| 125.1 | ¥99.1 | 6764 | 926.2 
Fic. 13. [Am,»], flexibility of redundant structure. 
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Fic. 14. Axial stress distribution at root compared with engi- 
neering beam theory. 


The remaining work is entirely mechanical. The 
solution of the equilibrium equations yields the [y¥;,,' 
vir] matrix of Eq. (2c), and this result is shown in Fig. 
10. Note that one’s are inserted where a member load 
is also a redundant—that is, for g1, gu, and giz. Follow- 
ing this, the matrix multiplications listed in the sum 


OCTOBER. 1952 


mary of formulas, Section (VI), may be carried oy 
The final result in the form of the unit load distributioy 
[I';,] and the flexibility influence coefficients [A.,,,] fo, 
the redundant structure are given in Figs. 12 and 13. 
It is interesting to see that even this simplified struc 
ture exhibits the characteristic behavior of a sweptback 
wing with built-in root. Perhaps the most significant 
feature is the concentration of axial loading in the rear 
beam at the inboard end. The variation from eng; 
neering beam theory distribution in this region is shown 


in Fig. 14. 
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A Second-Order Solution for an Oscillating 
Two-Dimensional Supersonic Airfoil’ 


ALEXANDER WYLLY? 
The RAND Corporation 


SUMMARY 


[his paper presents a second-order solution for the forces and 
moments produced by an oscillating two-dimensional airfoil of 
urbitrary cross section. This solution, valid to the second order 
in angle of attack and thickness, is obtained by means of an iter 
ition procedure. 

One part of this paper is concerned with the development of a 
linearized (or first-order) velocity potential, which is of suitable 
form for use in the iteration procedure. This velocity potential is 
developed as a power series approximation, in frequency, to the 
exact linearized solution. The first six terms of this series are 
shown to be within a few per cent of the complete linearized solu 
tion when the reduced frequency (@ = wC/ Bx) is less than 
1.3. 

The second part of this paper is devoted to the development of 
the second-order solution in angle of attack. This latter is de 
veloped with the use of the above-mentioned series approximation 
to the linearized solution, and, since the linearized solution valid 
to the second order in frequency was utilized in the iteration pro 
cedure, the final solution is valid to the second order in angle of 
attack and frequency 

The second-order terms were found to have a relatively strong 
influence on the final solution, particularly for the pitching mo 
ment. It will be seen that in many cases the second-order terms 
are larger in magnitude than the corresponding first-order terms 
and, thus, reverse the tendencies indicated by first-order theory 
In particular, it was shown that the theoretical instability pre 
dicted by linearized theory for an airfoil of zero thickness is 
completely eliminated for an airfoil having a thickness ratio as 
small as 4'/s per cent 

The damping moment for a 7'/2 per cent thick airfoil was com 
pared with that given by British experimental tests. The second 
order solution is shown to be in much better agreement with the 
experimental data than the linearized theory; however, the 
magnitudes of the damping moments predicted by second-order 
theory are appreciably lower than those indicated by experi 


ment. 
SYMBOLS 

( = speed of sound 

t = chord of airfoil 

CI = lift coefficient 

Cy = pitching moment coefficient 

Cp = pressure coefficient; subscripts u and / refer to upper 
and lower surfaces 

f = an arbitrary function of (x — y8), which is defined as 
f(t) = f(t) = O fore < O; f’(e) = Of(¢)/d¢%, and 
fand f’ are continuous 

g = an arbitrary function of (x — y8), which is defined as 
e(t) = g(¢) = Ofor¢ < C, and g and g’ are con- 


tinuous 


Received July 30, 1951. 

* Based on Doctoral Thesis,! California Institute of Tech- 
nology, 1951. 

t Research Engineer, Missiles Division. 


H = a function of the first-order (in @) velocity potential, 
which forms the nonhomogeneous portion of the 
wave equation to be satisfied by the second-order 


velocity potential 


p = pressure 

q = perturbation velocity = V u? + v? 

t = time 

i = thickness of airfoil 

u,v = perturbation velocities in the x, y directions 

Vi = free-stream velocity 

x, ¥ = coordinates, fixed in the leading edge of the airfoil with 
the positive x axis in the direction of the free stream 

Xe = nondimensional position of the point of rotation 

a = angle of attack = ay cos wl 

ay = maximum angle of attack, ay < <1 

8 = VM,? -1 

p = density 

? = complete velocity potential 

’ = complete perturbation velocity potential 

#,” = perturbation velocity potential correct to mth order in 
angle of attack a and nth order in frequency w 

gn” «= one term in series expansion of #,”, it includes the 
mth term in a and the nth term in frequency 

¢ = the portion of ¢,? which satisfies the nonhomogeneous 
wave equation; the “particular solution’”’ 

a = the portion of ¢,? which satisfies the boundary condi- 
tions on the surface; the ‘‘correction potential” 

y,@ = terms used in place of g and o for the discussion of 
more complex airfoils 

w = angular frequency in rad. per sec. 

@ = the reduced frequency = wl/ BV. 


( )o = (subscript zero) free-stream condition 


(I) INTRODUCTION 


— LINEARIZED THEORY for the supersonic, irrota 
tional, nonviscous flow of a perfect gas past an 
oscillating wing is at present in a reasonably complete 
state of development. Progressing from the first in 
vestigations of nonstationary flow, by Possio* and von 
Borbély,* many authors have considered the two 
dimensional oscillating airfoil problem. More recently, 
a number of linearized solutions*~* have been deter 
mined for nonstationary three-dimensional wings. 
Thus it seems probable that within a few years the 
first-order solution for almost any plan form will be 
available. The purpose of this paper is to extend the 
linearized solution for an oscillating two-dimensional 
airfoil to include the higher-order terms in angle of 
attack and thickness so that the relative importance of, 
and the tendencies indicated by, the second-order terms 


may be evaluated. 
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This paper is divided into three principal parts. The 
first part gives a general discussion of the problem. 
In this section the iteration procedure is applied to the 
complete nonlinear potential equation, and the partial 
differential equations that must be satisfied by the first- 
and second-order velocity potentials are determined. 
The first- and second-order boundary conditions are 
also discussed in this first part. The differential equa- 
tions and the boundary conditions then completely 
specify the problem. 

In the second part, a first-order solution, whose form 
and properties are satisfactory for use in the iteration 
procedure, is developed. In the iteration process each 
step of the iteration is dependent on the previous solu 
tion. Hence, to determine the second-order velocity 
potential, a first-order potential that is defined through- 
out the entire x, y plane and which is of a compara 
tively simple closed form must be available. A linear 
ized velocity potential that satisfied these requirements 
was determined by Guderley’® in 1946. In this paper 
Guderley produced the first two terms of the series ex 
pansion (in frequency?) of the linearized velocity po- 
tential. 
the first two terms and for the first six terms are pre 
sented. These throughout the 
whole x, y plane and are in a simple closed form. 

In the third section of this paper, the second-order 


In the present paper, the expressions for both 


solutions are valid 


solution in angle of attack is developed for the case of 
a flat plate, a wedge airfoil, a double-wedge airfoil, 
and an airfoil of general cross-sectional shape. The 
lift and pitching moment coefficients for several modi- 
fied double-wedge airfoils with different thickness ratios 
are compared with the coefficients as calculated by 
linearized theory. damping moments 
predicted by first- and second-order theory are com- 
pared with those given by British experimental data. 


Finally, the 


} There is a tendency for the expression ‘‘order’’ to become con- 
fusing in this paper, since order can refer to the maximum angle 
of attack ay, the thickness parameter e¢, and the frequency w. 
The maximum angle of attack and the thickness determine the 
order of the solution in the usual sense of the word—i.e., the 
first order (in ay or e) is the linearized solution. The maximum 
angle of attack and the thickness parameter are taken to be of 
similar magnitude, ay is of 0(e), and for brevity are referred to as 
order in angle of attack. Any order solution in angle of attack 
may then be expanded into a series solution in powers of the fre 
quency. Terms of this expansion are referred to as order in fre 
quency. 
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(II) GENERAL DISCUSSION OF PROBLEM 


(A) The Iteration Procedure 


The problem to be considered is the study of the syr 
face pressures produced by the uniform flow of a per 
fect gas past an oscillating airfoil of general shape (see 
Fig. 1). The airfoil is required to be slender and the 
maximum angle of attack ay is required to be small 
The airfoil rotates about a point (x = xoC located a 
nondimensional distance x) downstream from the lead 
ing edge. The angle of attack is defined by the symbol 
a = ay cos wt, where w is the angular velocity of the 
oscillation. The angular velocity is limited by requir 
ing that the reduced frequency be small. The flow js 
uniform upstream from the leading edge and is irrota 
tional, isentropic, and nonviscous. The assumption 
that the flow is irrotational and isentropic is valid in 
this problem, since the effects of rotationality and 
changes in entropy are of third and higher order in angle 
of attack." 

Under these assumptions we can derive the nonsta 
tionary Bernoulli equation [Eq. (1)| and the two 
dimensional nonstationary potential equation [Eq 
(2)] from Euler’s equation, the equation of continuity, 
the condition of irrotationality, and the isentropic ex 
pression for the speed of sound. 


= ®,”) = Sc = ®,”) oe ®;, oa 


2, P JP, y= 20.P,, — 28,?,, = () yy) 


®,,(c? 


The complete potential equation, Eq. (2), is nonlinear 
and cannot be solved directly. The approach to be 
used is that first suggested by Prandtl. It is assumed 
that the velocity potential @ may be approximated by 
the expression 


&® = Vox + ag! + a’g? + a®gi+...ta¢ 3 


The term Vox has been included in Eq. (3); hence, the 
n 


¢’s represent the perturbation potential, ¢ 
Since in this study only 


represent 
ing the nth-order potential. 
the second-order potential is desired, the terms ¢’” for 
n > 2 may be neglected. Putting the series expansiot 
for & [Eq. (3)] into Eq. (2) and replacing the speed o! 
sound c in Eq. (2) by its corresponding function of 4 
from Eq. (1), the potential equation becomes: term 


containing a, 


here o/c) has been replaced by the Mach Number 


/ ce 9 
Myand B = V M,? — 1; terms containing a’, 
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AN OSCILLATING TWO 


(Grr! ss Puy’) xX 
l l 
(e.! + V ¢:) T 2¢,! (ea + V, en') + 


I . 
2Z¢,' (en + V, Put! (9) 
0 


These are the differential equations that must be satis- 
fied by the first- and second-order potentials ¢' and 
e. Both of these equations are linear and may be 
solved by the usual methods applicable to linear dif- 
ferential equations. The method of solution is now 
apparent. First, one obtains the first-order solution 
y! that satisfies Eq. (4), the homogeneous wave equa- 
tion, and the proper first-order boundary conditions. 
Then, the function //(¢') iscalculated. Next, a second- 
order potential ¢*, which satisfies Eq. (5) (the non- 
homogeneous wave equation) and the proper second- 
order boundary conditions, is determined. 

The expression for the speed of sound c will be 
needed for the determination of the pressure coefficient. 
From Eqs. (1) and (3), 


c= oc — (y — 1) {Voh,? + B,? + 
(1/2) [Vo? + (@,')? + (@,')7]} (6) 


The pressure coefficient C, is determined by substitut- 
ing the ratio (c/co)? from Eq. (6) into the isentropic 
relationship between the pressure and the speed of 
sound. The expression for the pressure coefficient is 
found to be 


: Eee I ‘ ‘ 
bial r Vo? i 7 2 V0? | 0. - 


f 


( 0 


, | 2Me? a 
(4y1)? +" lbh + To (bE ]E (7) 


The third- and higher-order terms have, of course, been 


neglected. 


(B) Boundary Conditions 


The flow is required to satisfy two boundary condi- 
tions: 

1) The resultant velocity must be tangent to the 
surface of the body. 

2) All perturbation velocities must vanish up- 
stream from the plane x = 0. 

From the theory of hyperbolic differential equations 
itis known that these two requirements are sufficient 
to determine the solution. 

First, consider the tangency condition. If the sur- 
face of the plate is given as y = h(x, ¢), the tangency 
condition at the surface of the body is 


{[@,2 — (Oh/dt)]/(Vo + ,*)},-, = Oh/Ox (8) 


Here, the numerator on the left represents the differ- 
ence in vertical velocity between the flow and the plate, 
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while the term (Vo + ®,") is the resultant velocity in 
the direction of flow. If ®* is replaced by ag' + a’¢’, 
where ag!', 0h /Ox, and Oh /Ot are of the order a while 
a’g* is of order a’, the boundary condition at the sur 
face to first-order in ais 


(¢gy')y=90 = Vo (Oh/Ox) + (Oh/OF) (9) 


The second-order tangency condition may be deter 
mined by setting @? = ag! + a’yg’ and subtracting 
Eq. (9)]. 


out the first-order equation The second 


order tangency condition is found to be 


(¢y7)y=0 = (Oh/Ox) (¢2')yao + [ley )y=0 — (Gy) yaa] 
(10) 

Now consider the condition that the perturbation 
potential vanish for x < 0. This requirement may be 
stated and satisfied in two physically different, and 
yet mathematically similar, ways. Since both defini- 
tions are used later, both will be presented. 

The first approach requires that two conditions be 
imposed on the so-called timelike surface x = 0 and one 
The 
previously mentioned tangency condition is spacelike; 


condition along the spacelike surface of the body. 


hence, we have to impose two conditions on the plane 


x =0. These are 


(0, y, t) = 01 


(11) 
©,(0, y, t) = OF 


These conditions, along with the differential equation 
and the tangency condition, completely specify the 
problem. 

The second approach specifies the potential upstream 
by requiring that the potential be zero on the Mach line 
from the leading edge of the airfoil. 


¢g'(x, x/B, t) = 0 (12) 


As expressed by Eq. (12), the potential goes to zero on 
the Mach wave through the origin and not on the oscil- 
lating Mach line from the leading edge. However, 
since the movement of the airfoil is considered to be 
small, the movement of the Mach wave only produces 
effects of a second order in angle of attack. Therefore, 
since this criterion will only be used in the determination 
of the first-order solution, the second-order effects may 
be neglected. 


(III) Frrst-OrRDER SOLUTION IN ANGLE OF ATTACK 


In 1946, Guderley'® wrote a paper entitled ‘The 
Pressure Distribution on a Flat Plate Oscillating with 
a Small Frequency.’ In this unpublished paper he 
produced, without any development, the first two terms 
of the series expansion in frequency (of the linearized 
velocity potential) for the flow over an oscillating flat 
plate. A similar solution will now be developed. The 
method of attack is interesting because it allows an alge- 
braic determination of any desired number of the terms 


of the series expansion. The method is not restricted 
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to oscillatory motions but actually appears to be useful 
in the solution of any two-dimensional nonstationary 
supersonic airfoil problem. 

The location of the surface of the airfoil (see Fig. 1) is 
given by the relation y = a(x — xC) (for0 <x < C), 
where a = ay cos wf. The tangency condition, Eq. 
(9), is, thus, 


(gy')y=o = Voa + a(x — xoC) (15 


Now consider the series expansion of g! in powers of a 
nondimensionalized frequency @. The so-called re 
B°1y) is initially required 


duced frequency ® (@ = wC 
to be small (a < 1) in order that the higher order 
terms in frequency may be neglected. Later it will be 
seen that the sixth-order solution gives a good approxi 
mation to the exact solution for a < 1.4. The series 


approxinfation to ¢! is 


g! = gy! + ag! + wy3' +...+ 
®,4,;' (14) 


where ¢,,! is the nth-order term in the series. It should 
be noted that, while the so-called first-order term in 
frequency does not contain @, this term does contain 
the parameter a = ay cos wl. Thus, this term is due 
to the instantaneous position of the airfoil and is the 
first-order term in the frequency expansion. 

From Eq. (13) it was seen that ¢,' at y = 0 was of 
the form 


. 


(gy )y=0 = Voa + a(x — xoC) 


Hence, ®,' must be 


i 
#,'= f- (¢,'),-0dy + afAx + Dy? + Hxy + 
0 
Fx? + Jx? +...) + a{ Bx + Ey? + Gx? + 
Ixy + Kx? +... (15) 
where A, B,..., K are functions of x, C, Mo, and ve 
locity. The higher terms in frequency (i.e., &, a, etc.) 
are not included, since only a second-order solution in 
frequency is desired. ,' may now be determined 
algebraically by substituting Eq. (15) into Eqs. (4) 
and (12). When this is done, it will be seen that a set 
of n algebraic equations is derived which is exactly 
sufficient to determine the ” unknown coefficients A, 
B,...,n. To simplify the algebra, the unknown 
coefficients /, ..., will be taken to be zero immedi- 
ately, although the same result would be obtained if 
they were allowed to remain through the derivation. 
By integrating Eq. (15), it is found that ®,' is of the 


form 


@,' = yVoa + yxra — yaxoC + afAx + Dy? + 
Fx? + Hxy} + a{ Bx + Ey? + Gx? + Ixy} (16) 
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Now, substituting Eq. (16) into Eqs. (4) and (12), i, 
follows that 


O®.! = 0 = B2(2Fa + 2Ga) — (2Da + 2k a) -}. 
(2Mo?/Vo) (A& + 2xFa + 2Hya) + O(a 


7 ae. Voax xa xx Ca 
@,' |x, --,f] =0= + -+- ‘J 
re) B B B 


Dax? 
Aax + Bax + — oa 
B- GB 


Fax? ——— ; 
+ Fax? + Gax* a 
Tlax , lax 
ne 
ie} A 
Rearranging these two equations according to the vari 
ables a, a, x, and y, a set of eight simultaneous alge 
braic equations with eight unknowns is determined 
Again note that, if twenty unknown coefficients had 
been used, there would have been twenty simultaneous 
equations and the last twelve coefficients would be 
given as zero. Solving these eight simultaneous equa 
tions, the unknown coefficients are found to be 


B=nxC/6B: D=F=H=I1=0 


Hence, to the second order in frequency the linearized 


velocity potential is given as 


@,! = y*(a/28) (1 — 2My?) + vla(x — mC) 4 
aVo] + (x?a/2B8*) + (x/B) (xoCa& — &Vo) (17 


Eq. (17) is the solution given by Guderley for the po 
tential of a flat plate oscillating such that o < | 
As will be seen later, this solution approximates the 
exact linearized solution within a few per cent for @ < 
0.3 (this would include such motion as control move 
ments, stability, etc.). 

The previously determined second-order solution will 
now be extended to include the linearized solution valid 
to the sixth order in frequency. Since the method 1s 
exactly the same as that given previously for the 
frequency, only the final 

The determination of the 


second-order solution in 
equations will be presented. 
solution is further simplified by the two following 


rules: 


(a) No term Ax"y may appear, since ®, + 0 and 
the tangency y conditions would be changed. 


(b) The symmetry shown in the first two terms ol 
the series expansion is continued throughout the higher 
terms, and thus ¢,! = F(x”"y’), where [m + p = n,0S 
(mand p) <n]. 

The final sixth-order solution in frequency which 
satisfies the requirements of Eqs. (4), 9, and 12 tol 
(@*) is presented below. 











In t 
quit 


h,! 


The 
first 


in E 


gIVe 


The 


r 
pare 
ized 
the 

com 
(M, 

Ba. 
ten 


lax 


Vari 
alge 
lined 
> had 
1eous 
d be 
‘qua 


AN OSCILLATING TWO-DIMENSIONAL 


x '  /1-—2M,? {1 ~ x 
,' = Voa (» + aly ( Tx TX ( — Kel (y _ ) _ 
: B ' 28 2 B° B 


AIRFOLL 


I= Ga) -*G)+ 


689 


aM," 
Vo 


28? — | xo (x? R aM)? J 1 + 4M,? 2 — 12M,? + 8M! xy? xy? 
vy ( : ai Sa + - ye" saa FF ; a 
e $B” 28 B- Vo" 48 8" 1S 8° 6 SB 
ms se Moe + 1) + xy? 4 = aMy* { . (° + | Cy LOx*y*®  xy* U4 
") = (2 0 T ) P = a ale a — OF : : r : ( S sig 
as 128° 18° ' 6]f 192 V3 ' 8° a” B 
= =<" 6x-y" y! 1M! SA! 9M)? 1 
20 My? + 15) — 4x06 | - (2M + 3) — + (3 - 2,2) | —— | ( = r ) ' 
rey B® 38 f 1920 Vy! gt 
Saty? | Lx xyC , ; 10x9C. xy , L5x*y* 
— (15M)? + 1) + — (28Mot + 109M)? + 8) + ae (11M)? — 6) + 7 (1 — 3M)?) +4 
B° (pr (| : p' 
is Ye. ; 20x 9Cxy' nate 
l6xy> — 16xyCy® + — (15 — 99My? + 121My* — 48M.) + —S—— (25M? + 8) ] (18) 


In the calculation of the linearized pressure on the oscillating plate, only the potential on the plane y = 0 is re 


quired. This potential is presented in Eq. (19). 

| Voax ca ( x? “) aM? ( x? uc) ; aMeT . (' - Ms’) 
(d,*) ) + — : - 4 — x? - 4 
bs - B . 28 B | 0 + 6° 26° | 0” 1S 8° 


XoCx? 
| 2 B° 


aM,‘ 3 1M,? Lay Cx! 
seed Je? us ) - aes 3) | i 
192 Vo? 3? 3 | 


_(8Mot + 9M? + 1 fay Cx? 
x” gil + 


(2M,? + | - 


aM! 


(28.494 + 109A7,? + 8) (19 
1,920 Vo! . | 








ized 


nd 








76° 
[he fourth-order solution included in Eq. (19) has previously been obtained by Watkins,'' who expanded the 
first-order solution in powers of @ and then integrated to obtain the fourth-order solution in frequency. 
The linearized pressure coefficient (C,) may be determined by substituting the potential &,' [Eq. (18)] for ® 
in Eq. (6) after the higher-order terms in a have been dropped from Eq. (6). Then, the section lift coefficient is 
given by the expression 


a + 
te ( inet 1) ~ le + + Ga + [4M,° + 27M. + 4M? + 
12 960 
, a @ , j J , a ‘. i 
20x09 B?(41194 + Mo I ne wt ) (Lo 2 2x 987) + ry (Mot + 4Mo? + 12x08?M)?) —- 
a? | 
| zm + 16My® + 9Mo*) + 30x» 82?(19Mi® 4+ 10. | ) (20) 
20,160 f 


lhe pitching moment coefficient as calculated about the point of rotation Xp 1s 


| r | a 4 
Cy cos wf ) xo} — (My? + 2 + xX (4M? — 8) — 8xo?B?] + = | (20M? 4+ 
/ ) ) 16 9,670 


135Myt + 20My2) + x0(360Mo® —396Mo! 120.7) 120x97B?(4.Mo* + un) | 2 sin wt (200," 
) 


wo’? M, ; 
1 — 3x (28° 1) + 6x9?6?] + 0 [((4Mo? + 16) + x (4011.2? — 65) — 60 x9?8?] — 


otM 4 


6,720 


(24M, + 192M)? + 108 + x0(936.Mo* — 552Mo? — 563) 30x07? B2(19M)? 4 i]t) (21) 


each other. 


The solution for the lift, Eq. (20), will now be com 
pared with the lift as determined by the exact linear 


ized solution. The exact solution was computed from 


the tabular data presented by von Borbély.* The 
comparisons will be made at three Mach Numbers 
(Mp 10/9, 10/7, and 10/5) and plotted versus 
8*@. The reason for plotting versus 6°@ is that it 


tends to separate the results that otherwise lie close to 


As may be seen from Eq. (20), the lift is 
given in the form 


Cr = (4ay,/B) (P cos wt + Q sin oat) 


where P and Q are functions of (@, Mo, and xo). Fig. 2 
shows the function P as calculated from the second- 
order, fourth-order, sixth-order, and exact solutions. 


Fig. 3 shows the function Q under the same conditions 
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Both Figs. 2 and 3 are presented for the condition x, ~ 
Q—that is, the case when the airfoil is rotating aboy 
its leading edge, It will be noted that for relatively 
low-speed motions such as are found in stability calcula. 
tions or control deflections the second-order solution jy 
frequency is within a few per cent of the exact solution, 
For high-speed motions such as flutter, the sixth-order 
solution may give a sufficiently accurate result. 


(IV) THE SECOND-ORDER SOLUTION IN ANGLE of 
ATTACK 
(A) The Flat-Plate Solution 


When the first-order potential [Eq. (17)] is substi 
tuted into the right-hand side [//(¢')]| of the nonhomo- 
geneous wave equations [Eq. (5)], another difficulty 





becomes apparent. Consider the first term (¢,!) of 
Eq. (17). 


Voa; the second derivative is thus zero. 


The first derivative with respect to y gives 
The second 
derivative is then zero ahead of the line x = y, since 
¢i! is defined as zero ahead of the Mach line from the 
leading edge and is also seen to be zero downstream from 
this Machline. However, since.(¢,'), is a discontinuous 
function, (¢'),, is infinite along the Mach line from the 
leading edge. The second derivative is then given as 
zero everywhere except at the Mach line from the lead- 
ing edge at which point it is infinite. Since it would 
be difficult to work with 
another method of attack which eliminates the dis 


this velocity potential, 
continuous derivatives was chosen. The method was 
that employed by Van Dyke,'” and is also mentioned on 
page 365 of Courant and Friedrichs,'* for the station 
ary flow problem. 

In attacking the steady-flow problem, Van Dyke gen- 
eralized the airfoil shape so that the ordinate of the air- 
foil is given as y = ef(x). The function f(x) is unde 
fined, except that the surface and the slope of the sur- 
face are required to be continuous and to go to zero at 
Now 


there is no discontinuity in velocity occurring at the 


the leading edge (i.e., f(x) = /’(x) = Oatx = 0). 


leading edge or, for that matter, anywhere in the flow 
field. 
expression for the pressure was determined. 


The iteration process was then utilized and an 
This solu 
tion showed that the second-order pressure at any point 
on the airfoil depended only on the slope ¢f’(x) of the 
airfoil at that point. Therefore, suppose that ¢(x 
was required to represent a flat airfoil everywhere ex 
cept in the region 0 < x < 6 where the original restric- 
tions on f(x) are still enforced. The general solution ts 
still valid for this problem. Now the solution remains 
valid no matter how small 6 is allowed to become, and 
as 6 — 0 the pressure is known everywhere on the air- 
foil except for an infinitesimal distance aft of the lead 
In the limit where 6 = 0, the pressure is 
This is 


ing edge. 
known everywhere except at the point x = 0. 
the flat-plate solution. Van Dyke’? has further shown 
that, in the limit, the velocity components and the pres- 
sure are known everywhere in the flow field except in 
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an-shaped region y = [(x/8) + v] + vextending f = f(x — ys) = x — ys, must be determined. A 


the f 
fom the leading edge of the airfoil. The expression general potential that satisfies these requirements is 


was shown to be equal to the angle that the shock found to be 


for V 

s given by second-order theory) lies ahead of the | Ve? Cy 

: a Ry : Mo*x} Xo 

Mach wave, OF for an expansion 2v represents the ,'= —al i( ) + «| ihe: eee ae 
5° re} 


width of the Prandtl-Meyer expansion fan. The iter 


; is cet a ae eae re. ee ee M,? a 
ation process, as applied to the case of an airfoil of zero ( u ‘ ') f(¢) dt (22) 
0 


thickness by means of the “general solution” attack, B* 8 
was thus shown to be successful except within the pre 


: ee P . ares _ Eq. (22) reduces to Eq. (17) when f(x — y 8) is set equal 
viously mentioned fanlike region of order ¢ lying near : omen vi ’ 
: bn ap : to x — y@ and obviously satisfies the boundary condi- 
the Mach line from the leading edge. ; = ' “ ‘ 
uid tions [Eq. (11)] that @ and ®, are equal to zero at x = 
Asimilar procedure will now be applied to the problem : Pare ag ; 
0. Also, Eq. (22) satisfies the requirement ,' = 


an oscillating flat plate. A first-ordert velocity 


0 


‘ ‘ 0 + O(@), since all of the first-order terms cancel 
potential @,'(f, f), such that #,'(/, f’) becomes the pre 


j : : : , identically. Thus, #' [Eq. (22)] satisfies the wave 
viously determined first-order solution ®,'(x, y) when St ae : <a tee 

‘ equation [Eq. (4)] and the boundary condition | Eq. 
+ For the remainder of this report, order will refer to angle of (11)], and for the case f(f) = ¢ the boundary condi 


attack unless otherwise stated tion is that of an oscillating plate [Eq. (13)]. 


Now, having determined a satisfactory first-order general potential #,', a general second-order potential g* must 
be found which will satisfy Eqs. (5), (10), and (11). Substituting Eq. (22) into Eq. (5), it follows that 


” 


a Vily + 1) Mo 2MX(y+1) 2 . 2Mol(y + Imo 
O¢? = . rf" + ad xf'f” + . f's 
B? B 8° 
2Mo*(yv + 1) M,)*(3 — 2Mo?) (vy + 1),. 
- ff” + Ye prys (23) 
B- B 
Now let g? = ¢ + a, where ¢ is a particular solution that satisfies the nonhomogeneous wave equation, Eq. (23), 


and o is the correction potential that makes ¢? satisfy the boundary conditions [Eqs. (10) and (11)]. A function 


¢ that satisfies Eq. (25) is 


My *( 1) M,? ics M,y*( L)ayC 
e= ats| — 1% w(f’)? + Sie [ (f’)? dé + i A uy] + 
0 


12° . +34 4° 
M*(y + 1) an Mo‘(y¥+1) ,, 
ad = MC — xff (24) 
26° B4 
That y does satisfy Eq. (23) may be verified by direct substitution. It is also apparent that ¢ = g, = Oatx = 0- 


Now o must satisfy the wave equation, the boundary conditions [Eq. (11)], and the tangency condition 


(Fy), a —(Py)5 0 7 (af’®s,'), 0 + [(Psy"), a (Poy") sur (25) 
Let the surface of the general airfoil be given as y = af’(x — xpC). Now consider Eq. (25); the first term may 
be determined from Eq. (24) and the second term from Eq. (22). To determine %»,' at the surface [see the third 


term of Eq. (25)], expand f and f’ ina Maclaurin series—.e., 
af(x — YB)yaariz—mc) = of(x) — a? B(x — xoC) [f’(x)]? + 0(a*) (26) 


Thus the equation for (0,),—~9 is complete, and it only remains to determine o. The expression for @ is given be 
low [Eq. (27)]. As given, the equation is seen to be lengthy; however, the starred terms are only included so that 
the expression for o will be complete, since these terms go to zero in the case of a flat plate [i.e., when f(¢) = ¢ 
and, hence, f’(¢) = 0]. For purposes of brevity, the method of determining ¢ is not shown (for further details, 


see reference 1). 

















JOURNAL OF THE AERONAUTICAL SCIENCES—OCTOBER, 1952 


ee oe Tye ae L(y + (1) (Mot + Mo?) — 4M0?28? is 
ao = a’Vo wy aie | {() dé 438 + xCf’ - 
mn [7 oy a | Moly + 1)7*) f  (2Me ‘90 
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: 1) (M° M,') — 48°M¢? > - {2M ce ili M4 1) ]* 
1 = + P | ~ rca” ( = ) + wt ‘i id az | Y7 | \ 
2B B- 0 B* 
i 7 2M,?\* ym [ 2M,?\* me 2M? 
ae aa ’ " - — ww le f a oe n 
(x yp) 4 } ar ( =] | , jen) dn a ( 1 + y| f as ( 3? 
2Mo' Mo’ + My*]5 I aS ia ee 
! 38 cl oo dis ae 236 ii (x — yB) [ (f’)? d& — 4 Jy Lf} (n) |? dn de 
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The complete solution ®,” valid to the second order in frequency and angle of attack is 
$,/ = @' + eto 28 


where #,', g and a are given by Eggs. (22), (24), and (27). 


The second-order solution for the forces and moments on a flat plate can be determined from Eq. (28); however, 


since the second-order terms are zero for an airfoil of zero thickness, the solution is the same as that given in Eqs 
20) and (21). 


B) The Wedge Airfoil 


Consider an oscillating wedge that has a thickness ratio 2e [thickness = (2e) chord]; the ordinate for the upper 
surface is given as 


y = a(x — mC) + &x (29 


The general problem is similar to the zero thickness problem considered in Section (IV). In fact, the problem and 
the method of solution are so similar that only the major steps and the final results will be given for the case of the 
wedge. 


The first-order potential ®,' for the wedge is found to be 


; f M,*xf xoCf M2 | <—yp 
w * —-he + +e ae is ( ae ((¢) dé (30 
B §° B 33 B 0 


where the general first-order potential satisfies the homogeneous wave equations {O® = 0 + O(w*)| and the 
boundary conditions, Eqs. (9) and (11). 

By substituting ,' |Eq. (30)] into Eq. (5), the nonhomogeneous wave equation that is to be satisfied by the 
particular solution can be determined. Now let the second-order solution ¢*? be composed of a particular solution 
y and a correction potential 6. The particular solution that satisfies the nonhomogeneous wave equation 1s 


then, 
‘ VoM 94( ob 1) 2 a M,?( + 1) : a - 
y = —(a+ e€)’ z y(P)? + ala + 20 Vo] ; [ (f’)? dé + 
4B° . 1B J0 
Mi(y +1) w.. Mo * 1) Mo*(y + 1 be 
L ss nC | + ad(a + e) | : = xv(f’)? — : a (31 
4p° 26° _ B4 


As before y = ¥, = Oat x = 0; hence, the correction potential @ must satisfy the conditions 06 = 0, @ = @, 
OQatx = 0, and 


~ 
Ww 


(0,)y-0 = —(W)y=-0 + (a + ©)f’[(2!)r]y20 + [(Pa')y-o — (Pry), atef—xfe'] e 


By substituting Eqs. (30) and (31) into Eq. (32) and neglecting terms of order a’, the required correction potential 
6 is found to be 
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x vs 2 2 
rr\ 9 4,° met + 1) (M,! — M,”) . ’ rrr 
p= Vola t e)* [ Y )* dé [ 7 484 | 7 i Pine t ae % 


s— 98 M,2( _f. 1) r—yve XC r—yp 
ve hs,’ y | + a(a + e) ( —2 ff’ dé — : : / (f’)? de + 
0 he J : B2 Jo : : 


tallies 3 irene ie al 2Mo’ — 3Mo\ | Mo? — 2Mo! 
ye — YP), Cf (a) J? dn dge X | Cy + DV) 288 | Bs 
er M.? aes M,® + Mo 2M, (2M? 
ry 0 Jo / dé 32 + x a 5 ae i di (Y + 1) 2 36 ) a g + vit ( 8 = 
” , M? , 2M, M,° “— 2M’ i ee (2M? 
vc (Xx yB)/ + ad ; _ mC or i a. aes aC — + L} xCf + = ff’ de —f ( " ) - 
B? B82 Jo B? 


B* 
i — ik gO M,? 
wy2C2f’ + (x — vB)f’ 32 > + ae) )? dn dé — (x — yB) (f’)? dt Bt + 


M,’ , n Me ,, et ee 
T a —- as-—f +2 (f’)? dg (33) 
B? Jo 3° B? Jo : a 


[hus, the complete second-order velocity potential is 

?,” = ®,! + y + ial 
where ®,', y, and @ are given by Eqs. (30), (31), and (33). This completes the general solution. Now consider the 
y, /) to the oscillating wedge problem. 
yB, the potentiai &." may be determined 
| is found to be 


solution (in x, 
Setting {f(x — vB) =x - 
on the surface of the wedge [y = (a + ex — Ca 


y 2 a ‘a x 1 (y + 1) — 48° 
(, = (a -- €) —1 XC + (2 — M,’) rT (a + €) 
B* 18° 


a (x. . F Kel... ) : 
! = [y(Mo — 4Mo') — 3Mo® + 10Mo! — 10My?] — 5 [Moy + 1) - Ny (34) 


Then the pressure coefficient C, [Eq. (17) ] 


“= ff 


The pressure coefficient on the lower surface 


» of the wedge. 
The discussion of the second-order lift and 


lhis is the pressure coefficient for the upper surface 
is equal to, and 180° out of phase with, that of the upper surface. 


pitching moment will be taken up later. 


C) The Double-Wedge Airfoil 
[he problem now under consideration is the effect of a sharp bend or break in the upper surface of the wedg« 
\s a convenient reference term, this problem has been entitled the double-wedge problem, since the solution will 


certainly be applicable to such an airfoil shape. Consider the oscillating wedge as shown in Fig. 4. This problem 


must be solved in two parts, one part representing the solution for 


x/B< y< (x/B) + (C/B) 


ind one part representing the solution for 
8) + (C/B) < y < (x/B) + (1/8) (C+Z 


Of course, the previous wedge solution will take care of the solution for the diverging section of the wedge 
of a discontinuity in the slope 


$C). In addition, it will be necessary only to consider the effect, on the pressure, 


ofan oscillating flat plate; therefore, € will be set equal to zero. 
The problem for x < C is identical to that solved by the potential for the flow past the oscillating airfoil of zero 


thickness, Eq. (28). 
For x > C, the complete second-order potential is V2? (W.2 = yo! + y + 6). A new function g will now be de 
fined such that g =p, = Oatx = C. Later, for the specific problem of the broken plate, g(¢) will specifically be 


defined as g(¢) = g(x yB - C): thus, 0g/O0x = g’ and 0g/O0y = — 8g’. The first-order potential must satisfy 
the conditions 
Oy! = 0 
P ee eis M,* Mo*x ,, we 
(Yoy'),-90 = Volaf’ + 62°) ta] — mCf + 32 +1if- a I | Ist order (35) 


y.! = ye! and Yo,! = oz! at x = c 
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The last condition requires that there be no discontinuity in the potential or its . derivative anywhere in the flow 
In particular, there can be no discontinuities near the Mach line from the break in the plate. The second-orde, 
potential y? must satisfy the conditions 


Oyo? = Hy") 
| [Woy? - (Oy Ot] ( Vo + Yor”) a = Oo”, Ox 


2nd order 36 
Yo? = ge? and Yo? = go;z*atx = cy 


Again the boundary conditions allow no discontinuities in the flow. Now, note that this problem is similar ty 
the double-wedge problem in a and «. In fact, if 6g were replaced by ¢f, the problems would be identical excep 
for the requirements on ¢ and ¢,. However, it is apparent that this boundary condition would be satisfied j 
‘ both cases. The determination of the solution for x > C is much simplified by this fact. A first-order potential 


2.f als 9 "s vB 
didi a laf + a) + al 4 VMo?xf in XoCf _ Mo 1 ') (e) de y 
: B ; B° B B? B, AS0 i 


The method of solution which leads to the complete second-order potential is exactly the same as that previously 
employed for the flat plate or wedge airfoil. Since these equations are lengthy, they will not be presented; hoy 
ever, the equations and the steps involved in their development are presented in reference 1. 





is seen to be 


The pressure coefficient C, as determined from the second-order potential is given below. 
For x < C, the pressure coefficient Cp is given by Eq. (34) when « is set equal to zero. 


C 25 ( a )I C4 — Ve? | 1 E | (a )| v > Xe RB ) 
wy = —4, (a) — : Xo (2 — Mo) (a-) . —¢ t 0 
' i oe Vo B? 28 L/L 28° MS . 


where A = [(y7 + 1)Mo! — 487]/28", B = Mo?(y + 1) 1, and 





{ 


D = [y(4Mo* — Mo*) + 3M 1OMo! + 10.My*|/2° 


For x > C, the complete potential indicates that the pressure coeflicient Cp is again given by Eq. (35) if € is set | 


equal to —6. 


C 24 ( .) ( a JI C + Xx (2 Vv »| + | asin E | | a : | x D 4 xoC B|| 
~ = eh : Xo o— Ig C= = a — 0} ; ov 
a Vo B° 26 Vo ” tLe * 28° SS 


Thus the effect of a corner on the pressure has been determined. The third term [(A/28) (a@ — 6)?] of this equa 
tion is rather interesting, since it is the square of the instantaneous slope multiplied by exactly the same term first 
produced by Busemann for the second-order solution for the flow past a stationary surface. Thus, if the angle oi 
attack is not a function of time (@ = 0), this solution reduces to that given by Busemann. 


D) The General Airfoil 


Consider the airfoil shown in Fig. 1. From the results for the oscillating double wedge, it is apparent that bya 
series of m steps the pressure on the surface of an airfoil with ” breaks also is calculated. Then, independent of the 
shape of the airfoil ahead of the surface in question, the pressure at that point would be the same as that given by 
the equation for the pressure on the aft portion of double-wedge airfoil. Thus it is apparent that the pressure 
at any point is independent of the shape of the remainder of the airfoil and depends only on xo, Moy and the slope 
and the rate of change of slope at that point. Hence, the pressure on the general airfoil of Fig. 1 may be expressed 
inthe form: on the upper surface, 


ee 4( + “) («)] t+=>¢ Me) + ( + a) : ) 
mw BN dx =, — a a a 





( a )( 4 )/ x D 4 xoC B| { 40a 
VIN" * dx JLB 28° SI . 
on the lower surface, 


9 - : 2/ o | 
; “4( “) ( a )I - # , P | ( ra) A ) i 
ré = ‘ —_ ~_ Xo ve o M,? c —_ 
nm” B\ ol Ga * BS ‘de oa = | 
a oy x XoC | 
D B (40b 
( *\(< . dx | B T 2B f 
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Consider the specific problem of the modified double-wedge airfoil as shown in Fig. 5. 
= Cp, — Cp,) between the upper and lower surfaces of the modified doublé 


40b) the pressure difference (Cp; = 


wedge is found to be 


TWO-DIMENSIONAL 


AIRFOIL 695 


Lift and Pitching Moment Coefficients for a Modified Double-Wedge Airfoil 


From Eqs. (40a) and 


C med (1 or “) re i[2 249% (2 — M2) +2 =F D+ *B) C4 
rn : <COS W S Ww <Xp )~ eis ~— Mo an 4 ~, 19 P~ 
a ee a” oor Cc TS AS / 


where ¢ = J for << G 3: c 
Hence, the lift coefficient C),, is 


0 for C/3 < x < 2C/3; andft = 


—1 for 20/3 <x < C€. 


Cr = (4ay/B) {cos wt + (@/2) sin wt[2xoB? + (2 — Mo?) — (4€/9B)D}} (42) 


while the pitching moment coefficient Cy, may be expressed in the following form: 


2am} le @ ’ ‘ %e 
Cy = - cos wf | | 2x0 — A | + — sin of | 2[2 — M,? — — D} + 
B | 9,8 3 36 


Now consider the importance of the second-order 
terms in angle of attack. For the lift coefficient 
CL (4ay,/B) [P cos wt + (@/2) Q sin at | 


a plot of the function Q vs. Mach Number is presented 
on Fig. 5 for three thickness ratios (T/C = 0, 0.05, 0.10). 
It appears that the second-order effects in angle of at 
tack (see Fig. 5) are most predominant at the low Mach 
However, since the higher order terms in 
1.4, the relative impor 


Numbers. 
frequency are small for Mo < 
tance of the second-order terms in angle of attack is ap- 
proximately constant throughout the Mach Number 
range 1.2 < My < 2.0. 
For the moment coefficient where 
Cy (2ay/8) [R cos wt + (@/3)S sin of 

the functions R and S are plotted versus Mach Number 
for three thickness ratios (7./C = 0, 0.05, 0.10) on Figs. 
6Hand7. As before, the second-order terms have about 
the same relative importance throughout the Mach 
Number range 1.2 < My < 2.0. The second-order 
terms are especially important for the odd (sin) terms 
infrequency. In Fig. 7, it is seen that the second-order 
terms in angle of attack actually are of greater magni- 
tude than the first-order terms for most of the Mach 
Number range. It is this fact that causes the unstable 
tendencies, noted by Garrick and Rubinow for the case 
of airfoils of zero thickness, to be eliminated and the 
motion to be stable when the effects of thickness are 


considered. When the Mach Number is less than \/2 
for this case of x» = 0), the function S is seen to be 


positive for the case of a flat plate (7/C = 0). 
This term being positive corresponds to negative 
damping—i.e., an unstable motion. However, when 
the effect of thickness is considered (see Fig. 7), the 
function S is seen to be negative throughout the whole 
Mach Number range. This is further indicated in 
Fig.8. In Fig. 8, the curve labeled T/C = Ois the curve 


» 


€ 
Xo (os > + (2D 
3B 


-- Bs*)) — dna" | (43 
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Fic. 7. Comparison of the odd terms in the second-order 
solution for the pitching moment of a modified double-wedge 
airfoil at three thickness ratios 
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Fic. 9. Comparison of experimental and theoretical damping 


in pitch for a 7'/. per cent thick biconvex airfoil 


as presented by Garrick and Rubinow,'' which indi 
cates that an airfoil of zero thickness will be unstable 
for all points of rotation less than x» = °/; for some 
range of Mach Numbers. The curves for a thickness 
ratio (7'/C) of 0.03 and 0.0375 show that the region of 
instability rapidly disappears with any finite thickness, 
and in this case the instability disappears for 7/C > 
0.045. 


(V) COMPARISON WITH EXPERIMENTAL RESULTS 


In a paper” published in The Journal of the Roya 
Aeronautical Society, Beavan and Holder present the 
damping coefficients determined from supersonic tests 
of a 12-in. span, 7'/2-per cent thick, biconvex airfoil 
By utilizing Eqs. (40a) and (40b) and assuming a para 
bolic representation for the circular are airfoil, the 
pitching moment coefficient for a 7'/» per cent thick air 
foil is found to be 


2 0.05A 
i Om cos wt (: — 2xy — ) 
Pia B 
ao. 0.15D 
, Sinwt| 4 — 2M)? — 3 + Xxof (617° d+ 


0.15D 


- 0.07588) — custo" 44 
B 


The results given in reference 15 are given in the form 
where 


a 


of a nondimensional damping parameter M/ 


M . (Vo ou) (OCy Od) 


a 


(Concluded on page 704) 
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Partially Plastic Rotating Discs’ 


MELVIN ZAiDi 
Sperry Gyroscope Companyt 


ABSTRACT 


The solution is obtained for the flat rotating disc in the plastic 
region Strains are assumed to be small, and elastic compressi 
bility is considered in order to avoid discontinuities at the elas 
tic-plastic boundary. 

The distortion energy theory is used, and results are obtained 
for both strain and nonstrain hardening materials. 

Comparison is made to previous solutions based on the maxi 
mum shear theory and also to the numerical solutions obtained 
by Nadai and Donnell, all of which neglected elastic compres 


sibility. 
NOMENCLATURE 
D = flow function for tension 
E = modulus of elasticity 
H = 1+D 
[or’ — (0,'/2)} (3/2)|(o, — ot)/r| 
P; - 
or (o,/2) 
{(1/2) v|Tr 
Os = 
; C1 (or 2) 
rtats = (7) (KZ) (XE) 
g g So 2 So 
= radius to plastic-elastic boundary 
b = radius to outer rim of disc 
e = effective strain (+/ 2/3) xX 
V (e €2)? + (@2 — @3)? + (E €;)? 
é:,@; = principal strains (elastic plus plastic) 
é,é@z = radial, tangential, and axial strains (total 
r = radius 
Ni, 2 = r/a,r/b 
= effective stress (+/ 2/2) x 
V (51 Se)? + (te — $2)* + (s 51)? 
= yield stress 
» Sz, S = principal stresses 
,%,S2 = radial, tangential, and axial stresses 
Y = weight of disc per unit volume 
€ = Ee/so 
€r. €, € = Eer So, Ee, SO, Ee, So 
v = Poisson’s ratio 
) = radial displacement 
@ = §/s, 
Or, Of = Sr/So, St/Se 
w = angular speed of disc (radions/unit time) 


Primes denote differentiation with respect to radius r, 1, or 


(1) INTRODUCTION 


T HE PROBLEM OF THE THIN DISC is frequently en 
countered in the design of rotating parts. 


Received March 8, 1952 
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The 


engineer often desires information on the speed neces 
sary to cause complete yielding of the The 
stresses and strains in the elgstic region have long been 
known, but it was not until 1929 that solutions were 
In a paper by Nadai 


disc. 


presented for the plastic region. 
and Donnell,'! solutions were presented which were 
based on the maximum shear theory and on the dis 
tortion energy theory. The results, assuming the 
latter theory, led to numerical integrations and had to 
be presented graphically. 

In the contents of this paper the distortion energy 
theory is used, and an exact solution in the form of an 
infinite series is presented. A further exact but simple 
solution for ‘‘almost constant” yield stress is developed 
and extended to the arbitrary stress-strain law. The 
latter of the two solutions is to be recommended be 
cause of its simplicity and flexibility. Inasmuch as the 
stresses are statically determinate, any later assump 
tions with regard to strain do not affect the agreement 
between the stresses obtained from 
the article by Donnell! this 


paper. 


solutions for 


Nadai_ and and 


In solving for strains it was assumed, in the above 
mentioned paper, that the volume remains constant. 
In effect, this means that the elastic strains are ex 
tremely small as compared to those in the plastic re 
gion. It can be shown that in the partially plastic disc 
this is not the case; moreover, in almost any plastic 
elastic problem such an assumption will lead to serious 
discontinuities at the plastic-elastic boundary. To 
avoid such an occurrence, the more general stress 
strain relations are utilized. 

Since the paper of Nadai and Donnell presents the 
solution of strains based on the maximum shear theory, 
no comparison is made beyond the fact that the dis 
continuities the boundary. 
However, Fisher* has, by numerical means, evaluated 


exist at plastic-elastic 
the plastic-elastic coefficient and obtained strains based 
on the assumptions of Nadai and Donnell, using the dis 
tortion energy theory. A comparison is made on the 
basis of these results. 

Infinitesimal strains are assumed, and, based upon 
this, exact solutions for stress and strain at any point 
in the dise are obtained for a general stress-strain law. 
Since this often requires considerable labor in the form 
of the numerical solution of differential equations, it has 
been found necessary to introduce specified stress-strain 


laws. 
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For the case of the thin disc, 











S; = S;, q =e 
ma = Sy ey = €; 
3 =5s,=0, =e, altho 
man 
From the assumptions and the above definitions, one 
can write the equation of effective stress, strain, and the Case 
stress-strain laws in dimensionless variables as Pl 
ao = (V2/2)V (oe, — o,)? + o 2? + o?? { to be 
” repre 
e= (72/3) V (e, — €,)? + (e, — €,)? + (e é,)? an el 
5 the s 
e- = o, (1 + D) — a,[v + (D/2)] 6 
«, = o,(1 + D) — a, [v + (D/2)] (7 
é,= —\(e, + 0;) [v + (D 2)] (S Fron 
9/5 ” /< ince 
e = (2/3)o[1 + » + (3/2)D] f . 
(II) THEORY Eqs. 
Cp — i £ i ae 3 € have 
As the rotational speed of a thin disc is gradually in- = ee => (10 a 
; : i Bn Cr — 0; O1 — 3d; 20a ing 
creased, yielding starts at the center. The plastic- will! 
elastic boundary then progresses outward with in- The equilibrium of the disc is 
increase of speed until complete yielding of the disc 7 Fr 
ives! : oO, — 0, — 1. 2 (da,/dn 2») — Ti 2m 22 = 0 (11 a 
occurs, as in Fig. 1. it 18 
. — For compatibility of strains 16) 
It is assumed that (a) the material is homogeneous ; om 
and isotropic; (b) Hooke’s law applies in the elastic re- ri, 2 (de,/dn, 2) = & — & (12 of tt 
rion; () > sum of the total strains is . ' . , ‘ : 
gion; (c) the sum of the total strains is Introducing Eq. (10) and the stress-strain relations exatl 
_" (6) (7), and (8), one obtains will 
@ + 2 + €3 = E (Si + So + 53) to lis 
7 (dH/drn, 2) + P:\H = Q (13 
(d) the ratios of the principal strain differences to the [{ should be noted that 7, » depends on whether the 
corresponding stress differences are equal; and (€) radius of the plastic-elastic boundary (a), or the radius 
there is a unique relationship between effective stress of the outer rim (/) is used. It will be found that the 
and effective strain. Assumption (c) takes into con- — gojytion is most easily obtained in terms of r;, but, fora 
sideration the compressibility of the material even in meaning to be attached to the solution. it is al 
the plastic region. It is assumed that the total strains ways necessary to refer to a/b; here, ro will be neces 
in the plastic region are made up of two components. gary. In this connection it should be noted that, for 
Phese satisfy the condition that the change of volume the elastic solution, a = 7; = 0, leaving angular speed 
> ‘ . a Py ‘ 2c > . > ‘ MT) . ei a° . . 
equals (1 2v)/E times the sum of the principal (4) as an indefinite value. Since w is an important 
strains. Further, the plastic component satisfies the — yarjable and is often desired for comparison, it is neces 
relation that it has no volume change; hence, the elastic gary to introduce ro and, hence, 7>, which is never 
component must produce the entire volume change. zero ynlessw = 0. 7 
This leads to assumption (c). 
. None _ . 
The effective stress at a point is defined as (IIT) SOLUTION FOR STRESSES ) = 
6 
¢ = (4/2/2)X Case (1) Elastic 
° 
Ue = a) A Sa = Sa)? Fe ie —— Sy)® C1) Here D = D’ = 0, H = 1, H’ = 0. Using Es 
, : (11) and (13), the solution is easily obtained as : 
and the effective strain as , ‘ Eqs 
ee : o, = Ci — Core? — (72/8) (8 + vy)? (14 ical 
e = (V2/3)V (e1 — @2)? + (€2 — €3)7 + (e3 — e:)? (2) seri 
mS ao, = CQ + Carn? — (T2/8) (8v + 1)r.? (10 hej 
where 5), S2, S3 are the principal stresses and @, és, e; are er 
the principal strains. Eqs. (14) and (15) will be used to determine the stresses exp 
‘ } Soviet ; in the elastic region. We can prescribe o,, a, at 7 = I 
Assumption (e€) can be derivec m a tens > F 
lied I \ erived Irom a tension test a/b and o, at ro = 1}, and there are three unknowns, nT 
in the form . ' np ; B in 
Ci, Co, and a/b. The usual case wil! be take a,],,-1 = 
; or 
s = fle (3) O,and 
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PARTIALLY 


ore I, =a/b orp a =a/b 


Ft, I, =a/b = Ftp I, =a/b 


although the zero radial stress condition is by no means 
mandatory. 
Case (2) Elastic-Plastic (Nonstrain-Hardening Material) 


Here, the plastic region is assumed 
The case of 7; = 72 = 1 


Plastic Region. 
to bein the region 0 <7 < 1. 
represents complete yielding, whereas 7; = 0 indicates 
an elastic disc. Eqs. (4) and (11) completely specify 
the stresses; these are 


0,7 + ao,” — 0,0, = | (16) 
Te 2m Be. ri(da; ‘dry) - Tir: = (0) (17) 
From Eq. (17) it is easy to see that, for 7, = 0, o, = o,, 


since stresses are finite at this point. The solution of 
Eqs. (16) and (17) can be obtained numerically and 
have been utilized to check the accuracy of the follow- 
ing analysis. However, a far more powerful method 


will be formulated. 


From the Cauchy-Kowalewski‘ existence proposition, 
itis known that, having the values at a point, the Eqs. 
16) and (17) can be expanded in a Taylor Series which 
will converge in a sufficiently small region. The extent 
of this region is not easily determined, but a series of 
examples show that in all cases two or three expansions 
will cover the entire disc. Before proceeding, it is best 
to list the required derivatives of Eqs. (16) and (17). 


ov, = (1 11) (o, — o, — Tyr?) (1Sa) 
go,’ = (1/n) (o,’ — 26,’ — 2Ty7) (18b) 
go," = (I 11) ta,"" — 3e,’’ — 27)) (1Sec) 
Oe me (1/7) Ce,’ — 46,'"") (18d) 
og,” = (1/n) (of—' — no") (18e) 
| = g,* + 6;° — 60: (19a) 
0-= 2a,0,' + 20,0,’ — o,0,' — o,0," (19b) 
0 = 20,0,” + 20,'2 + 20,0," + 20," — 
g,0,/' — 2a,'a,' — 9,0," (19c) 
(0) = Gag ttl + 6a,’c,"" + ce Td + 
Ge.’ e,”" Z a ae Se. a: i oe: 6.” P=2 ee ed 
e (19d) 
Eqs. (18) and (19) are best solved for o;”, o,” numer 


There are now two likely starting places for a 
These are at 7; = O or 1, the former 


ically, 
series expansion. 
being discarded because of the resulting indeterminate 


expressions. Hence, the procedure is as follows: 


Let r, = 1 and assume a value of 7}, ¢,. Use this 


in Eq. (19a) to evaluate o,; then proceed to Eq. (18a) 


lor o,’. Using the known values of o,, ¢,, ¢,’, and 7) in 
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Eq. (19b), it is easy to evaluate o,’. In this manner 


one proceeds back and forth between Eqs. (18) and (19) 
to evaluate o,”, ¢,",atr; = 1. 
The solution is then presented in a Taylor expansion 


about this point: 


6, = o;(1) — 


(20) 


There still remains one boundary condition ¢, = | at 
r, = 0. In this connection it will be noted that we had 
assumed both 7}, ¢,(1), but these are not independent. 
Instead, they are related through the boundary condi- 
tion at 7; = 0, which will only be satisfied when the cor 
rect values of this combination are chosen. 

The values of o,, ¢, at 7, = 1 are then used as bound 
ary conditions for the elastic region, where this solution 
is carried out as described in Case (1). 

Fig. 2 presents some results and serves as a guide for 
estimations of 7), o,(1) for other solutions. From this 
figure it is seen that 7, = 2.43, start of yield, and 7, = 
3.35, complete penetration of yield. Therefore, the 
increase of speed for complete penetration is w, 
This value should be com 


1.175w,, or 17'/» per cent. 
pared with the 11 per cent based on the theory of maxi 
mum shear and the 17 per cent based on distortion 
energy obtained by Nadai and Donnell. 

Fig. 3 presents a graph of 7) vs. a/b, so that a definite 
radius of penetration may be assigned to each speed. 
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Case (3) Elastic Plastic (‘‘Almost’’ Constant co) 


By allowing the stress-strain curve to vary very little 
from that of a nonstrain-hardening material, a simple 
solution can be obtained. This method is capable of 
extension to more general types of materials. 

Let us assume that the stresses can be written as 


k 
o, = >) A,r,” — (1/2) Tyr? + C (21) 
1 


k 
Oo, = z Gn" — (1 2)7T in" +C (29 
1 


or, as a special case, 


CO, = Ar,” —(l 2)7Tir?® + C } 


te 


a, = Gr,” — (1 2)7Tir72= +C (* 


! 


Eqs. (21) and (22) will lead to a ‘‘perfect material” jj 
enough terms are taken, but the work becomes too 
great. Eqs. (23) and (24) lead to variations of ap 
proximately 2 to 3 per cent from the ideal; the latter 
is used in this paper. 

Using Eqs. (23) and (24) in the equilibrium Eq. (17), 


there results 


G- A(n+ 1) =0 (95 
Let us further specify o = 1, at 7, = 0, | and da /dr, = 
0 at 7, = 1, and we will automatically find do/dr, = 0) 


at r, = Oif we stipulate that » > 1. The three condi 
tions, when applied to Eq. (16), lead to 


= (26 


2(n> +n + 1)A2 + (n + 2) (2 — T))A + 
7{(1 2)7 | = al =( (27 


In(n? +n + 1)A2 + 2n(n + 2)A — 
Ti\(n? + 4n + 4)A — 47, — 27:7 0 (28 


Fig. 4 presents n, A vs. 7). Eq. (25) leads to G once 
these values are known. Thus, by assuming an arbi 
trary 7}, one can compute o,, a, in conjunction with 
Fig. 4 and Eqs. (25), (24), (25), and (26). Eq. (16) is 
then checked, and o vs. ” is plotted to see the degree 
of “‘perfect’’ material obtained. Typical results are 
shown in Fig. 5. By changing m slightly but still satis 


istying o = | atr, = 1, 0, a slight improvement can be 
obtained. This is also shown in Fig. 5. The elastic 
region is solved as in Case (2). Typical results for 
stress vs. 7; are shown in Figs. 6 and 7. A comparison 


with those obtained by methods of Case (2) shows al 


most zero deviation. 


Case (4) Variable « (an Inverse Method) 


Let us take the solution in the form 


Cr, = A 171" —_ | l 2)Tir,? + Fr, + s (29 
o, = Gyr" — (1/2)Tir,? + Ein + C (30 
also, we will assume that 7; = 7 or that the disc is com 


pletely plastic. This is not a necessary restriction but 
leads to some simplification. The boundary conditions 


are 


= 0 te, 31 


mi=1 71 = UV (0 
a, = af ' o, = Cf' 


Using these and the solution given above in Eqs. (4) and 


(17) leads to 
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PARTIALLY 


- = Lore 0 

F =m 2/3620" 

ky = 2F _ 
. (92) 

Ai = (G¢=9)/") — (F/n) 

7, = 2(A, + C+ F) 

Gi = Ai(n + 1) 


Eos. (32) can be solved in two ways: 
(a) Specify o(,=0), %,=1» 7 (=o) and solve for vari- 
ous parameters of n. 
(b) Specify 7, and only two values of the effective 
stress or its derivatives such as o(,~»)’, 0, =9) OT Oy=1) 
and then solve for n. 


Both methods are useful. Method (a) will indicate 
the relation between speed and the stress strain curves, 
while Method (b) will give a family of curves (c-e), 
which gives complete yielding at a specified speed. 

In all cases after the values of the coefficients are de 
termined from Eqs. (32), they will be used in con- 
junction with Eqs. (4), (29), and (30) to determine the 
curve o vs. ”%. Figs. 8 and 9 show such results for the 


/ 


case of o;,=9)’ = 0, —O0.2, respectively. In both cases 


the value o¢ —»)/o~=1;,)= 2 was used. It is quite feasible 
and not too difficult to obtain a complete set of such 
curves, so that any arbitrary stress-strain curve can be 
almost exactly represented and results obtained by 


interpolation. 


(IV) SOLUTION FOR STRAINS 


In order to solve for strains, it is necessary first to 
solve Eq. (13) for the modified plastic elastic coefficient 
H. The solution is 


Hel Pi 4 = S Que fPidrdy +. ¢ 
which is generally too complicated to evaluate. 


A second method is to use the same technique as in 
Case (2) for stresses. From Eq. (13), together with Eq. 


(17), we obtain 


(>, — a + (0, +o,’ + = iri) = 


, » 


| 0, a o >) 
(«, - a + (20,' + ms se Tin; )ur + 


3 ] 
1) HT = n( — ) 
9 9 


33b) 


(o,” +0, + 


0, 3 
(«, salker ) HH’ + (301 + = rn) H” + 


a1 + (¢,/" + ¢,°")H=0 
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, 


O, C, 3 
(«, — °) HW’ + (40, -— + = rin) HH!" + 


3 9 
(6o." +~ a," +- 1) H"’ + (40,'” + 


» 


0") A + (0 + 9," )H = 0 (38d) 


and all higher derivatives as required. From Eqs. 
(33a)—(33d) it is seen that, having the stresses and their 
derivatives at a point, // can be expanded about that 


point. This is, 


d ry sy 
(rn) =d) — “" W'd) + 


(d - r)? ii , 
a H''(d) — (34) 
where usually d = 1. Atd = 1, it is easy to see that 
IT = 1, since this is the plastic boundary. Now Eq. 


(34) will converge within a sufficiently small rgion, 
and it is possible, although not necessary, that this re 
gion will include the entire disc. However, to calcu- 
late the required derivatives for Eq. (34) is much too 


tedious a job. To simplify this procedure, it 1s neces- 
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e & | | | | | | | | | | d = 1, which is sufficiently accurate for this case. Ii 
BEeeea5 4 | | |_| ‘ et higher accuracy were desired, an expansion about d ~ 
| Bm + t +++ 0.2 would be sufficient. 
J pace ee 7 : Fig. 11 shows results of a series expansion for T ~ 
| a a [ io is 3.55 and almost constant o. Here, we see that for 
|_| Za ‘22 ie | ; 
| | : } } accurate results it was necessary to expand a second 
#8 | Sa ere time about d; = '/2. By expanding back to 7, = |. 
| | rl the error at 7; = O can be estimated. If the results are 
o6-- +++ ++ + | | a used for the expansion about d 1 for 1/2 <7, < Land 
| Ovsr | | » : , 
Prt eee TTT i — that about d; = '/, for0 < ™ <1/2, a maximum error 
— 0 ACCURATE NUMERICAL INTEGRATION | of about 4 per cent at 7; = O will result. 
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sary only to calculate the first few derivatives and values 
at a new point, where the convergence is sufficiently 
rapid. Using these values, derivatives are found about 
the new point, 
tinued over the entire disc. 

The solution can be divided into two main groups. 

Group (I).—Constant and almost constant o: The 
case of constant and almost constant o has been shown 


and in this way the solution is con- 


to have, for any degree of accuracy, the same a,, o,, and 
so. Because of this fact, it is apparent that they will 
have this same correspondence in strains. Hence, we 


will take, from Eqs. (23) and (24), 


o, = Ar,” — (1/2)Tir,? + C (35a) 

o,/ = nAr;— — Tir; (35b) 

a,” = n(n — 1)Ar,"—" — 7) (35c) 

ok = n(n —1)...(n —K +1) Ar," *® (35d) 
o, = Gr,” — (1/2) Tir? + C (36a) 

o,, = nGr;"—' — Tin (36b) 

oa,” = n(n — 1) Gr,"-? — T; (36c) 

o*&§ =n(n —1)...{n K + 1) Gr,"-* (36d) 


Group (II1).—Variable o: We obtain the values of 


the stress and derivatives from Eqs. (29) and (30); 
they are: 

o, = Ay" (1/2)Tin?® + Fr + C (37a) 

o,/ = nA"! — Tin, + F (37b) 

o,” = n(n — 1)Ayr,"-? — T, (37c) 

go~= =n(n—1)...(n K + 1) Ayr’ (37d) 

o, = Gir” — (1/2)Tin? + Ain + C (38a ) 

o, = nGir,"—'! — Tin + E, (38b) 

a,” = n(n — 1) Gyr”? — Ty; (38c) 

of} = n(n —1)...(n — K +1) Gir (38d) 


Thus, /7 can be found by use of Eqs. (35) and (36) or 
(37) and (38), together with Eqs. (33) and (34). 
Eqs. (6) to (9) inclusive, the strains are computed, and, 
in the case of variable o, the o-« curve can be established. 

Fig. 10 compares D obtained by numerical integra 
tion and by series expansion for constant and almost 
The series expansion is about 


From 


constant o, 7; = 1. 


Fig. 12 shows a comparison of the strains for ty 
individual numerical integrations with constant ¢ fo; 
the series expansion based on approximately constant 
o and for the theory of maximum shear neglecting elas 
tic strains in the plastic region. Here, the discontinyj 
ties of the latter theory become apparent. 
especially that of e,, are greatly in error. 
that the results of the series expansion agree with the 


The strains, 
We also note 


numerical integrations. 

Typical solutions for a strain-hardening material are 
13 to 16, inclusive. They were only 
A comparison with 


presented in Figs. 
investigated for Gjnside = 2outside: 
Figs. 8 and 9 will lead to the conclusion that the more 
easily obtained stress-radius diagram can be of great 
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value in determining a good approximation to the re- 


quired stress-strain curve. 
(V) EXAMPLE 


Given: Disc: (1) 24 in. radius, flat; (2) steel (4 
0.283); and (3) yield stress (so) = 50 X 10%. 
Required: (1) Speed necessary to catise a plastic 
region of 2S-in. diameter. 
the plastic region to extend throughout the disc. 


Stress distribution in the disc with a plastic region of 


(2) Speed necessary to cause 


(3) 


28-in. diameter. (4) Strain distribution in the disc 
with a plastic region of 28-in. diameter. 

Solution: (1) a = 28/2 = 14 and b = 24 in.; 
therefore a/b = 0.5833. From Fig. 3 we find 7) = 
2.94. Since 

T2 = (y/g) (w?/so) b* 
we have w = 590 rad./sec. or w = 5,640 r.p.m. to cause 


a penetration of 14-in. radius. 
(2) For full penetration, a 


Fig. 3, T> 


banda/b = 1. From 


3.35 and, as above, 


635 rad. = 6,050 r.p.m. 


om 
(3) 7, = (a®/b®)T», or for 14-in. penetration a/b = 
0.5833, 7; = 1. From Fig. 4, A = 0.228, n = 1.82, and 


Re 


using these values in Eq. (25) leads to G = 0.642. 


turning to Eqs. (23) and 


have 
S,/So = 0.228r,!-82 — (1/2) (1) nm? +1 (Ex 1) 
$1/Sq9 = 0.6427,!-8? = (1 2) (1) r3? + ] (Ex 2) 


(24) for the plastic region, we 
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O< 4 < I. 
(15), which give 


S,/So9 = Ci — Core? — (2.94/8) (3.3) ro? (Ex 3) 
$,/So = Cy + Core? — (2.94/8) (1.9)ro? (Ex 4) 
0.583 < re< 1. Now, atr; = 1, re = 0.583, we have 


(S; S50) plastic (S; 50) elastic 


(S;, S50) plastic (S1/ So) elastic 


using these facts in the above equations results in 


C,; = 1.259 and C. = 0.041. Thus, 


For the elastic region we use Eqs. (14) and 
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S,/So = 1.259 — 0.04172.-? — 
1.259 + 0.0417r2~? — 0.699722 


5;/So = (Ex 6) 

We can check Eq. (Ex 5) by the fact that, at ro = | 
Ss, = 0.008s9, but we expect s, = 0. Considering that 
the values of 72, n, and A are read from graphs, the 
Recalling that r; = 


accuracy is high. r/Q, ro = p/h 


we use Eqs. (Ex 1), (Ex 2), (Ex 5), and (Ex 6) to fing 
S;, 5; VS. ras Shown in Fig. 17. 
(4) The strains can be obtained from Fig, 12 ; 


terms of €,, €;, €,, VS. 7. 
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A Second-Order Solution for an Oscillating Two-Dimensional 
Supersonic Airfoil 
. (Concluded from page 696) 


Hence, the theoretical second-order solution for the 


airfoil rotating about its midchord (x) = ! 


to be 


l 0.15D 
M, = 2— M*? - — 0.075858 
5 8 


(45) 


The damping coefficients as given by the experimental 
investigations, linearized theory, and the present paper 
[Eq. (45)]| are presented on Fig. 9 as a function of Mach 
Number. 
reproduces the trends of, but is somewhat lower in 
Some of the 


As may be seen, the second-order? solution 


magnitude than, the experimental data. 
discrepancy is probably due to the effects of flexure in 
the model which would tend toward increasing the 
magnitude of .1/,. 
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On Heat Transfer over a Sweat-Cooled 
Surface in Laminar Compressible Flow 
with a Pressure Gradient’ 


MORRIS MORDUCHOW? 
Polytechnic Institute of Brooklyn 


SUMMARY 


A simple expression is derived for the normal injection velocity 
distribution theoretically required to maintain a given uniform 
temperature along a porous surface in the laminar boundary 
layer region of a compressible flow with a given velocity distribu 
ses outside of the boundary layer. This expression is valid for 
any given free-stream Mach Number but is based on a Prandtl! 
Number of unity and on the assumption that the viscosity co 
efficient varies linearly with the temperature. By using the 
Dorodnitsyn type of transformation, the variation of fluid prop 
erties even in the case of zero Mach Number is taken into ac 
count. This study is of particular practical interest in connec 
tion with the sweat-cooling of turbine blades and of airfoil surfaces 
in high-speed flow. The method of analysis consists of applying 
the Karman-Pohlhausen method to both the momentum and 
energy boundary-layer equations and of using an additional heat 
A closed 


form approximate solution of the equations is then derived. Nu 


balance equation, involving the coolant temperature 


merical examples for flow in the immediate vicinity of a stagna 
tion point and for a typical type of flow over a turbine blade are 


given 
INTRODUCTION 


A PROBLEM OF PARTICULAR INTEREST in connection 
with the sweat-cooling of turbine blades or of air 
foil surfaces in high-speed flow is that of maintaining 
a given uniform wall temperature. This can be done 
by designing the porosity of the surfaces so as to ob 
tain a certain distribution, along the walls, of normal 
mass flow of the coolant. The main purpose of the pres 
ent analysis is to derive this required distribution 
The 


required mass flow distribution will, in general, be a 


theoretically for the general case of laminar flow. 


function of the coolant temperature, the desired uni- 
form wall temperature, the free-stream Mach Number, 
and the velocity distribution outside of the boundary 
layer, which depends on the shape of the surface profile. 
The 
coolant will be assumed to be the same as the external 
gas flowing over the surface. Radiation effects are 
neglected, while the temperature of the coolant is as- 
Centrifugal and 


All of these factors are here taken into account. 


sumed as constant along the surface. 
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Coriolis effects due to possible rotation of the airfoil 
surfaces will not be considered. In view of the results 
of Fogarty,' this may not cause any serious errors in 
actual application of the theory to rotating blades. 

Heat-transfer analyses have already been carried out 
for a number of special cases. Yuan** has investi 
gated in detail the heat-transfer characteristics for in- 
compressible and compressible flow over a sweat-cooled 
flat plate, using the Karman-Pohlhausen method. 
Klunker and Ivey' have also investigated this case. 
They used various simplifying assumptions but in- 
cluded, in their heat-transfer calculations, effects of 
radiation and (if the coolant is a liquid) of vaporization. 
Eckert,® using the transformation methods of Schlich- 
ting® and of Falkner and Skan’ has analyzed the case in 
which the velocity outside of the boundary layer is pro 
portional to a power (especially the first) of distance 
along the wall from the stagnation point. Eckert’s 
analysis was based on incompressible flow and con- 
stant fluid properties. 

The method to be used in the present investigation is 
essentially the same as that of references 8 and 9 and is 
based on an application of the Karman-Pohlhausen 
method to both the energy and momentum boundary 
layer equations for compressible flow. In reference 8, 
however, a normal velocity at the wall was not con- 
sidered. In reference 9, a normal velocity at the wall 
was considered, but a heat balance involving the cool 
ant temperature was not taken into account. Hence, 
in reference 9 the coolant injection velocity and the wall 
temperature were considered as independent parain 
eters, instead of as being related through the temper 
ature of the coolant; the separate effects of the injection 
parameter and the wall temperature parameter on the 
boundary-layer characteristics could then be investi 
gated there in some detail for general cases. In the 
present analysis, a heat balance equation is added to 
the boundary-layer equations, and the primary em- 
phasis is now not on the boundary-layer characteris- 
tics as such but on the heat-transfer properties at the 
wall and on the injection velocity distribution required 
to maintain a given uniform wall temperature. 

By making use of the Dorodnitsyn type of transfor- 
mation for compressible boundary-layer flow and by 
making the approximating assumption that the ratio 
of dynamical to thermal boundary-layer thickness is 
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constant along the flow, it is possible to obtain a simple 
approximate solution to the ordinary differential equa- 
tions developed here. For any given main-stream 
velocity distribution and Mach Number, this solution 
permits a rapid determination of any pertinent heat- 
transfer characteristics and of the injection velocity 
distribution required for a given uniform wall temper- 
ature. 

The results obtained here are based on laminar flow, 
which is the type to be expected at least near the lead- 
ing edge of a surface where the pressure gradient is 
favorable. The results are also based on a Prandtl 
Number of unity and on the assumption that the vis- 
cosity coefficient is proportional to the absolute tem- 
perature. 

As numerical examples, results for the flow near a 
forward stagnation point are given in some detail, since 
this is the type of flow which may be expected in the 
immediate vicinity of the leading edge of a round-nosed 
object in subsonic flow. Numerical results are also 
given for what might represent the flow over the suc- 
tion side and over the pressure side of an actual turbine 


blade. 


SYMBOLS 


ay = coefficient of 7; in velocity profile 

a), a2 = constant ‘“‘average’”’ values of a; and a,” 

bd; = coefficient of 72 in stagnation enthalpy profile 

c = coolant-temperature parameter, defined by Eq. (16b) 

Cp = specific heat of gas at constant pressure 

g ratio of dynamical to thermal boundary-layer thick- 
ness in (x, ¢) plane (g = 67/6 = x/X) 

H = stagnation enthalpy, defined as the quantity 


[(u?/2) + cpT] 
h = H,/M,; for a Prandtl Number of 1, 2 = 7) 


of actual wall temperature to equilibrium zero 


T, (ratio 


heat-transfer wall temperature ) 


k = coefficient of heat conductivity of gas 
l = characteristic length 
M = Mach Number 
Re = Reynolds Number (p,,%,,//n., ) 
te = absolute temperature 
Te = temperature of the coolant 
t = variable defined by Eq. (7) 
u = velocity of flow along surface 
v = velocity of flow normal to surface 
x,y = distance along and normal to surface, respectively 
a = ratio of specific heats (¢p/cy); y = 1.4 for air here 
6n, 52 = dynamical and thermal boundary-layer thicknesses, 
respectively, in (x, ¢) plane 
m, 72 = (t/5n) and (t/5,2), respectively 
k = (5,/1)Re'” 
d = (62/1)Re'/? 
v = kinematic viscosity (u/p) 
é = x/l 
¢ = [(potr/pu,, )Re'/*], injection parameter 
Subscripts 
0 = values at wall—e.g., 7», vo, po 
1 = values at local outer edge of boundary layer—e.g., J/,, 
T1, Uy 
© = values at a convenient reference point outside of bound 


ary layer—e.g., in undisturbed flow at infinity in ab- 
sence of shock waves 


A prime (’) denotes differentiation with respect to & 
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Basic EQUATIONS 


The momentum and energy equations for compres. 
sible two-dimensional flow in the laminar boundary 
layer of a gas with a Prandtl Number (uc,/k) of unity 
can be expressed in the form 

Ou 4 Ou duy fe) ( Ou 
i= pv = pil + m ) 

Ox Oy dx Ov \ Oy 
ol oH a) ( =} 
d pu = fm 
Ox Oy oy Oy / , 


pl 


pl 


The continuity equation is 
[O(pu)/Ox] + [O(pv)/Oy] = O 9 


Since the pressure is assumed as constant across the 
boundary-layer thickness, the ideal gas law implies 
p/p. = 1,/T { 
It will be assumed here that the coefficient of viscosity 
is proportional to the temperature—.e., 
M/Mi = g 7; (3 
The heat-balance equation, stating that the heat 
transferred by the hot gas to the wall is absorbed by 
the coolant, can be written as (cf. references 2, 3, 4, and 
5) * 


[R(OT Oy) lo _ pooCp( To ‘il Te) 0 


Introducing the Dorodnitsyn variable ¢ defined by 


y = / (T/T) dt 7 


« 


integrating both Eqs. (1) and (2) with respect to ¢ from 
t = 0 to t = 6, or 62, whichever is larger, and using 
Eqs. (3), (4), and (5), the following equations expressed 
in nondimensional form, are obtained, under the bound- 
ary conditions of u = 0, v = v(x) at the wall and smooth 
transition of the velocity and stagnation enthalpy 
profiles to the local main-stream values at the outer 
edges of the boundary layer: 


iy ‘ . pL oy es, u 
Fixx’ + x? | FP, + FF, + P, = x 
Pi u te 
p Vy 
Kg + f 5 
pi v 
P é “ur > Uy! p’ 
FX - F, + FP, —T 
1 pi / 
u p V} 
rect — h) + F; d 
in Pi V 


* In references 2 and 3, Eq. (6) was used in integrated form over 
the length of the plate, so that the heat balance was satisfied ther 
only in the average rather than at each point along the plat 
Consequently, the constant surface temperature appearing 1 
references 2 and 3 must be interpreted as a kind of average tem 
perature over the length of the plate, since according to Eq 6 


in the present form the surface temperature will actually vary 


along the plate if the injection velocity is uniform. 
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where plane: At y = 0, 
tas’ u u/u, = 0, A/i, =t 
F, = ( ( il — dn 
0 \M My Ou du, oO Ou 
F, = [(0/Om) (w/t) Jo pe ay), ae * Lay\" dy/ | 
Fy = [(0/On2) (7/FD)) Io : : : 
er oll o/ Ol 
: l 7 u Povo - Ke 
Fy = 2F, + —— dy. + Ov /o Ov\ Ov 0 
g/0 7) uy 4 : : 
— 
/ (; _*% ) dm (g> 1) The last two conditions above follow from the original 
Jvo\T,  W/ ny/m=1 - partial differential Eqs. (1) and (2). 
=n OR 4 [ ’ i _u ) in 4 At the outer edge of the boundary layer, the follow- 
ee ors Sy —< om * (10a) ing conditions are satisfied by the profiles (11) and (12): 
. : I 


(g < 1) 


[ rT =) 
g Jo e i WS u/us ‘aia 
[()0 — 7 dy. (g = 1) 
Jo \uy IT; ae eee 
“0/1 I] 

ie J (‘JC . = re 


From the definition of the stagnation enthalpy J//, it 


u\?/y — 1 
we) —(“)(™S) ae 
My ( 2 ; 


(10) 


follows that 
1 (= )( > 
1 + 
T, HH, 2 


Eqs. (S) and (9) are essentially extensions of the von 


Karman integral equation to the momentum and 
energy laminar boundary-layer equations in compres- 


sible flow. 


The following fourth-degree profiles in m and m2 are 
now assumed for the velocity and the stagnation en- 
thalpy : 


u/Uy = aym + (6 — 8a,)m" + 


(3a, = S)m.” + (3 — a;)mi3 (11) 


H/H, h + bine + (6 — Gh — 3d1)n2? + 
(Sh — S + 30,)no*® + (38 — 3h — by)no* (12) 
where 
1 my’ y-— 1 v 
* i aide h (1 +- ui) a 
Ou ' 2 VY} 
ay = 
1+ (1/6)xg(T../T)) 
a = 2 hA)y/’l+ a HAG(T T;)] (14) 


lhese profiles have been chosen to satisfy the following 
boundary conditions at the wall in the physical (x, y) 


Att = On (m = £}, 


ui ' O(u/u)) - O"(u/ Uy) 

uy ; On) On" 
Att a 6,2 (ne = 1), 

H | OH/th) _ o*(H/th) 

Tt; . ; One On” 


Since the Prandtl Number has been assumed as unity, 


the equilibrium wall temperature 7), for zero heat 


transfer (b,; = 0) will be 


Moreover, 


Hy tpl 


= T,/T. (15b) 
TI, + Col 1 te 


k= 
(1) 2) 


Thus, for a Prandtl Number of unity, / is the ratio of 
actual wall temperature to the equilibrium zero-heat 
transfer wall temperature. | is that of a 
heat-insulated wall, while the cases h < 1 and h > | are, 
respectively, those of an artificially cooled and an arti 


The case h 


ficially heated wall. 
Eq. (6), with the use oi Eqs. (10b), (12), (14), (15a), 
and (15b) and with the relation \ = «/g, can be written 


in nondimensional form as a quadratic in (gx). Solv- 
ing this equation for (¢x) yields 
gx = Co(T,/T.) (16a) 


where C is a constant defined as 


C a(qi4 ‘| l—h 4 ') ” 
= «) D) 
{ 3LA—(T./T.) 19 = 


Thus the required injection parameter g(£) varies in 
versely as the dimensionless boundary-layer thickness 


xin the (x, ¢) plane. 
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With the profiles (11) and (12), the expressions for the /; will be 


F, = 0.1143 + 0.00953a, — 0.00397a,? 
F, = Q 
Fy = b 


l 0.4 _ (by y- 
y «oF, +@011 ~ +. ( h + 0.05 ( —a)+ = 
g g g 2 
0.4 0.4 0.05 _ a 
jot — + ( )n + ( Jo — 0.0405a, — 0.003071" 
g ‘2 g 


I 


F, = (1 — h) [0.4 — (04g) + 0.22862* — 0.1429g4 + 0.0286g°] + ai(1 — h) (0.05¢ — 0.0429¢* 4 

0.0286g4 — 0.00596g°) + b,(—0.05 + O.1g? — 0.1143g? + 0.0536g! — 0.00953g*°) + 

ayb; (— 0.016672 + 0.0214g* — 0.0107g4 + 0.00199g°) (g < 1) 
’ 0.2286 0.1429 0.0286 0.1 0.1145 0.0536 0.009538 
F, = (1 — hy ( — — ee ) +a — ny ( — — + -- 

g g" e Ff. g 4 r * 
f =: 0.0429 0.0286 0.00595 0.01667 0.0214 0.0107 0.00199 
by (- . +. = = ; + ah; | — t = ” ah ; (g 
 - g' g g g* g" g 


With Eq. (16a), the expressions (13) and (14) for a 


and ); become 


I uy’ — | 
4-3 (: = uy’ ? 
Ou. Y V} 


= (18) 
6 Cg 


bh = 211 — h)/[1 + (1/6)C] (19) 

From the relation (u?/2) + c,7 = constant, and the 
assumption of isentropic flow, the following relations 
can be derived for the flow, outside of the boundary 


layer :* 
7,\/T. = 1+ [(y — 1)/2])M.2[1 — (u/u.)?] (20a) 
My = 1. (h/#.)* (11/1) (20b) 
p/p. = (7;/T.)-” (20c) 


The problem in this investigation can now be formu- 
lated thus. For a given constant ratio (4) of wall tem- 
perature to equilibrium temperature, a given value of 
the coolant-temperature parameter (C), a given Mach 
Number (M.), 
(u1/u..)(&), itis required to solve Eqs. (8), (9), and (16a) 
for x(£), A(E), and ¢(€). 
of these equations will be derived in the next section. 


and a given velocity distribution 


A simple approximate solution 


It should be noted that in the present analysis the 
case of zero Mach Number is not equivalent to incom- 
(i.e., Zo 7) the 
density p, as well as the temperature and viscosity 


pressible flow, since when h + | 


* In the case of a shock-wave at the leading edge, the quanti 
ties with subscript ‘‘ ©’’ refer to a convenient reference point be 


hind the wave. 


coefficient, will still vary across the boundary-layer 
thickness in accordance with Eqs. (4), (5), (10b), (12), 
and (14). 


(GENERAL APPROXIMATE SOLUTION 


To solve Eqs. (8), (9), and (16a) in a simple manner, 


the following two approximating assumptions aré 
made: 

(a) The ratio (g) of velocity (dynamical) to stag 
nation enthalpy (thermal) boundary-layer thickness 
in the (x, ¢) plane is constant along the flow (although 
it may be different for different types of flow). 

(b) In the expressions for F; and F» (and only there 
the quantities a; and a,* may be replaced by constant 
‘average’ values a, and a,’. 

Assumption (a), aside from the fact that it may ap 
pear reasonable physically, is based on noting that it 
is exactly valid for the flow over a flat plate (m/u. = 
|) with a uniform wall temperature, regardless of the 
Mach Number and of the injection velocity distribu 
tion. In that case, in fact, Eqs. (8) and (9) imply 


that g = | everywhere. Moreover, as will be shown 
subsequently, this assumption is also exactly valid for 
flow at zero Mach Number near a_ stagnation 


point. 

Assumption (b) is justified by the fact that the a 
terms appearing in F; and F, [Eqs. (17)| are relatively 
small and that, consequently, errors made in evaluat 
ing these terms will have only a relatively small effect 
there. As a check on the accuracy of assumptions (a 
and (b), the cases m/v. = 1 + ~ and m/u | =% 
were solved numerically and also by the use of thesé 
assumptions. The details are shown in reference 9; 
satisfactory agreement between the two sets of solutions 
was obtained. 
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By substituting for g, as given by Eq. (16a), into Eq. (8),* using the expression (18) for a; in the right side of 


Eq. (8), and making use of assumptions (a) and (b), as well as of relations (20a), (200), and (20c), Eq. (8) becomes 
the following linear differential equation in x*: 


cy u;" awe 2 u 2 aia i : 2 
(x?)’ + (J, — JoM,?)«x? = — — Cg + (21) 
it; PF, i 7 ] + (Cg 6) 


where Fj, J1, and J, are constants given by 


F, = 0.1148 + 0.00953a, — 0.003974," 


2 (0 , (: ') ., 0.1 j jO4 1/6 0.1 j ane ) 
J, =4 . = : a = : = : 0; 
= T F, g gil + (C/6) lg 1 + (Cg/6) g{1 + (C/6)}h ™ 


y¥-—1\/... 0.4 0.1 ; 
j=2- ( 2 (0.514 — + * — 0.0405a, 
i, g g| 1 - (C6) 


" 0.4 1/6 0.1 \ 
0.00397 a,2 + h = : s 7R\1¢ 
lg 1+ (Ce/6) gfl + (C/6))f 


(22) 


An integrating factor J of Eq. (21) is, by the use of relations (20a) and (20b), found to be 
[ = (%/uU )/ (7; ‘ig 7 i 


[he solution of Eq. (21), satisfying the condition of zero or finite? boundary-layer thickness at the leading edge 


t = ()), 1s, then: 


¥ Iimt par pa \ Geert se 

A [ (uu. ) (Del F ad dé 
= S70 9 
(2 


» 
x? = — | Cg + 
aI S 1 + (Cg/6) (u,/u.)"' (T,/T.) ; 


The distribution of the injection parameter ¢ required to maintain a uniform wall temperature will, according 
to Eq. (16a), be 
g(f£) = Cel (7) T.)/« (24) 
where «(£) is given by Eq. (23). The magnitude of the required injection parameter ¢ is evidently proportional to 
Ce(Cg + (2/[1 + (Cg/6)]}) 
and thus will increase as.(Cg) increases. 
For a given velocity distribution (# ..)(€) and a given ratio (/) of wall temperature to adiabatic equilibrium 
temperature, { the effect of Mach Number is given simply by the (7\/7’..) terms in Eqs. (25) and (24). From 
q 


Eqs. (23) and (24) it follows that, whereas in the flow over a flat plate (w/u.. = 71/7 |) the injection param 


eter ¢(£) required to maintain a uniform wall temperature ratio h will depend only on the parameter C and not on 


h itself, nor on the Mach Number ./.., the required injection parameter in the case of flow with a pressure gra 
dient will depend on the value of h (which appears in J; and J:) and of V/.,, as well as on C. 

In cases in which x = 0 at the leading edge (£ 0), x will increase with £, and ¢g(£) according to Eq. (24) will 
then diminish along the wall from an infinite value at & = 0. Such is the case for flow over a flat plate or for 
supersonic flow over a thin sharp-edged airfoil. In the case of a blunt-nosed object in subsonic flow, however, the 
leading edge will be a stagnation point, and « will then have a finite value there. In that case, Eq. (24), with 


) 


application of L’H6épital’s Rule to Eq. (23), implies that 


—_ Ce (t ) 2(y—1 (“) ee 
rang = (Fs351/2) Fe <0) 
= = 2 1 t=0 u r=() 

6) 


Co is. - » g 
| . 1 + (Cg 


Thus, for flow near a forward stagnation point, ¢ will be finite at the leading edge and will be proportional to the 


square root of the velocity gradient (1 u..)’ there. 


* ‘The use of Eq. (16a) by means of this substitution actually simplifies Eq. (8) [and would also simplify Eq. (9)] considerably. 
t Incase u,/u,, = Oaté = 0 
t For a given ratio (7)/T.,, ) of wall temperature to free-stream temperature, a Mach Number effect will also appear in the value of 


i, Since, for a given value of (7\/T7..), A diminishes with increasing Mach Number 
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By means of Eqs. (23) and (24), the injection velocity 
distribution required to maintain a given uniform wall 
temperature ratio / with a given value of the coolant- 
temperature parameter C can be quickly calculated for 
any flow with a given velocity distribution (w/w...) (&) 
and given Mach Number M.. The calculations will 
involve, at most, numerical integration (e.g., Simpson’s 
Rule). The present results can be used for the case 
of heat transfer with zero injection velocity by putting 
C = Oin Egs. (22) and (23) and by keeping / arbitrary. 
A value of g in any givefi case can usually be obtained 
aes. 2 «xk =} = Oat 
0 and 2,’ is finite there, 


from leading-edge conditions. 
the leading edge* while u + 
Eqs. (8), (9), and (16a) will be ‘found to imply g = 1 
there; this value may then be assumed for the entire 
flow (cf. reference 9). On the other hand, if the leading 
edge is a stagnation point, the value of g to be chosen 
is usually less than | and can be obtained in the manner 
shown in the first numerical example below. 

The Nusselt Number Nu, which is a measure of the 
local rate of heat transfer at the wall, can be calculated 
by means of the expressions 


eh g ' 
as 2 z Re (26) 
T./«{1 + (C/6) 


From Eqs. (26) and (23) it will be found that an increase 


: [R(OT/Oy) |o 
Nu = — — 

R(T, — To) 
in the parameter C—-i.e., an increase in the value of 
(1 — h)/(h — (T./T,)] 
selt Number. 


causes a decrease in the Nus- 


Eqs. (24) and (26) imply the following simple relation 
between the local Nusselt Number Nu (£) and the injec- 
tion parameter ¢ (&) required to maintain a given uni 
form wall temperature: 


Nu (&) Re~ “*/g (&) = 2/C[1 + (C/6)] (27) 


Eq. (27) states that, in general, the Nusselt Number 
at any station & will be proportional to the injection 
parameter there, the proportionality factor being a 
constant dependent only on the coolant-temperature 
parameter C. This relation will be valid for any Mach 
Number and any main-stream velocity distribution 
(w/u.) (€). Thus, the Nusselt Number distribution 
and the injection parameter distribution along the 
wall will, in general, be exactly similar if the wall tem- 
perature is uniform. 

The local skin-friction coefficient C, follows readily 
from the expressions 


0u/OY)o u fp 
Cy OY ( N(. ; ) = Re (28) 
pats" ors ar « 


where a, (£) is given by Eq. (18). Because of the x? 
term in the numerator of a, it will be found that the 
value of the parameter C has a smaller effect on the 


value of C; than on the value of Nu. On the other hand, 


* Such will be the case, for example, if 


= 0), 


U,/U. = Co + OE + cot? +...(& 
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for a fixed C, a decrease of the wall temperature ratio 
h will tend to lower the skin friction in a favorable pres- 
sure gradient (#' > 0) but will have only a small effect 
on the Nusselt Number. 

Velocity profiles can be calculated by means of Rg, 
(11) in conjunction with the transformation (7) to the 
physical (x, y) plane. The temperature field can pe 
similarly obtained from Eqs. (10b) and (12) 


NUMERICAL EXAMPLES 


kixample (1) Flow Near a Stagnation Point.—The 


case (;/u..)(~) = mé, where m is a positive constant, 
represents the flow in the immediate vicinity of a for 
ward stagnation point, such as at the leading edge of g 
round-nosed object in subsonic flow. It may be oj 
particular practical importance to analyze such a flow, 
since, thereby, leading-edge conditions for the subsonic 
flow over most blunt-nosed objects can be obtained. In 
particular, as will be illustrated by the second example, 
the value of gy can thus be determined. 


For “,/u., = mé, the change of variables 


” 1 1/o 
& = m& k= Km", X, Am’, ¢1 = om 


. (9Q 


will leave Eqs. (8), (9), and (16a) unchanged, except that 
£ is now replaced by &, « by «, A by Ay, and ¢ by ¢ 
However, 4;/u.. = & now, while the constant m no 
longer appears in the equations. For u,/u, f 
and M., = 0, it can be easily verified that Eqs. (8), 
(9), and (16a), or Eq. (23) with a; = a, and a,? = a’, will 
be satisfied exactly by the simple type of solution «x, 
A, and ¢ = constant. In this case, in fact, the follow 
ing simultaneous algebraic equations in «x; and \, result 


from Eqs. (8), (9), and (16a): 
Fox" = Cg + ay 30) 
Fd? — i(1 _ h)C) of by 31 


where 2, Fy, a, (which contain the unknowns x; and 
\i), and ; are given by Eqs. (17), (18), and (19). 

Eqs. (30) and (31) can be solved numerically for 
given values of Cand h. The values of ¢ are then cal 
culated from Eq. (24). 
Nusselt Number are given in Figs. 1 and 2. From 
Eqs. (29) and (26) it is evident, as would also follow 


Solutions for xm, g, ¢i, and 


from Eq. (25), that ¢ and Nu are each proportional to 
m’*. Thus, the velocity gradient m plays an important 
role here. 

E:xample (2) Flow over a Turbine Blade.—Fig. 3 repre 
sents, at least qualitatively, the velocity distribution 
(W1/U)(€) over the forward part of a turbine blade 
airfoil in subsonic flow. Both the pressure and suction 
sides of the airfoil are represented, and the length / is 
here taken to be the chord length of the blade profile 
The dimensionless distance £ is measured along the 
surface from the forward stagnation point. In_ this 


example, the value 1/7. = 0 will be assumed. 
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HEAT TRANSFER OVER A 

For this example, Eq. (23) can be used. Since M/, = 
ghere, Te = 11 = Tao. Let h (T/T) = '/3 and 
T/T, = 0.2; then, C = 5.31. A value for g 
obtained from conditions at the leading edge (stag- 
)\(£) behaves 


can be 


nation point here) by noting that (m/w. 
Consequently, the results of the 
From Fig. 1, for h = !/3 
and this is the value 


like (mé) as & > O. 
frst example can be used. 
and C = 5.31, one finds g = 0.81, 
From Figs. | and 2, one also 


that will be used here. 
9 = O.86m' Here 


finds (k)g;-9 = 4.96m 
the value of m, which is equal to (m'/u..)¢— 
taken as 34.0 for the suction side and as 9.50 for the 
pressure side.* with a, = 1.95 and 
a: = 3.80, it is found that PF, = 0.1178 J, = 
3.80. These used for both the suction 
and pressure sides, since at § = 0 they are independent 
«°(£) can now be obtained from Eq. (23) by numer- 
Simpson’s Rule was used for the pres- 
p (£) [= Cg/x(é)] are 


‘and (¢) 
, has been 


From Eqs. (22), 
and 


values can be 


of m. 
ical integration. 
ent calculations, and the results for ¢ 


shown in Fig. 4 
CONCLUSIONS 


The main results of this analysis for a sweat-cooled 
surface with a uniform wall temperature in laminar 


* In cases where actual accurate data for (u,/u,, )(é) are given, 
9 can be evaluated by fitting a polynomial to several 


reference 10). It 


ul u. jt 
points near £ = 0 and pr roy om. , e.g, 
may be noted that, in applying Eq. (23), it is not only sufficient 


but actually convenient to have (“,/u,, )() given as a set of nu 


merical values. 
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FIG.3 TYPICAL VELOCITY DISTRIBUTION OF 


SUBSONIC_FLOW OVER FORWARD PART 
OF A TURBINE BLADE 


compressible flow with a Prandtl Number of unity can 


be summarized as follows: 


(1) For any given velocity distribution outside of 
the boundary layer and free-stream Mach Number, the 
normal injection velocity distribution of the coolant 
required to maintain a given uniform wall temperature 
can be readily calculated by means of the equations 
Eqs. The Nusselt 
Number and the skin-friction coefficient can also be 
rapidly calculated by means of the equations derived 


be used for the 


(23) and (24)] developed here. 


here. The equations can, moreover, 


case of heat transfer with zero injection. 
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FIG.4 INJECTION PARAMETER DISTRIBUTION 


REQUIRED TO MAINTAIN A UNIFORM 
WALL-TEMPERATURE RATIO h OF '/3 
FOR EXAMPLE OF FIG. 3 


(2) For flow over a flat plate, the required normal 
injection velocity distribution depends only on a single 
parameter (C), which is a function of the ratio of the 
difference between the equilibrium wall temperature and 
the actual wall temperature to the difference between 
the actual wall temperature and the coolant temper 
ature. For flow with an axial pressure gradient, how 
ever, the required injection velocity also depends on 
the ratio of actual wall temperature to adiabatic equilib 
rium temperature. 

(3) <A simple relation exists [Eq. (27)| between the 
local Nusselt Number distribution (Nu) and injection 
parameter distribution (¢) required for a uniform wall 
temperature. This relation, which is valid for any main- 
stream velocity distribution and any Mach Number, 
states that Nu is proportional to g, the proportionality 
factor being a constant dependent only on the coolant- 
temperature parameter C. The the 
Nusselt Number along the wall will thus, in general, be 


distribution of 


CAL SCIENCES OCTOBER, 1952 
exactly similar to the distribution of the required nor 
mal injection velocity. 

(4) For objects in a flow with a forward stagnatioy 
point, the required injection parameter at that point 
will be proportional to the square root of the maiy 
stream velocity gradient there. 

(5) As the value of the coolant-temperature param 
eter C is increased, the required injection velocity js jn 
creased while the Nusselt Number is decreased. Fo; 
a fixed value of C, decreasing the ratio (h) of wall 
temperature to equilibrium zero-heat-transfer, tem. 
perature will decrease the skin friction in a favorab 
pressure gradient but will have relatively little effect oy 
the Nusselt Number. 
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On the Lift and Circulation of Airfoils in Some 
Unsteady-Flow Problems 


Nelson H. Kemp 

Graduate School of Aeronautical Engineering, Cornell University 
Ithaca, N.Y 

June 13, 1952 

I’ CONNECTION WITH SOME INVESTIGATIONS of unsteady flow 
in compressor cascades, it has been found useful to apply cer 

tain results of von K4rm‘*n and Sears" ? and to extend these re 

sults in two directions—(1) to calculate the circulation of airfoils 

in sinusoidal motion and (2) to generalize Sears’? results for an 

described 


airfoil flying through a sinusoidal gust in a manner 


below. 

1) If an airfoil in two-dimensional incompressible flow is sub 
jected to steady periodic variations of angle of attack or of verti- 
cal velocity, the ‘‘quasi-steady”’ circulation can be written as 


I'y = Goe'” 


where Gp is a constant and »y is the circular frequency. The 
quasi-steady”’ circulation is that produced by the instantaneous 
motion of the airfoil relative to the air without regard to the ef 
fects of the wake of trailing vortices produced by the nonsteady 
motion and is calculated from ordinary thin-airfoil theory 
Using the results of von Karman and Sears! for this case, it is 
easy to show that the total circulation I'(¢) is given by 


r(t) = Toe “S(w) (1 


where 

S(w) = 1/tw|Ko(iw) + Ki(tw)] 
vce/2U,c and U being the airfoil 
Ky and kK, are modified 


with the reduced frequency w 
chord and forward speed, respectively 
Bessel functions of the second kind 

Eq. (1) shows that S(w) bears somewhat the same relation to 
the circulation as the well-known Theodorsen function bears to 
the lift. (It must, of course, be remembered that the Kutta 
Joukowski relation between lift and circulation does not hold in 
unsteady motion. ) 

Now S(w) has never been calculated, although it is useful for 
determination of the energy in the wake. In fact, it is easily 
shown that, if x is measured rearward from the center of the air 
foil, then the jump in the velocity potential @ across the wake is 
given by 


— bottom = Toe 


ae ivx/l Sle) 
However, Professor Sears has called my attention to an interesting 


reciprocal relationship; he showed* that the /ift on an airfoil 
entering a sharp-edged gust is proportional to the circulation in 
that is, the total circulation due to a sud 


Now it is well 


the ‘Wagner problem” 
den unit change in the quasi-steady circulation. 
known that the steady-state response of any linear system to an 
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input of the form e'” is found from the Laplace transform of its 
response to a unit step-function input by replacing the variable 
of the Laplace transform by iy, dividing by vy, and multiplying 
by the input function. Applying this rule to the above-men 
tioned reciprocal relationship, we immediately see that the /ift 
due to an infinitely long sinusoidal gust must be proportional to 
the circulation of an airfoil in sinusoidal motion and, thus, pro- 
portional to S(w). 

The lift on an airfoil in such a sinusoidal gust was calculated 


by Sears? and found to be 


it | JA@)AKi(tw) + tJi(w)Ko(tw 
Lit) rpc l We" : (2) 
K (tw) + Kol tw) 
where |W is the amplitude of the gust velocity, J») and /; are the 
usual Bessel functions, and p is the air density 
The following identities regarding Bessel functions can be 
deduced :* 
Koliw) = —(r/2)i|Jolw) — 1Vo(w) 
A (iw —(2r/2 | Ji(w) iV (w 
Thus, 
Jil@o)AK (tw + 1 J,(w) Aol tw) (w/2)1 Jil@a)Vi((w) — JSi(w) Volo) 
But * 
J i(w)Vuailo) — Junsilw) (ew) —(2/rw 
2) is seen to be 


so that, finally, the expression in brackets in Eq. (2 


just S(w), and 


rpc UWe' S(w) (3) 


L(t) = 


is the result we expected. 

Sears? has plotted S(w) in his Fig 
S(w) = Siw) + iS;(w), is given in Table 1 of the present note. 
2) The result given in Eq. (3) is the lift due to a gust velocity 


3, and S,(w) and S;(w), where 


of the form 


ox, t) = Welrl! B/E] 


In connection with the effect of the wakes of stator blades on a 
rotor blade moving through them, it is necessary to calculate the 
lift due to a gust of the form 


V(x, t) = WeiMe—tua/t 


with »y ~ uw. This is easily done by a generalization of the de- 
velopment given by Sears,? and the result is 


L(t) = rpcUWe'" K(w, d) (4) 
where A pc/2L , With 
K(w, ») [Jol A) — TSA) IC(@) + (w/A)VSTi(A (5) 


is the Theodorsen function, defined as 


K,(iw)/|Koliw) + Ky(iw)| 


Here, C(w 


C(w) = 
and 77 of Reference 4 


* See pages 78, 75, 73 
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TABLE | 
i 
w S, S, w Ss, s 

0 1.000 0 2.5 0.048 0.247 
0.02 0.963 0.075 3.( 0.145 0.178 
0.04 0.924 0.114 3.5 0.197 0.081 
0.06 0. 887 0.139 1.0 0.198 0.021 
0.08 0.852 —0.154 1.5 0.155 0.106 
0.10 0.821 —0.164 5.0 0.081 0.159 
0.20 0.702 0.160 6.0 +0.081 0.141 
0.30 0.624 -0.126 7.0 0.150 0.008 
0.40 0.568 0.085 8.0 0.082 +0.114 
0.50 0.525 —0.044 9.0 —0.049 0.124 
0.60 0.488 0.005 10.0 0.124 0.002 
0.70 0.456 +0.032 11.0 -~0.084 0.086 
0.80 0.426 0.066 12.0 +0.026 —0.112 
0.90 0.397 0.097 13.0 0.104 -0.037 
1.0 0.369 0.126 14.0 0.084 +0. 065 
1.2 0.312 0.176 15.0 0.009 0.1038 
1.4 0.255 0.215 x 0 0 

1.6 0.196 0.242 

1.8 0.138 0. 260 

2.0 0.982 0.268 


While v is a real number, it should be noted that uw (and there 
fore \) can be an arbitrary complex number. é 

It is readily seen that K(w, w) = S(w), so that Eq. (4) reduces 
to Eq. (3) if uw = v. If, on the other hand, yp = 0, A(w, 0) = 
C(w) + (1/2)iw, and the lift from Eq. (4) agrees with that given 
in Eq. (20) of Reference 1 (where W = AoW, ¢ = 2). Thus, 
K(w, \) is seen to be a generalized function, which connects the 
familiar Theodorsen function for oscillating airfoils with the 
Sears function for a sinusoidal gust 

The circulation in this case can be easily found using Eqs 


(4), and (5), and(1}. The result is 


l= rcWei"e IT o(A) — iJ\(A)]S(@ 


The application of these results to the unsteady flow through a 


compressor stage will be published in the near future 
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Note on Creep Buckling of Columns 


George Gerard 

College of Engineering, New York University, New York 

June 30, 19592 

T' APPEARS FROM LIBROVE’S INTERESTING ANALYSIS! that, for 
the case of creep buckling of columns, the initial imperfec 

tions contained in ordinary columns provide the mechanism by 

In fact, 

it can be concluded from this analysis that a theoretically perfect 


which failure due to creep occurs after a period of time 


column that is initially loaded below the time-independent critical 
load will not buckle at all. 
case of static buckling where small initial imperfections play an 


his is an interesting contrast to the 


insignificant role, since the failing load of an initially imperfect 
column is substantially the same as that of a theoretically perfect 
column. It is of interest, therefore, to conjecture whether there is 
any possible mechanism by which a column containing no initial 
imperfections can fail as a result of creep when the initial load is 
less than the theoretical buckling load 
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First, it is pertinent to consider a column of a viscoelastic m 
ild 

terial. Such a model demonstrates rather vividly the significan, 
ms < t 


of initial imperfections in creep buckling. 
the application of an infinitesimal lateral disturbance to deter 


For a perfect column 


mine the point where an exchange of equilibrium configurations 
occurs merely results in a temporary increase in stress on the con 
cave side and decrease on the convex side of the column jf the 
application and removal of the disturbance is instantaneous 
Since the path of disturbance follows a line of slope E, the in. 
stantaneous bending rigidity is the same for all values of E in th, 
time-dependent phase, and, consequently, creep buckling cannoy 
occur. 

For the imperfect column, a small amount of bending is presen 
as a result of the initial loading. Consequently, in the tiny 
dependent stage of deformation a stress reversal occurs, with 
an increase in lateral deflection of the column. After some tim 
this deflection may become excessive, essentially constituting 
failure in a restricted sense. 

Instead of considering a viscoelastic material, the creep law 
used by Libove appears to be more realistic for aircraft materials 

€ = (o/E) + Ae**t! 

If a mechanical equation of state concept is applied to the creer 
law used by Libove, Eq. (1), then a series of ‘‘stress-strain”’ curves 
If the 
lateral disturbance is applied instantaneously, then the concave 


under conditions of constant strain rate can be found 


side would follow a loading path of slope £, and therefore the 
local bending stiffness would always have a value of £, thus pre- 
cluding creep buckling for the theoretically perfect column 
However, if it is now assumed that a lateral disturbance is applied 
at a strain rate corresponding to the local strain rate of the stress- 
strain curve, then the loading path on the concave side would 
follow the local tangent modulus and the bending stiffness would 
then be a function of the local tangent modulus. Since, with this 
assumption, the local bending stiffness would decrease with time, 
creep buckling would be possible after a ‘‘critical time.” 

The assumption that the lateral disturbance is applied at the 
local strain rate is based on the room temperature behavior as 
sociated with plastic buckling of a column. For such a column 
there is an implicit assumption that the loading on the convex 
side follows the stress-strain curve pertinent to the strain rate at 
which the axial loading is applied. If this assumption were not 
implied but, instead, an instantaneous lateral disturbance wer 
assumed, then at room temperature, as well as at elevated tem 
perature, the local bending stiffness would have a value of E and 
buckling would not occur. Thus the assumption concerning the 
strain rate of the disturbance appears to have a somewhat realistic 
basis 

With use of this assumption, the solution for creep buckling of 
a theoretically perfect H-section column reduces to determining 
the “critical time’’ at which the following equilibrium equation-is 
satisfied 


scAw = 0 y 


E,(t)I(d?w/dx?) 4 


In this case @ is the applied stress and is constant, and the lateral 
disturbance produces a deflection w essentially independent of 
time ¢. Hence, when £,(t) has reduced sufficiently to satisfy Eq 
(2), the critical time is reached. The solution of Eq. (2) is the 
usual column equation 
PE (t)/(L’ /p? 3 


Since the average stress & is held constant, the buckling process 
must occur with a strain reversal, and therefore the correct value 
With this 


suggested mechanism, a finite critical time is predicted for creep 


of the effective modulus is the reduced modulus, EF; 
buckling of a theoretically perfect column 


REFERENCE 
Creep Buckling of Column Journal of the Aeronautical 
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Transonic Flow Past a Wedge with Detached 
Shock Wave 


AH. Armstrong, M. Holt, and C. K. Thornhill 
Armament Research Establishment, Fort Halstead, Kent, England 


R centi and Wagoner? 


transonic flow past a finite wedge (Fig. 1), using Tricomi’s equa 
tion. In both cases treated (AZ = 1 and M > 1, respectively), 
consideration is restricted to wedges such that flow downstream 
of the shoulder does not influence the subsonic region of flow up 
tream of the sonic line CF so that a Tricomi-Frankl* boundary 
In this note the condition is 


ECENT PAPERS BY GUDERLEY AND YOSHIHARA! and by Vin- 


present approximate treatments of 


yalue problem can be formulated 
determined for this to be true 

The condition clearly depends on the extent of the expansion 
Provided @ < (1/2) (@ + a), then, in the hodo 
graph plane (Fig. 2), the backward characteristic CD defining the 
centered wave at the shoulder, will extend far enough to meet 
the forward characteristic through F at some point G; FG will, 
in fact, then be one of the forward characteristics of the centered 


A Tricomi-Frankl boundary-value problem can then be 


at the shoulder. 


wave 
formulated in the hodograph plane, bounded by FABCG 

If, however, (1/2) (@ + a), the forward characteristic 
begins on the wedge surface downstream of the 


through F 
the centered 


and the last forward characteristic of 
waveis DE. In this case the boundary FABCD is rot complete, 
and, to render the problem determinate, the conditions on the 


These are evidently 


shoulder 


portion of the sonic line EF must be known. 
determined by the wedge contour downstream of the shoulder. 








Physical plane 


Fic. 1 


FORUM 715 











Hodograph plane 


Fic. 2 


Acknowledgment is made to the Chief Scientist, Ministry of 


Supply, for permission to publish this note 
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Remarks on ‘‘An Analytical Approach to the 
Fuel Sloshing and Buffeting Problems of 


Aircraft’’ 


M. Z. Krzywoblocki 

Professor of Gas Dynamics and Theoretical Aerodynamics 
Department of Aeronautical Engineering, University of Illinois 
Urbana, III 

June 15, 1952 


wrote, on page 


TT AUTHORS OF THE PAPER (see reference | 
227: “Some new ideas on buffeting were recently reported 
And next: ‘The central points in Liepmann’s 

suggestions are that (1) the wake of a stalled body is turbulent 
ind is described mathematically by The 
writer takes the liberty to mention that the idea of a turbulent 


by Liepmann.’ 


mean values 


wake is probably as old as the concept of Bénard-Karman vortex 
street. Previous investigations were concerned mainly with the 
neasurements of characteristic dimensions and frequencies of a 
Bénard-K/ rman vortex street, as well as with the stability of its 
configuration. But here and there one may trace some indications 
that Bénard-Kérman vortex street is in some relation to a turbu 
lent wake In the Fage and Johansen experiments,’ the varia 
tions of the spacing ratio in the wake was thoroughly investi 
In a Bénard-Karman vortex range, the fluctuations have 
Further downstream, the existing vorticity 
3 and 4) “‘in 


gated 
definite frequencies 
appears as “isotropic turbulence’ (see references 
which discreet vortex frequencies cannot be detected but a more 
or less continuous spectrum of eddy frequencies of much smaller 
scale and intensity than the Bénard-Karman vortexes exists.’’ 
An analysis by Hooker’ attributes the downstream variation in 
the spacing ratio principally to the viscosity. The change in the 
configuration of the wake in the vicinity of the critical Reynolds 
Number from a distinct periodic motion to a motion in which defi 
nite periods are lacking (i.e., random motion, remark of the 
writer of the present note) was discussed by Goldstein (.Wodern 
Developments in Fluid Dynamics, 1938) and others.” But for a 
certain length of the wake with no distinct Bénard-Karman vor- 
tex street, the frequencies of the most prominent disturbances in 
the wake could be recorded Krzywoblocki,* following the opin 
ion that to obtain a better review of this phenomenon one has to 
look for a solution in the theory of turbulence, applied certain 
results of Schlichting’s calculations of the flow in the boundary 
layer along a flat plate (1935) and obtained the values of the 


Strouhal Number in the turbulent wake which agree with the 








716 JOURNAL OF THE 


results cited previously by Goldstein (Vol. II, p. 419, Fig. 149). 
The idea of Professor Liepmann to apply some results of the the- 


ory of isotropic turbulence to buffeting calculations seems to be 


interesting, but for the time being, at least, it is not known ex 
actly how far behind a body the turbulence in the wake is fully 
isotropic. It seems that the transition from a distinct Bénard 
Karman vortex street to an isotropic turbulence is a gradual one 
This would mean that in the turbulent part of the wake between 
the Bénard-Kaérman vortex street and fully developed isotropic 
turbulence there are some recordable frequencies of the most 
prominent disturbances signalized by Goldstein and others long 
ago. It may be that the wake in its entire length is a Kolmo 
gorofi’s system of locally isotropic turbulence 

The recent investigations and hypotheses of the structure of a 
turbulent wake with no distinct vortices by Townsend, Kovasz 
nay, and Batchelor seem to be of great interest to both theore 


ticians and designers. 
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Supersonic Nozzles with Continuous Wall 
Curvature 


R. C. Frost and H. P. Liepman 

Aerodynamic Engineer, Consolidated Vultee Aircraft Corporation, 
Fort Worth, Tex., and Director, Supersonic Wind Tunnel 
Engineering Research Institute, University of Michigan, 
respectively 

June 15, 1952 

fluid 


Crown DEVELOPED A METHOD for designing a_ perfect 


nozzle with continuous wall curvature, since the elastic 
contour of a flexible wall nozzle will admit only this type of cur 
vature. His analysis further assumes that the point of inflection 
at the wall lies on the ‘‘patching’’ Mach line, which divides the 
region of curved intersecting characteristics from the simple wave 
region with straight nonintersecting characteristics 

The requirements of continuous wall curvature can also be 
shown to exist on the basis of a physically desirable Mach Num 
ber distribution along the centerline of the nozzle. In this case, 
however, the point of inflection must be upstream of the end 
point of the patching Mach line 

The existence of this requirement was discovered in 1950 by 
Frost, then a Research Assistant at the Engineering Research In 
stitute, during the design of a Mach 2.5 nozzle for the 8- by 13-in 
supersonic wind tunnel of the University of Michigan In this 
design the Nilson method* was modified to eliminate discontinui 
ties in the first derivative of the assumed Mach Number distribu 
tion along the nozzle centerline Evvard and Marcus, of the 
N.A.C.A., have recently given the same relations between con 
tinuous wall’ curvature, axial Mach Number gradient, and loca 
tion of inflection point with respect to patching Mach line 

Since the discussion of reference 2 has not been published and 
is based on physical reasoning instead of on equations as done by 
Evvard and Marcus, an abbreviated version may be of interest 
1 was calculated from an as 
A sinus 


The nozzle contour shown in Fig 
sumed Mach Number distribution along the centerline 
chosen to assure continuity of d.\//dx 


a final value of d\//dx = 


oidal distribution was 


along the nozzle with dat point P’ 
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(Fig. 1). It was reasoned that, if the nozzle were designed for 
continuous changes in flow, its actual performance would be least 
affected by manufacturing errors and variations in stagnation 
conditions 

Referring to Fig. 1, the flow in the nozzle can be discussed jn 
terms of the strength of the two families of characteristics, th, 
strength of a characteristic being proportional to the change of , x 
pansion angle or Mach Number. Since d\//dx is made to g0 to 
zero at point P?’ continuously, the strength of the reflected or 
right running family of characteristics R must do the same. [py 
the initial portion of the nozzle starting at point 7”, dM/dx js 
nearly constant, and thus the strength of the incident or left 
running characteristics 1 must be nearly constant there and syh 


sequently decrease continuously to zero also at point 7?’ 








l" 
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Design 





¥ 


Figure |. Typical Characteristics and Mach 
Number Distribution on Nozzle Center Line. 


Along the wall, between points 7 and P, the change in slope 
I g | 


of the contour—i.e., the curvature—is proportional to the differ 
ence in strength between reflected and incident characteristics 
since the flow must be tangent to the wall. To provide the re 
quired expansion, the slope of the contour increases downstrean 
of 7 and, hence, the strength of family R must be greater thar 
that of family L. Since the strength of family FR is continuously 
decreasing, it will be equal to the strength of family L at som 
point /, the point of inflection, where the curvature is zero 
Downstream of this point /, the strength of family R continues 
to decrease smoothly until it is zero at the patching Mach line 
point P? 

Downstream of the patching Mach line, the contour between 
points ? and E is entirely determined by the continuously dé 
creasing strength of family L which vanishes at £, thus ensuring 
continuous curvature at E also 

Thus, since all changes in slope of the wall contour betwee! 
points 7 and E have been accomplished continuously, the curva 
ture is continuous along the wall, and it has been shown that the 
point of inflection must be upstream of the patching Mach line 
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Direct Measurement of Supersonic Skin 
Friction 


Donald Coles 
Research Engineer, Jet Propulsion Laboratory, California Institute 
of Technology, Pasadena, Calif 


arab EXPERIMENTS AT THE Jet Propulsion Laboratory have 
established the effect of compressibility on local skin friction 
for Mach Numbers up to 4.5 Friction coefficients were obtained 
imultaneously at three stations on a flat plate, for both laminar 
and turbulent boundary-layer flow, for Reynolds Numbers from 
+X 10° to 1 X 10%. The floating element technique used is 
jmilar to that of Liepmann and Dhawan.! 

Figs. 1 and 2 show typical measurements with a clean leading 
edge and with a tripping device. The Reynolds Numbers in the 
figures are based on the distance from the leading edge and on the 
free-stream velocity and kinematic viscosity. No corrections for 
local pressure gradient or Mach Number variation have been 
made 

Fig. 3 shows the dependence of turbulent local friction on Mach 
Number at a Reynolds Number of 8 X 10°. It is believed that 
these data will remain substantially unchanged by any corrections 
to be applied during the analysis of the experiments 
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Fic. 3. Variation of turbulent local skin friction coefticient with 
Mach Number at a Reynolds Number of 8,000,000 


The boundary-layer research summarized here was discussed 
in detail at the Annual Summer Meeting of the Institute of the 
Aeronautical Sciences in Los Angeles in July, 1952 
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On Theodorsen Function in Incompressible 
Flow and C-Function in Supersonic Flow 


Chieh-Chien Chang 
Research Contract Director, The Johns Hopkins University, and 
Consultant, Air Research and Development Command, U.S.A.F 


June 6, 1952 


Freee INTERESTING PAPERS!~! in the Reader's Forum of the 

March JouRNAL discuss the validity of generalizing the Theo 
dorsen function to the stable airfoil oscillation. To readers not 
too familiar with flutter problems, it should be emphasized that 
these discussions are not on the validity of the original Theodor 
sen® function. The current confusion lies entirely in the extension 
of the Theodorsen function to stable oscillations 

In the notation of reference 1, k = ko + ik; (where &p is real, 
0, respec- 
tively for unstable, harmonic, and stable oscillations. A sketch 
shows airfoils in flight at constant velocity V, leaving wake vor 


positive and & is real). The references'~* take ;, 


tex sheets corresponding to the three cases, k; 5 0. It is as- 
sumed that all the airfoils started from rest at past time —s/V 
and produced wakes of equal finite length s. Observe that the 
maximum magnitude of vortex strength U/(s)* in the wake sheet 
exponentially decreases with s for ky < 0, is constant for ky = 0, 
The limit of case 
Actually, 


and exponentially increases with s for k; > 0. 
is the Theodorsen function c(k) 


x s $ + ] 
f etks as/f ; 
1 W 52 1 j 1 Vs52 1 


* The author is fully aware that for a finite wake the vortex strength also 


(b) ass — 


c(k) = tS ds 


depends on the time U(s) is purely for demonstration 
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which is neither finite nor unique in value unless we first impos 
k; < Oand then let k} ~ 0. From the sketch for k; 
Fig. 1), it is obvious that the maximum magnitude of ls) goes 
From the physical point of 


> 0 (bottom, 


to infinity exponentially as s > @ 
view, it is not surprising that c(k) is divergent for k; 
out,'~* because the vortex strengths at the extremely large values 


> 0 as pointed 


of s contribute more to the value of the integral. 


Furthermore, the Theodorsen function generalized to the 
stable airfoil oscillation has no physical reality. 


of flutter analysis is to determine the critical frequency at which 


The primary aim 
airfoil oscillations occur. The theory was developed for vanish 
ingly small oscillations. To apply the generalized Theodorsen 
function to calculate the aerodynamic behavior of an airfoil, a 
few limitations must be considered. 

(a) The amplitude of the airfoil oscillation present or past as 
shown in the wake vortex strength must be small in order to 
satisfy the original assumption of Theodorsen 

(b) In incompressible flow, time derivatives do not enter into 
the differential equation. However, the transient effect appears 
in the vortex and source strengths determined from the boundary 
conditions. Judged in this sense, the vorticity shed by the os 
cillating airfoil during the whole past history of flight will influ 
ence the present flow. Since the viscosity in the fluid takes time 
to dissipate the kinetic energy into heat energy, the vortex 
strength in the wake must deteriorate progressively with the past 
time interval or s. Even if exponentially increasing vorticity 
with s could be produced, the viscosity will dissipate all the vor 
ticity at © into heat in the infinite time interval 

It is impossible physically that vorticity produced a long time 
ago still maintains its strength and influences the flow now. The 
Theodorsen model of constant oscillation magnitude must de 
teriorate in wake vortex strength with s. This suggests that 
k,; < 0 may be assumed to take care of the viscosity dissipation 
It is worth while to investigate what 
This also 


in the Theodorsen model 
value of k; is equivalent to the viscosity dissipation 
offers a physical justification for the method of evaluating c(k) 

A finite length wake has physical meaning in all three cases 
k, = O, because, the magnitude of U(s) is always boundled If 
the time interval is very small, the dissipation may be unimpor 
Unfortunately, the Wagner function cannot be expressed 
T. Jones’ approximation for the Wagner 


tant. 
in a closed form. R. 
function should be applied cautiously and should be compared 
against exact solutions wherever possible 

In supersonic flow, the corresponding Wagner function has been 
found. By means of a convolution integral, the effects of the 
transient and quasi-steady harmonic oscillation have been found 


and expressed as incomplete and complete C-functions, respec 








a(t) 
Yo 
Fic. 1. (Yop) Unstable oscillating airfoil, k, <0. (Center) 


Harmonically oscillating airfoil, k, = 0 (Theodorsen’s case, s > 
). (Bottom) Stable oscillating airfoil, k; >0 
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takes 


It corresponds to van der 


ively. [he incomplete C-function IntO account the 
transient effect of abrupt oscillation 
Vooren’s results of finite wake in incompressible flow 
Miles 


function 


Recently 
identified the complete C-function with th« Schwart, 
Fortunately, the supersonic wake has no influence or 
the aerodynamic behavior of the airfoil, and, consequently, the 
viscosity dissipation in the wake negligibly influences the airfoi 
pressure. Thus the aerodynamic load is valid for any oscillatio, 
stable (k; 
0 has been investigated,® the corresponding functions for }, + | 


> 0) or unstable (Rk; 0). Although only the case }, = 


can be deduced easily with the same principles 


REFERENCES 

1 van der Vooren, A. [., Generalization of the Theodorsen Function to Siy} 

Oscillations, Readers’ Forum, Journal of the Aeronautical Sciences, Vol. | 
No. 3, pp. 209-211, March, 1952 

2 Laitone, E. V., Theodorsen's Circulation Function for Generalized M Li 


Readers’ Forum, Journal of the Aeronautical Sciences, Vol. 19. No 2 pI 
211-213, March, 1952 

Jones, W. P., The Generalized Theodorsen Function, Readers’ Forum 
Journal of the Aeronautical Sciences, Vol. 19, No. 3, p. 213, March, 1952 


and Luke, Y. I Votes on the Calculation 


Readers’ Forum 


Dengler, M. A., Goland, M 
Stable Ae 


Aeronautical Sciences, Vol. 19, No. 3, pp. 213-214 


Journal of the 
March, 1952 


the Response of odynamic System 











Theodorsen, Th General Theory of Aerodynamic I) 
Vechanism of Flutter, N.A.C.A. T.R. No. 496, 1935 

Chang, Chieh-Chien, 7vansient Aerodynamic Behavio an Airfotl Du 
to Different Arbitrary Modes cf Non-stationary Motion in a Supersonic I 
N.A.C.A. T.N. No. 2: April, 1951 

Miles, John W., On Chang's Function for Nonstationary Flow, Readers 
Forum, Journal of the Aeronautical Sciences, Vol. 19, No. 2, p. 138, Febrn 


iry, 1952 


Radii of Gyration of Aircraft 


J. P. Chawla 

Aerodynamics Department, Hughes Aircraft Company, Culver City 
Calif. 

June 9, 1959 


Momests OF INERTIA OF AIRPLANES are needed in the prelim! 
nary design stages of an aircraft in order to estimate dy 
namic loads on various components of its primary structure and 
to evaluate its dynamic stability characteristics. Various 
authors have put forth simple formulas for the radii of gyration 
of aircraft, about the body (or principal) axes through its center 
of gravity, expressed in terms of its overall dimensions of height, 
length, and span. These studies have invariably been from 
statistical point of view, using data on airplanes built and success 
fully flown, and they have proved valuable to the preliminary 
design engineer. This note is intended to report significant re 
sults of an analysis recently made at Hughes Aircraft Company 
in which the same basic approach was used 
Data on 101 airplanes, 73 of American design and 28 of British 


design, were analyzed. The airplanes were subdivided into four 


namely, (1) fighter, (2) bomber, (2 
Thirty-three of the air 


planes analyzed were fighters, 22 bombers, 26 cargo transports 


main functional categories 
cargo transport, and (4) flying boat 
The fighters range in design gross weight 
the bombers, from 6,000 to 400,00 
and the 


ind 20 flying boats 
from 5,000 to 40,000 Ibs.; 
Ilbs.; the cargo transports, from 24,000 to 300,000 Ibs 
flying boats, from 8,000 to 400,000 Ibs 

Radii of gyration about the three axes were nondimensionalized 


in the following manner 

fg = kz/b, ry = ky/L, 1, = he/e 
where k;, ky, and k, are the radii of gyration in roll, pitch, and 
yaw, respectively; 6 is the span, L is the overall length, and ¢ 1s 


the arithmetic mean of 6 and L—that is, 


=(/)(b+ L 
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The criterion used for obtaining a mean value for these nondimen 
jonal ratios was to find a number, ./, which would give the 
greatest number of points in the band between 0.95.7 and 1.05. 
The results indicated that only unconventional designs may b« 
expected to vield more than 10 per cent deviations from the mean 

Mean values of rz, 7,, and rz for each category are listed in 
Table 1 Mean values of ry, and r, for flying boats are not in 
cluded because of paucity of data. We note that values of 7 
for different functional categories of aircraft lie in a narrow 


ind rz : 
A single 


band, while the values of r, extend over a wide range 
mean value of r, = 0.165 confirms a previous study? in which the 
airplanes were not subdivided into different categories. 


TABLE |! 
Mean Values of Nondimensional Radius of Gyration Ratios 


rype of Aircraft le , Ts 

Fighter 0.120 0.165 0.200 

Bomber 0.155 0.160 0.208 

Cargo Transport 0.145 0.165 0.211 
0.165 


Flying Boat 


A detailed analysis of the radii of gyration in roll showed that, 
with the mean values of r, chosen for each category, 35 per cent 
of the airplanes analyzed gave a deviation of less than 5 per cent 
ind 45 per cent gave a deviation of less than 10 per cent. This 
bears out the fact that the radius of gyration in roll varies consid 
erably with the design of the airplane—in particular, the type 
ind spanwise location of power plant, wing taper (both in plan 
ind elevation), and the presence of tip pods, if any 

The selected mean values of r, put 48 per cent of the airplanes 
in the 0 to +5 per cent band of deviation and 74 per cent in the 0 
to +10 per cent band. It may be inferred that the value of (2, = 

is fairly constant 

From analysis of the radii of gyration in yaw, it was found that 
56 per cent of the airplanes had less than 5 per cent deviation and 
83 per cent had less than 10 per cent deviation from the mean 
value. Constancy is again manifest in the value of the parameter 
k./e 

In closing, it may be remarked that statistical studies have to 
be revised from time to time as design trends undergo significant 


hanges 
REFERENCES 
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paper presented to the Society of Aeronautical Weight Engineers, Inc., at 
their Eleventh National Conference in Buffalo, N.Y., from May 5 to 8 
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On Folds Occurring in the Mapping from the 
Physical Plane to the Hodograph Piane 


Keith C. Harder and E. B. Klunker 
Langley Aeronautical Laboratory, N.A.C.A., Langley Field, Va 
June 20, 1959 


N THE QUALITATIVE study of two-dimensional transonic flow 


problems, it is often enlightening to consider sketches of the 
flow in the physical and hodograph planes simultaneously. Such 
considerations yield considerably more information than studies 
restricted to The starting point of such 


either plane alone 


studies is usually the physical plane. One troublesome point in 
sketching the corresponding hodograph is locating the folds in 
the hodograph plan Che following rule has proved useful in 


this regard 


S’ 
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RKule._-A line in the physical plane, separating expansion and 
compression waves of the same family, corresponds to a fold in 
the hodograph plane. 

Proof._-A characteristic A BC (see Fig 
il? ps into a unique line in the hodograph plane 
characteristic from A to C, at first only expansion waves are 
crossed, and the velocity increases until point B is reached. From 


1) in the physica! plane 
Traversing this 


point B to C, the velocity decreases, since only compression waves 
are crossed. Since we must stay on the same characteristic in 
the hodograph plane, part of this characteristic is covered twice 
in the hodograph plane. The locus of points in the hodograph 
plane (such as 8), where we first start to retrace our path, is 
termed a fold. Such a retracing always occurs as we move from 
ul expansion field to a compression field (or vice versa) 


Compression characteristics 


— —— Expansion characteristics 


c,c; 


Hodograph plane 


Physical plane 


he rule gives the following information concerning the map 
ping of flows with shock waves: 

(a) Flows with a shock wave in an expansion field will have one 
fold along the initial conditions of the shock and one along the 
final conditions. 

(b) Flows with an expansion field immediately upstream of the 
shock followed by a compression field immediately downstream 
of the shock will have a fold along the initial conditions of the 
shock 

(c) Flows with a compression field immediately upstream of 
the shock followed by an expansion field immediately downstream 
of the shock will have a fold along the final conditions of the 
shock 


Determination of Modes and Frequencies 
Above the Fundamental by Matrix Iteration 


P. Cicala 
Professor of the Politecnico di Torino. Now at the E.S.A. and 
University of Cordoba, Cordoba, Argentina 


July 2, 19592 


Moats ITERATION PROCEDURES tay be used to calculate the 

second mode and frequency if the dynamical matrix is 
postmultiplied by a factor eliminating the first mode component 
from input deformations. A general expression for this factor 
is here indicated, from which the sweeping matrix of references 
1 and 2 is deduced as a particular case. It is shown that, in 


general, inaccuracies in the first modes markedly influence sub 


~ 
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sequent results, the errors being magnified by factors containing 
the square ratios of higher to lower frequencies. By a particular 
matrix factor, sweeping out the first mode component from the 
rows of the dynamical matrix, the effects of first mode errors are 
considerably reduced. 

These concepts lead to general formulas for the calculation of 
the rth mode which reduce the influence of the errors affecting 
the previously calculated r — 1 lower frequency modes to second- 
order effects. 

The dynamical equations for the free vibrations of a system 
with » degrees of freedom may be written in matrix form 


Dig} = Ati (1) 


where D is the dynamical matrix, product of the flexibility in- 
fluence matrix by the inertia matrix VW; {gq} is the column of 
parameter amplitudes; and A = 1/w? (w = natural frequency) 
Let Q be the matrix of modal columns whose jth column {q;} is 
formed by the parameter values representing the jth mode, ap 
propriate to the eigenvalue 4 = \,;. Let P be the transpose of 
MQ; its ith row [p;]| is the characteristic row appropriate to 
Ai. Let the modal values be normalized so that PQ = I, where 
I is the unit matrix of order n. 

It can be shown that the matrix ) may be expanded as fol- 


lows :* 


D= > Av tges [Pe | 
k=1 


Any system of » values, in particular the column jg}, may be 


expanded in the series of normal modes; thus we write 


tq} = Q{A} (3) 
The 7th element A; of the last column is the coetticient of the 
ith term in this modal expansion. 

Iterative solution of Eq. (1) yields approximate values for the 
first mode which we denote by symbols with superposed bars; 
these are expressed as follows: 

1M} = VU for+ er; (4) 

lpr] = [6 + ¢,|P (5) 
with 6; = 1 for the first element and 6; = O elsewhere; the (small) 
coefficients in the column {e,! and its transposed |[e;] represent 
in modal expansions the errors by which the calculated mode is 
affected. We denote by e the square root of the quantity 
fer} fe}. 

According to Eqs. (3) and (5) we get 


laliq) = 41+ 
le |{A{. We put 


[— 


where s; = 

10} [pr | (7) 

where {)} is a column of arbitrary coefficients which may be ex 

pressed by 

{b} = Q{B} (8) 

in terms of the new coefficients B), > Be Thus, from Eqs. 
(3), (6), (7), and (8), we obtain 

S{q} = Q{A‘} (9) 

The ith element of the column |A’} is given by 

A; — B;A,; — Bys; (10) 


1 , 


* This equation shows an interesting analogy with the expansion for the 


kernel of the corresponding integral equation 
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If error terms are neglected and the following condition js satis. 
tied 


Bi = [piltbj = 1 (11) 


then A,’ = U—i.e., premultiplication by S sweeps out from the 
column j{q} the first component in the modal expansion, There. 
fore, if the matrix D in Eq. (1) is replaced by the product DS, 
the new equation gives, through iteration, the second mode, 
whatever the values of Bo, , By are. 

The above equations yield, for S, the expression of the sweeping 
matrix indicated in references 1 and 2 if we give zero values to 
all the coefficients 6 but one, which, according to Eq. (11), must 
be equated to the reciprocal of the corresponding , factor. Ip 
this case, B; is represented by the ratio p;/p; at this same station, 

In order to evidence the effects of errors in the first calculated 
mode, the asymptotic relations between the coefficients A, after 
infinite iteration cycles, are considered, neglecting terms of the 
order of «7. Denoting by p2 the ratio of the values of Az in two 
successive steps of iteration, from Eqs. (2), (9), and (10), with 
B, = 1, we get 


2 L + é€:2Bo(Ar — Az2)/A (12) 


p2/A2 = 


where ¢j2 is the second element in the row [e,]. The error term 
is considerably magnified by the factor (A: — A2)/As; it affects 
the frequency calculation. This and similar errors appearing in 
the mode shape would be nullified if we could make B, = 

.= B, = 0. The best possible approach to this ideal condi- 


tion is attained by taking {>} = {G,} and, consequently, 


DS = D — vyf{Hn} lpr) (13) 


having replaced \,{4:} for D{g,}. Thus the error in Eq, (12) re- 
duces to the order of «,?. 

In more general terms, for the calculation of the rth mode, after 
the preceding r — 1 have been determined with certain error ¢, 
we replace the matrix D by 

1 
D® =D— YY rihduh pel (14) 
k=1 

With this matrix, the iteration process gives in the limit an 
rth mode {g,} which contains errors of the order of e* in the coef- 
ficient A,. Consequently, from it the rth eigenvalue can be cal- 
culated exactly up to e® terms. The calculated mode contains 
impurities of the order of €; these can be eliminated by evaluat- 
ing {gr} as follows:t 


Gr} = tGr5 — ,> 


If main interest is placed upon certain overtones, for other fre- 


[Pelhdrt (es (15) 
1 Ak — 


quencies the iteration process can be less accurate and the cor- 
rection by Eq. (15) can be omitted, since only squares of mode 
errors enter in subsequent calculations. 


REFERENCES 
! Flomenhoft, H. 1., A Method for Determining Mode Shapes and Frequew 
cies Above the Fundamental by Matrix Iteration, Journal of Applied Mechan- 
ies, Vol. 17, No. 3, pp. 249-256, September, 1950 
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+ If the modes are not normalized, the kth term of the summations @ 
Eqs. (14) and (15) must be divided by [peliq, } 
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